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1.  LITERATURE  SURVEY 


1.1  Introduction 

The  most  widespread  commercial  use  of  the  laser  is  in  cutting, 
welding,  drilling,  and  heat  treatment.  Many  applications  involve  lasers 
of  sufficient  power  (Intensity)  such  that  a  laser  Induced  plasma  can 
form.  The  plasma  is  defined  as  the  collection  of  vaporized  material  and 
ambient  species  which  are  excited  to  a  high  enough  energy  state  such 
that  the  gases  can  interact,  via  absorption,  with  the  laser  energy.  A 
review  of  the  literature  will  demonstrate  that  the  plasma  effects  are 
generally  circumvented  by  mechanical  means.  With  the  advent  of  high 
peak  power  pulsed  lasers  (In  excess  of  10 J  per  pulse)  and  CW  lasers  In 
excess  of  IkW,  research  has  focused  on  the  quantitative  measurement  of 
laser-plasma- target  interaction.  The  above  efforts  have  been  primarily 
Interested  in  measurement  of  "enhanced  coupling"  phenomena  at  fluxes 
greater  than  uAf/cm2.  Enhanced  coupling  is  the  effect  where  more  laser 
energy  has  apparently  coupled  to  the  target  than  can  be  explained  via 
simple  radiation  absorptivity  and  thermal  conductivity. 

Along  with  the  experimental  work  several  models  of  plasma 
initiation  and  maintenance  were  developed.  As  laser  powers  increased 
most  of  the  studies  evolved  Into  the  areas  of  fusion  and  damage.  Most 
of  the  work  was  subsequently  classified  about  1982.  Some  of  the  studies 
were  concerned  with  laser  propulsion  which  received  a  renewed  interest 
In  the  early  1980's.  From  this  work,  a  general  understanding  of  plasma 
phenomena  can  be  developed.  The  goal  of  this  effort  is  to  parameterize 
laser-plasma-target  Interactions  In  the  sub-107W/cm^  Intensity  range  and 
utilize  the  results  to  accurately  predict  heat  affected  zones  (HAZ)  in 
targets. 


1.2  Material  Processing 

One  of  the  earliest  exoerimental  studies  of  laser  coupling  to 
targets  was  presented  by  Locke,  et  al.  [1).  The  authors  were  concerned 
with  the  non-uniform  penetration  depths  in  laser  welding.  Half  Inch 
thick  304  stainless  steel  plates  were  welded  at  CW  laser  powers  from  8 
to  20  kW,  f/8  optics,  1.00mm  spot  diameter.  The  results  indicated  that 
the  edge  fit  had  to  be  less  than  1/10  the  thickness  for  a  successful 
weld.  The  authors  were  interested  In  efficiency  of  the  laser  process 
and  defined  an  arbitrary  value  of  30*.  The  efficiency  was  defined  as 
the  ratio  of  the  actual  coupled  power  to  the  power  needed  for  a 
hypothetical  case  of  zero  conduction.  The  results  Indicated 
efficiencies  of  22  to  45  percent  for  laser  powers  of  3.8  and  20kW, 
respectively.  The  values  represented  a  slight  Improvement  over  E-beam 
welding  which  Increased  with  Increasing  power.  Finally,  a  qualitative 
study  of  gas  assists  was  made.  A  gas  assist  with  a  component  normal  to 
the  target  surface  resulted  in  a  deeper,  less  uniform  penetration.  A 
transverse  gas  assist  was  seen  to  develop  a  "more  luminous...  concen¬ 
trated  Ionization  cloud",  but  most  of  the  laser  energy  appeared  to  be 
transmitted  to  the  target.  The  transverse  gas  assist  resulted  In  a 
uniform  weld.  With  no  gas  assist,  the  plasma  was  reported  to  shield 
most  of  the  laser  energy. 

Plrrl,  et  al.  12)  Investigated  the  momentum  transfer  to  a  target 
when  Irradiated  by  a  pulsed  C02  laser,  peak  flux  of  3  to  6xlO®W/cm2, 
25  vs  pulse.  The  results  of  ballistic  pendulum  measurements  indicated 
that  momentum  decreases  with  Increasing  pulse  energy.  The  effect  was 
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Independent  of  material  for  carbon,  aluminum,  and  tungsten.  The  impulse 
also  increased  with  target  area.  This  result  led  to  the  hypothesis  that 
the  momentum  transfer  was  due  to  a  surface  expansion  wave  and  that  the 
plasma  must  therefore  be  radlatively  driven.  Pirri  assumed  all  of  the 
laser  energy  was  deposited  In  a  point  above  the  target,  forming  the 
center  of  a  blast  wave.  Blast  theory  was  then  applied  to  show  that  the 
maximum  momentum  transfer  occurs  if  the  target  is  large  enough  for  the 
blast  wave  to  equllllbrate  prior  to  expansion  off  the  target.  A  video 
system  was  used  to  monitor  the  expansion  velocity  at  a  few  microseconds 
into  the  laser  pulse.  The  expansion  velocity  was  approximately  60 
percent  below  the  predicted  blast  wave  velocity  of  1.8xl06cm/s.  By 
increasing  the  spot  size,  the  authors  bounded  a  transition  between 
combust Ion-type  and  detonation- type  waves  to  be  107  to  4xlO®W/cm2,  which 
Is  close  to  the  classical  detonation  value.  At  107W/cm2  most  of  the 
laser  energy  was  assumed  to  be  absorbed  Into  a  thin  surface  plasma,  but 
this  was  a  consequence  of  the  observed  plasma  shielding  at  higher 
fluxes.  A  detailed  treatment  of  combustion  and  detonation  phenomena  Is 
presented  In  Chap.  2. 

Lowder,  et  al.  {31  Investigated  target  Interaction  with  a  pulsed 

C02  laser.  The  pulse  shape  consisted  of  a  1  us  spike  followed  by  a 

lower  plateau  of  10-30  ys.  The  authors  were  one  of  the  first  to  admit 

that  the  experimental  spot  size  was  not  diffraction  limited.  Using 

streak  photography  the  expansion  velocity  was  measured  to  be  105  to  106 
cm/s  and  constantly  decreasing.  Implying  combustion  wave  phenomena.  The 

camera  results  Indicated  that  laser  absorption  occurred  throughout  the 

plasma  volume  and  the  blast  wave  calculations  accounted  for  the  long 
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absorption  length.  The  conclusion,  as  In  Plrrl's  work  [2],  was  that  the 
combustion  to  detonation  transition  occurred  at  an  Intensity  of  4xl07 
W/cm*  (±25  percent)  for  a  beam  radius  of  0.2  to  0.3  cm.  Using  an 
electrostatic  probe,  the  authors  found  that  the  air  Instantly 
photolonlzes  from  the  breakdown  radiation  and  that  the  Ionized  air 
expands  more  slowly  than  the  blast  wave.  From  the  absorption  length 
measured  from  the  streak  camera,  the  blast  wave  energy  was  found  to  be 
of  an  equal  magnitude  to  the  laser  energy  In  the  detonation  regime.  The 
Impulse  was  again  found  to  be  Independent  of  material  and  due  primarily 
to  the  high  pressure  gas  volume  formed  at  the  surface. 

Fowler  and  Smith  14]  reported  the  use  of  a  6kW  CW  C02  laser  In  a 
study  of  plasma  Initiation  on  metal  targets.  High  speed  photography 
showed  that  nearly  complete  decoupling  of  the  laser  radiation  occurred 
as  the  plasma  expanded  Into  the  6kW,  f/7  beam.  The  results  also 
Indicated  that  the  Initiation  In  air  occurred  analogously  to  an  arc 
breakdown  and  that  the  plasmas  extinguished  at  f/10  and  higher.  At 
lower  f/#s  where  the  plasma  could  be  maintained,  a  calorimeter  was 
placed  behind  a  small  hole  In  the  target.  Calorimetry  Indicated  that 
approximately  50  percent  of  the  laser  radiation  was  transmitted  by  the 
plasma.  Combined  with  the  photographic  studies,  a  volume  averaged 
absorption  coefficient  was  calculated  which  Inferred  an  average  plasma 
electron  temperature  of  14,000  K.  Minimum  plasma  maintenance  power  was 
less  than  2kW  (d  *  .07  mm)  which  agreed  with  predicted  values. 

The  paper  also  reports  on  the  threshold  Intensity  of  various 
materials.  For  spot  sizes  less  than  0.3mm  diameter,  the  threshold  was 
material  Independent.  At  spot  sizes  In  excess  of  1mm  diameter,  the 
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carbon  threshold  was  slightly  higher  than  aluminum,  a  result  that 
appears  to  be  counterintuitive  but  no  explanation  was  given.  Transverse 
gas  assists  had  little  effect  on  initiation  for  spot  sizes  smaller  than 
.5mm  diameter.  At  larger  spot  sizes,  the  transverse  gas  assist 
monotonlcally  Increased  the  threshold  with  spot  size,  between  25  and  200 
percent.  The  effort  also  utilized  a  two  wavelength  Interferometer  to 
measure  the  electron  number  density  in  the  air  plasma.  The  calculated 
number  densities  were  sufficient  to  suspect  refraction  or  blooming  was 
an  effect  on  the  delivered  intensity  to  the  target.  A  detector  behind  a 
hole  In  the  target  measured  the  energy  during  plasma  formation.  The 
results  showed  that  as  much  as  90  percent  attenuation  occurred  which 
disagrees  with  the  50  percent  value  from  calorimetry.  The  conclusion 
was  that  blooming  indeed  further  reduced  the  Intensity  and  altered  the 
spatial  distribution  of  energy  deliverable  to  the  target.  This 
conclusion  thus  contradicts  some  of  the  previous  effort's  conclusions  as 
the  various  diagnostic  techniques  did  not  necessarily  account  for 
blooming. 

Engel  [51  reported  the  status  of  COg  laser  cutting  of  thin 
materials  with  a  lkW  TEMgg  laser  with  Intensities  on  the  order  of 
10®W/cm2.  An  empirical  result  was  that  optimal  cutting  (minimum  kerf 
width  and  lateral  damage)  occurred  when  the  laser  focus  was  located  1/3 
the  material  thickness  below  the  surface.  In  metals,  Engel  concluded 
that  surface  finish  and  oxides  had  little  effect  on  cutting  rates.  The 
report  provides  one  of  the  first  empirical  relationships  for  cutting 
rate: 
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Rate  (cm/sec)  3  KP/(EaAt)  (1.1) 

where  K  Is  an  empirical  constant  (cm/sec),  P  Is  the  laser  power  (W),  Ea 
Is  the  ablation  energy  (W/cm3),  A  Is  the  spot  area  (cm2),  and  t  the 
material  thickness  (cm).  Engel  also  investigated  the  use  of  assist 
gases.  The  use  of  oxygen  increased  cutting  rates  by  40  percent  over 
inert  gases.  The  effect  was  attributed  to  exothermic  reactions  which 
aided  material  removal.  For  non-metals,  the  cutting  rates  demonstrated 
were  as  much  as  two  orders  of  magnitude  greater  than  conventional 
techniques  for  quartz.  The  laser  was  also  shown  to  improve  the  cut 
quality  in  difficult  to  machine  materials  such  as  Kevlar  and  the  lack  of 
tool  wear  made  these  processes  competitive. 

Robin  and  Nordln  (6)  Investigated  the  possibility  of  improving 
penetration  depths  by  tailoring  the  laser  energy  profile.  It  was  noted 
that  the  laser  energy  mainly  went  into  the  heat  of  vaporization  once  a 
melt  formed.  Mechanical  means  had  been  used  to  remove  the  melt  and 
maximize  the  solid/liquid  propagation  velocity.  From  the  momentum 
transfer  work  cited  previously,  the  authors  proposed  that  the  impulse 
delivered  by  a  pulsed  laser  could  be  used  to  drive  the  melt  out  of  the 
kerf  during  CW  cutting.  The  melt  was  assumed  to  escape  via  the  pulsed 
laser  Impulse  In  a  hollow  cylindrical  column  around  the  CW  focus 
geometry.  The  melt  mass  removal  was  assumed  total  and  the  rate  was 
integrated  to  determine  pulse  duration  and  thus  the  impulse  required. 
The  minimum  Impulse,  Im,  was  calculated  to  be  (6): 

Im  «  1/2  kpA(R/S)(l-s).  (1.2) 
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Where  k  is  the  thermal  conductivity,  A  is  the  spot  area  (=*R  ),  o  is  the 
melt  density,  S  Is  the  cylindrical  ejection  path  wall  thickness, 
and  8(=k/[k  +  VQh 1 )  is  a  measure  of  the  percentage  of  the  maximum 
solid/melt  interface  velocity.  The  approximation  is  valid  for  most 
metals  where  8  Is  25  percent  below  unity  for  a  given  melt  thickness, 
h.  The  result  for  aluminum  was  that  an  impulse  of  2  dyne-sec  or  a 
microsecond  pulse  at  108  to  lO^W/cm^  (20-100  dyne-sec)  would  be 
sufficient  for  a  substantial  Improvement  in  cutting  rates.  A  reference 
to  an  experimental  paper  showed  that  the  applied  impulse  was  sufficient 
to  open  a  hole  In  a  thin  target  which  remained  closed  in  the  case  of  no 
superimposed  pulse.  Other  experiments  also  Indicated  a  doubling  of 
cutting  rates  when  a  small  pressure  difference  was  maintained  across  a 
thin  target. 

Duley's  classical  text  1 7 J  on  carbon  dioxide  lasers  in  1976  marked 
the  first  encompassing  treatment  of  practical  laser  systems.  The  text 
covers  laser  physics,  heating,  drilling,  welding,  diagnostics,  and 
communication  aspects  of  COg  lasers.  Plasma  formation  is  treated 
qualitatively  In  the  various  applications.  It  was  noted  that  a  laser 
supported  combustion  wave  (LSC)  could  completely  shield  a  workpiece  from 
10.6  pm  radiation.  It  was  further  noted  that  plasma  effects  were 
adequately  circumvented  through  mechanical  means.  The  author  suggests 
the  use  of  a  pulsed  laser  system  with  a  pulse  separation  long  enough  to 
allow  the  plasma  to  extinguish.  Ho  report  of  experimental  verification 
was  made.  The  suggestion  will  be  shown  to  be  limited  as  several 
subsequent  Investigations  demonstrate  plasma  formation  to  be  desirable. 
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Metal  drilling  processes  were  treated  by  von  Allmen,  et  al.  [81.  A 
pulsed  Nd-YAG  laser  was  utilized  and  the  decrease  in  material 
reflectance  with  temperature  was  measured.  The  change  in  reflectance 
was  attributed  to  the  metal-dielectric  transition  during  melting  and  the 
formation  of  a  stationary  plasma  layer  at  the  surface.  The  sample  was 
placed  in  one  focus  of  an  ellipsoid  and  the  detector  placed  in  the 
other.  Three  phases  of  reflectivity  were  monitored.  First,  the 
reflected  signal  remained  constant  until  the  melt  began  to  form. 
Second,  a  fast  decrease  in  reflectivity  occurred  which  was  inversely 
proportional  to  laser  Intensity.  The  final  range  was  a  constant  low 
value  with  a  slight  decrease  as  the  Interaction  zone  Increased,  l.e., 
drilling  effects  became  significant.  The  reflectivity  recovers  to 
nominal  In  10  to  20  nanoseconds  after  the  laser  pulse,  depending  on 
target  damage. 

At  times  in  excess  of  720  nanoseconds  ,  the  authors  noted  the 
beginning  of  an  explosion-like  event  and  development  of  a  surface 
absorption  layer.  From  the  reflectivity  data  the  radial  expansion 
velocity  was  determined  to  be  104cm/s.  Enhanced  coupling  (reduced 
reflectivity)  was  always  noted  upon  the  initiation  of  a  plasma.  When 
the  mass  removed  was  compared  with  previous  drilling  models,  the  results 
Indicated  that  70  to  80  percent  coupling  occurred  which  agreed  with  the 
reflected  signals.  The  transition  temperature  for  reduced  reflectivity 
was  determined  to  be  a  few  hundred  degrees  above  the  vaporization 
temperature  for  Cu,  Al,  N1,  and  Fe.  The  authors  concluded  from  the 
plasma  volume  that  conduction  was  the  primary  transport  mechanism 
between  the  plasma  and  target  as  a  surface  temperature  conduction  model 
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adequately  explained  the  experimental  results.  From  the  conduction 
model  it  was  determined  chat  the  plasma  formation  corresponds  to  a 
relatively  low  evaporation  pressure  of  3  to  10  atmospheres.  This  range 
Is  significantly  below  the  ablations  pressures  proposed  by  other  authors 
which  are  on  the  order  of  ore  hundred  atmospheres. 

McKay  and  Schriempf  [3|  proposed  a  thermal  transient  analysis  of 
the  heating  profiles  of  thin  targets  as  a  method  of  measuring  an 
effective  coupling  coefficient.  The  authors  used  a  5  us  pulsed  C02 
laser  with  a  maximum  pulse  energy  of  140J.  The  analysis  was  considered 
for  pulse  times  and  target  sizes  such  that  lateral  heat  flow  was 
negligible  compared  to  the  normal  component.  With  this  assumption,  the 
target  thermal  respor.se  could  be  modeled  via  Laplace  transforms  which 
depended  on  the  material  properties  and  beam  parameters.  The  thermal 
response  was  measured  at  two  locations,  one  on  focus  center  and  the 
second  off  axis.  The  two  responses  could  then  be  applied  to  adjust  the 
thermal  response  model  parameters  to  effective  values.  At  low  fluences 
where  no  plasma  formed,  the  effective  fluence  deposition  and  hence 
absorbtlvlty  could  be  determined.  The  thermal  response  agrees  very  well 
until  late  In  the  profile  (greater  than  a  few  seconds)  when  the  finite 
size  of  the  target  made  the  modeling  assumptions  Invalid.  When  a  plasma 
was  formed,  both  the  fluences  and  beam  area  had  to  be  adjusted  to  fit 
the  thermal  responses.  The  conclusion  was  that  plasma  spreading  and 
beam  blooming  were  considerable  effects.  One  shouTd  note,  however,  that 
the  thermal  lag  of  the  thermocouples  was  apparently  not  considered,  an 
effect  that  would  be  very  critical  when  transient  slope  matching  Is  the 
basis  of  the  diagnostic  tool. 
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McKay,  et  al.  [Ill  continued  the  pulsed  interaction  studies  with 
the  implementation  of  ether  diagnostics.  Care  was  taken  to  renew  the 
surface  of  Al  targets  between  shots  by  abrasion.  The  result  was  that  4 
micrometer  dig  and  1  micrometer  Al  particles  were  left  after  abrasion. 
The  authors  admit  the  crltlcallity  of  the  "defect"  heating  in 
initiation,  but  were  more  interested  in  comparison  of  the  results  from  a 
separate  investigation.  Walters,  et  al.  [Ill  had  proposed  the  defect 
heating  theory  of  plasma  initiation  in  an  attempt  to  explain  the 
variation  in  reported  threshold  intensities  from  author  to  author.  The 
theory  is  generally  accepted  but  being  stochastic  is  difficult  to  verify 
or  model . 

The  results  of  McKay's  effort  [101  was  that  the  effective  spot 
radius  increased  50%  as  discussed  In  Ref.  9.  Below  the  plasma  threshold 
the  absorptivity  of  Al  was  9  percent  and  jumped  to  18  percent  upon 
plasma  initiation.  Spectroscopic  studies  of  the  gas  volume  near  the 
target  surface  revealed  that  the  unabraided  surfaces  spectra  were 
dominated  by  the  ambient  air  N  and  0.  For  the  abraided  surfaces  N  and 
Al  dominated  the  spectra  which  may  imply  the  removal  of  oxides.  The 
time  averaged  temperature  for  the  Al+  to  Al++  ratio  was  12,000  K 
compared  to  16,000  K  for  a  nitrogen  line  pair.  The  discrepancy  was 
attributed  to  the  target  lines  appearing  after  the  ambient  plasma 
formation  when  cooling  processes  within  the  plasma  have  begun. 
Interferometric  measurements  showed  the  electron  number  density  to  be 
10*®cnf^  and  the  blast  wave  expansion  velocity  on  the  order  of 
105cm/s. 


The  authors  concluded  from  a  simple  analysis  that  a  blackbody 
plasma  at  16,000  K,  xmax  =  220  nm,  and  a  target  absorptivity  of  40 
percent,  could  not  sufficiently  explain  the  enhanced  coupling  as 
measured.  The  33*  discrepancy  had  to  be  made  up  by  conduction  which  is 
possibly  the  dominant  effect.  The  conclusion  was  that  the  HAZ  doubled 
at  threshold  intensities  and  Increased  up  to  nine  times  the  area  at  peak 
powers. 

Maher  and  Hall  [12]  were  also  interested  in  pulse  tailoring.  The 
authors  noted  that  the  thermal  profile  radius  increased  with  increasing 
pulse  energy  but  did  not  elaborate  if  the  energy  was  in  the  pulse  spike 
or  in  the  tail  region.  From  previous  experiments,  the  authors  concluded 
that  optimal  laser  coupling  occurs  when  a  plasma  forms  near  the 
surface.  For  this  case,  the  primary  transport  mechanism  was  assumed  to 
be  UV  radiation  from  the  plasma.  Thus  a  pulse  with  an  initial 
spike  (1  us)  to  initiate  the  plasma  followed  by  a  lower  power 
(5xl05W/cm2)  plateau  to  maintain  the  plasma  should  result  in  optimal 
coupling.  The  results  Indicated  that  no  optimal  pulse  shape  could  be 
found  in  direct  contrast  to  the  results  of  Ref.  [10].  The  authors  also 
looked  at  the  production  of  A1  lines  during  the  plasma  formation.  The 
increasing  strength  of  A1  In  the  plasma/volume  with  increasing  laser 
power  was  Interpreted  as  an  Increased  amount  of  laser  energy  going  into 
vaporization  and  hence  a  higher  plasma  absorption. 

Carlhoff,  et  al.  [13]  Investigated  continuous  optical  discharges  in 
pure  gases  at  very  high  pressures  (1-210  bars).  The  experiments 
involved  a  2kW  CW  CO2  laser  with  intensities  In  the  range  of  106  to  107 
W/cm*.  Spectroscopy  of  A r  and  He  continuum  and  a  Mach-Zender 
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Interferometer  in  the  finite  fringe  mode  were  used  to  characterize  the 
plasma  electron  parameters.  The  qualitative  results  indicated  that  the 
peak  temperature  volumes  lie  upstream  from  the  laser  focus.  Maintenance 
power  increased  with  pressure  for  low  f /#  systems  and  decreased  for  high 
f/#  systems.  The  plasmas  were  stable  with  He  being  more  stable  than  Ar 
at  pressures  above  100  bar.  The  effort  demonstrated  that  high  pressure 
plasmas  could  be  maintained  relatively  easily  in  contrast  to  the  then 
current  theory.  No  further  investigation  was  presented  and  the 
direction  of  the  research  was  not  clear. 

Rosen,  et  al.  (14J  investigated  the  coupling  of  ultraviolet  laser 
radiation  (XeCl  at  .35  urn)  to  aluminum  targets  in  vacuum.  A  model  of 
the  plasma  was  developed,  including  the  effects  of  photoionization. 
Other  mechanisms  found  to  be  important  were  the  ionization  of  excited 
aluminum  and  the  elastic  collisional  deexcitation  leading  to  a 
temperature  difference  between  the  electrons  and  heavy  particles.  The 
model  results  indicate  that  the  aluminum  number  density  drops  in  10"6 
seconds  and  in  the  case  of  an  ambient  gas  presence,  inverse 
bremsstrahlung  absorption  would  be  governed  by  the  non-aluminum 
species.  The  model  agreed  with  the  measured  initiation  intensity  of 
pulse  durations  of  10~7  to  10_®  sec,  demonstrating  photoionization  being 
important  at  UV  laser  wavelengths.  For  2024  Al  the  as  received  coupling 
coefficient  was  determined  to  be  0.4  and  ground  and  polished  2024A1  was 
seen  to  be  0.2.  For  99.999  percent  pure  Al,  the  coupling  coefficient 
was  0.08  to  0.11  over  the  range  of  20  to  5  J/cm2.  The  plasma  threshold 
was  measured  to  be  in  excess  of  6xl07W /cm2  and  decreased  with  increasing 
pulse  duration,  as  predicted  by  the  model. 
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At  this  point  in  time,  late  1981,  much  of  the  American  literature 
was  classified  or  had  been  redirected  into  the  area  of  laser  fusion. 
Several  non-American  efforts  were  made  and  new  groups  began  to 
investigate  the  concept  of  laser  propulsion. 

One  of  the  foreign  efforts,  and  now  one  of  the  leading  Institutes 
of  industrial  laser-target  Interaction  was  presented  by  Miyamoto,  et  al. 
[15J  from  Osaka  University.  This  paper  attempted  to  quantify  the  role 
of  the  gas  assist  in  laser  welding  noting  the  poor  reproducibility. 
Several  gas  assist  parameters  were  investigated  beginning  with  nozzle 
angle  of  attack.  It  was  noted  that  the  penetration  depth  increased  with 
increasing  nozzle  angle  (Increasing  component  normal  to  target 
surface).  The  effect  saturated  at  approximately  30  degrees  and 
subsequent  experiments  were  performed  at  40  degrees.  Three  regimes  of 
pressure  dependence  were  found.  First,  the  plasma  bead  regime  at  low 
gas  pressure  was  found  to  result  In  a  shallow  bead  deemed 
unacceptable.  In  the  second  regime  (sub-atmospheric  pressure)  a  sound 
and  uniform  bead  resulted  with  an  increased  penetration  depth.  The 
third  range  resulted  in  the  bead  being  blown  off  the  weld  and  large 
porosities  were  formed  during  solidification.  A  smaller  nozzle  diameter 
and  placing  the  nozzle  closer  to  the  surface  resulted  In  a  larger 
operational  range  (second  regime).  The  cavity  in  region  two  was  located 
at  the  leading  edge  of  the  pool  and  periodic  changes  In  cavity  shape  led 
to  eddy  behavior  In  the  molten  flow  behind. 

The  effect  of  using  a  gas  with  a  higher  Ionization  potential  was  an 
increase  In  penetration  depth  for  the  case  of  zero  dynamic  pressure 
exerted  on  melt  pool.  Eventually,  pressure  overwhelmed  the  ionization 
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potential  effects  at  pressures  above  half  an  atmosphere.  The  effort 
also  demonstrated  the  use  of  phototransistors  for  the  monitoring  of  the 
weld  regime.  Conclusions  were  that  suppression  of  plasma  formation  was 
critical  and  at  low  speeds,  vacuum  welding  was  preferable. 

Another  text,  edited  by  Bass  [161,  treats  laser  material  processing 
on  a  more  general  level.  The  various  applications  of  welding,  cutting, 
etc.  are  treated  individually,  to  some  extent  comparing  the  use  of 
different  lasers.  Again,  however,  the  plasma  formation  is  treated  on  a 
semi-quantitative  basis.  In  the  chapter  on  laser  cutting  by  Steen  and 
Kamalu,  efficient  nozzle  design  is  noted  as  a  prerequisite.  The  authors 
discuss  the  effects  of  shock  wave  formation  from  the  nozzle  exhaust 
resulting  in  poor  cut  quality.  Little  or  no  reference  of  plasma 
formation  other  than  intermittent  plasma  formation  may  further  degrade 
cut  quality.  Only  now  have  studies  begun  to  investigate  the  gas  dynamic 
effects  of  shocks  on  target  transport.  Also  in  Bass,  the  chapter  by 
Mazumder  on  laser  welding  reports  qualitatively  on  scattering  and 
absorption  of  laser  radiation  as  the  beam  propagates  through  a  plasma 
volume.  Effective  assist  gas  systems  were  shown  to  be  capable  of 
circumventing  the  problem.  It  was  further  suggested  that  a  shield  gas 
with  a  higher  ionization  potential  would  serve  to  inhibit  plasma 
formation,  as  was  indicated  in  several  experimental  efforts. 

Another  effort  on  CW  CO2  material  processing  was  presented  by 
Oonati,  et  al.  [17J.  The  experiments  involved  the  use  of  a  15kW  C02 
laser,  with  an  f/7  telescope  resulting  in  a  focus  spot  diameter  of 
0.7mm.  High  speed  photography  and  C02  laser  transmission  experiments 
were  used  to  monitor  plasma  formation.  The  transmission  experiments 


Indicated  that  plasma  formation  was  Intermittent,  appearing  to 
extinguish  and  re-form  in  response  to  temporal  laser  and  gas 
fluctuations.  From  the  authors  data,  the  plasma  appears  to  decouple 
from  the  target  at  a  frequency  between  300  and  4000  Hz.  However,  the  CW 
interaction  times  may  have  been  relatively  short  (less  than  1  sec)  and 
the  transmission  data  appears  to  approach  a  DC  value  near  the  end  of  the 
interaction  time.  The  effect  of  the  plasma  coupling  data  was  not 
evident  in  the  absorbed  energy  measurements.  If  the  plasma  is  indeed 
decoupling,  it  is  doing  so  at  a  relatively  high,  random  frequency  with 
the  net  target  effects  being  the  result  of  a  smooth  integrated  effect. 

1.2  Modeling  Efforts 

Ouring  the  course  of  the  previously  cited  experimental  efforts, 
several  attempts  were  made  to  model  initiation  and  plasma  propagation. 
Most  of  these  efforts  expanded  upon  the  early  work  of  Raizer,  whose 
classical  treatment  of  blast  wave  theory  will  be  treated  in  Chap.  2. 

Wu  and  Plrri  [181  treated  the  breakdown  of  contaminated  air  at 
intensities  from  10®  to  107W/cmz  at  10.6  ym.  The  process  was  modeled  in 
three  temporal  stages.  First,  heating  and  vaporization  of  the 
contaminant  material  in  air.  The  contaminant  material  was  necessary  for 
absorption  to  occur  at  room  temperature  and  pressure.  Second,  the 
absorbed  energy  Is  transferred  to  the  surrounding  air  and  finally  an  air 
plasma  bubble  forms.  To  obtain  a  solution,  several  assumptions 
regarding  vaporization  of  the  contaminants  were  made  and  a  perturbation 
analysis  of  the  energy  and  propagation  equations  were  solved.  As  the 
plasma  expands  at  later  times,  a  snowplow  model  was  developed  in  which 
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the  plasma  front  pushes  away  the  ambient  air.  Any  mass  transfer  across 
the  front  was  assumed  to  Immediately  equillibrate  with  the  higher 
pressure  plasma  volume.  The  results  indicated  that  plasma  formation  Is 
nearly  complete  In  45  ys  which  agreed  with  an  experimental  effort.  The 
plasma  peak  temperature  rose  and  stabilized  at  17,800  K  after  about 
1  ys,  and  peak  temperatures  rose  with  increasing  intensity.  Finally, 
plasma  growth  rate  (volumetric)  Increased  with  increasing  Intensity  for 
both  the  early  stage  and  snowplow  models. 

More  recently,  Muller  and  Uhlenbusch  (19)  developed  a  model  for 
continuous  optical  discharges  In  conjunction  with  Ref.  (13J.  The  paper 
considers  two-dimensional  heat  flux,  one  dimensional  Imposed  flow  field 
with  a  constant  mass  flow  rate.  The  beam  absorption  of  a  focused 
Gaussian  profile  was  also  taken  Into  account.  The  authors  made  the 
simplification  of  Cp/K  is  constant,  which  is  not  valid  at  these 
temperatures,  and  also  utilized  a  volume  averaged  absorption 
coefficient.  Even  with  the  questionable  assumptions,  the  results  showed 
reasonable  agreement  with  the  measured  temperature  fields  of  Ref.  [13] 
at  5  bar  argon.  The  results  Indicated  a  weak  dependence  on  focused  beam 
quality  and  a  stronger  dependence  on  flow  velocity. 

The  brief  treatment  of  the  early  models  Is  somewhat  warranted. 
Almost  all  the  models  made  various  simplifying  assumptions  in  order  to 
obtain  analytical  solutions.  In  the  late  1970' s.  Interest  was  renewed 
In  the  concept  of  laser  propulsion.  With  the  experimental  work  came 
several  numerical  models  which  Incorporated  real  gas  properties  and  more 
recently  full  two-dimensional  analysis.  The  literature  survey  Is 
completed  by  a  review  of  the  field  of  laser  propulsion. 
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1.3  Laser  Propulsion 

Many  of  the  works  on  pulsed  laser  interaction  and  momentum  transfer 
evolved  Into  the  study  of  laser  propulsion.  A  separate  group  of 

Investigators  began  to  experimentally  study  continuous  wave  laser 

propulsion  In  the  early  1980's.  The  concept  Is  that  a  remote  laser 

source,  either  satellite  or  ground  based,  be  directed  and  focused  into 

and  absorption  chamber  onboard  a  payload.  Direct  heating  of  the 
propellant  occurs  upon  plasma  formation  and  the  hot  gas  subsequently 
expanded,  converting  thermal  energy  to  kinetic  energy  or  thrust.  The 
relatively  simple  absorption  chamber  eliminates  much  of  the  mass 
Involved  In  classical  combustion  processes.  Since  no  combustion  Is 
Involved,  low  molecular  weight  propellants  can  be  used,  further  reducing 
propulsion  mass.  The  net  effect  is  that  moderate  thrust  levels  at  high 
specific  Impulses  are  feasible.  The  current  application  appears  to  be 
low-earth  to  geosynchronous  orbital  transfer.  A  comprehensive  review  of 
history  and  status  of  laser  propulsion  was  made  by  Glumb  and  Krier 
(201. 

Part  of  the  research  to  be  discussed  In  this  study  and  In  Ref.  [21] 
deals  with  the  use  of  a  lOkW  CW  COg  laser  (f/2.4-3.4)  for  initiation  and 
maintenance  of  atmospheric  pressure  argon  plasmas  at  low  flow  rates. 
The  results  Indicated  that  the  plasma  volume  grew  with  Increasing  power 
resulting  In  a  net  absorption  of  80  percent  of  the  Incident  laser  power 
at  7kW.  The  absorbed  power  was  measured  via  a  copper  cone  calorimeter 
which  Intercepted  the  transmitted  laser  beam.  Thermocouple  measurements 
of  the  downstream  temperature  profiles  were  used  to  calculate  the  net 
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power  retained  by  the  gas  as  thermal  energy.  The  results  Indicated  that 
thermal  efflclences  In  excess  of  25  percent  were  obtained.  The  trend  of 
the  absorption  and  efficiency  curves  were  corroborated  by  results 
calculated  from  the  spectroscopically  measured  temperature  fields  at  the 
corresponding  laser  powers. 

The  Impact  of  this  study  and  similar  research  at  the  University  of 
Tennessee  Space  Institute  which  will  be  discussed  subsequently,  were 
that  plasma  initiation,  an  early  concern,  was  in  fact  routine.  Both 
studies  utilize  metallic  targets  injected  into  the  focal  region  to  lower 
the  breakdown  Intensity  by  the  creation  of  the  first  free  electrons. 
The  electron  densities  created  in  this  manner  have  been  found  to  be 
sufficient  for  cascade  breakdown  In  the  ambient  gas  to  occur.  A  second 
early  concern  was  the  ability  to  stabilize  a  plasma  In  the  systems  which 
has  again  proven  to  be  routine.  The  two  efforts  are  In  the  process  of 
mapping  larger  ranges  of  flow  and  pressure  parameters  to  attain  thermal 
efficiencies  up  to  50  percent  where  laser  propulsion  appears  feasible. 
The  University  of  Illinois  Is  also  looking  at  more  classical  means  of 
Improving  efficiency  such  as  regenerative  preheating  of  the  propellant 
and  multiple  plasma  formation  to  minimize  radiative  view  factors  to 
chamber  walls. 

A  simultaneous  modeling  effort  was  begun  at  the  University  of 
Illinois  by  Glumb  and  Krler  (221.  The  model  Is  numerical  with  a  one¬ 
dimensional  flow  field  with  two-dimensional  heat  transfer.  The 
converging  beam  geometry  and  real  gas  properties  for  argon  Incorporated 
In  that  model  were  at  the  time  unique.  The  results  agreed  with  the 
experimental  results  of  Ref.  [21} ,  and  continue  to  do  so.  The  model  has 


been  used  to  look  at  more  realistic  powers  and  pressures  In  an  attempt 
to  develop  scaling  laws  for  full-scale  propulsion  systems.  More 
Importantly,  the  model  has  been  used  to  direct  the  current  experimental 
work  In  the  Investigation  of  regions  of  optimal  thermal  efficiencies 
within  the  physical  limitations  of  the  current  laboratory  facilities. 

Welle,  Keefer,  and  Peters  (231  recently  reported  the  Tennessee 
effort  at  lower  laser  powers  (sub-kW)  but  higher  flow  rates  (.4  to  4.5 
m/s).  Earlier  results  demonstrated  similar  trends  as  the  U  of  I  effort 
but  particular  Importance  was  placed  on  beam  geometry.  The  latest 
results  Indicate  a  lesser  dependence  on  beam  geometry  versus  flow  rate 
as  noted  In  Ref.  [19j.  All  of  the  Tennessee  data  has  been  calculated 
from  the  continuum  emission  of  argon  at  626. 5nm,  with  no  independent 
diagnostic  technique  having  been  reported  to  date.  The  results  appear 
reasonable  with  peak  temperatures  in  excess  of  15,000  K  and  thermal 
efficiencies  approaching  40  percent.  The  Increase  In  thermal  efficiency 
Is  as  expected  due  to  the  smaller  volume  and  generally  cooler  plasma 
formed  at  the  higher  flow  rates,  resulting  In  smaller  radiative 
losses.  The  study  also  notes  that  the  experimental  error,  assumed  to  be 
±10%,  resulted  In  a  net  difference  of  the  global  properties  of 
absorption  and  radiative  losses  of  a  few  percent,  as  will  be  shown  in 
this  effort. 

The  University  of  Tennessee  has  recently  presented  the  results  of  a 
numerical  modeling  effort  by  Jeng  and  Keefer  [24|.  The  model  currently 
Incorporates  laminar.  Incompressible  Navler  Stokes  equations  in  two- 
dimensions  and  geometric  ray  tracing  of  the  converging/diverging  laser 
beam,  as  well  as  real  gas  properties.  The  model  agrees  well  with 
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experimental  flowing  argon  temperature  fields  for  static  pressures  In 
excess  of  two  atmospheres.  The  model  was  compared  with  the  results  of  a 
quasi  two-dimensional  model  similar  to  that  developed  in  Ref.  [22J.  The 
results  Indicated  that  the  optically  thick  radiation  at  pressures  below 
two  atmospheres  was  not  accurately  modeled  by  the  diffusion 
approximation.  Apparently  by  trial  and  error,  it  was  found  that  by 
using  50  percent  of  the  radiative  conductivity,  agreement  between  the 
model  and  experiment  was  satisfactory.  Assessment  of  the  solution 
procedure  remains  to  be  performed  as  well  as  the  applicability  of  the 
radiative  conductivity  criteria  In  the  quasi  two-dimensional  models. 
The  model  Is  currently  being  reworked  to  account  for  beam  diffraction 
and  refraction  In  an  attempt  to  clear  up  some  minor  discrepancies.  A 
useful  outcome  of  this  model  may  be  the  ability  to  Impress  a  more 
realistic  flow  field  onto  the  quasi  two-dimensional  model  grid.  The 
predicted  effects  are  that  the  radial  component  Is  small  compared  to  the 
axial  velocity.  However,  the  axial  velocity  approaches  zero  at  the 
centerline  near  the  peak  temperature  contour  and  local  variations  in  the 
radial  velocity  may  be  significant. 

It  has  been  demonstrated  that  a  considerable  effort  has  been 
devoted  to  laser-plasma- target  interaction.  The  results  of  various 
authors  have  varied  significantly  In  threshold  values,  optimizing  laser 
parameters  and  dominant  transport  mechanisms.  The  regime  of  CW  laser 
interaction  In  Industrial  applications  has  been  largely  Ignored  except 
for  some  quantitative  studies  of  optimal  gas  assist  configurations. 
Many  of  the  efforts  that  quantitatively  defined  experimental  plasma 
parameters  did  not  calculate  transport  phenomena.  Due  to  these 


discrepancies  and  the  availability  of  accurate  material  transport  codes 
{25.26J  this  Investigation  proposes  to  experimentally  monitor  plasma 
parameters  and  perform  the  subsequent  plasma-target  transport 
calculations.  The  goal  of  this  effort  Is  complete  a  priori 
determination  of  nominal  laser  parameters  required  to  effect  a  given 
material  transformation. 


2.  PLASMA  PHYSICS 


When  a  laser  of  sufficient  energy  is  focused  into  a  gaseous  medium 
(lO^W/cm^)  or  onto  a  solid  target  (105W/cm2),  a  plasma  will  form.  The 
plasma  is  a  collection  of  particles  which  can  range  from  a  few  percent 
to  complete  ionization.  In  this  Investigation,  the  concern  is  for  low 
temperature  plasmas  where  the  particle  temperatures  range  from  5,000  to 
20,000  K.  The  plasmas  formed  are  at  or  near  atmospheric  pressures  which 
imply  equilibrium  number  densities  on  the  order  of  1016  to  10*7  cm~3  for 
the  heavy  particles  and  electrons  (27]. 

In  this  chapter,  the  physics  of  low  temperature  plasmas  is 
discussed.  Although  the  study  Is  of  steady  state  CW  plasma  formation,  a 
discussion  of  the  Initiation  and  breakdown  mechanisms  Is  warranted  and 
Is  presented  first.  Next,  the  mechanisms  of  laser  energy  transport 
within  the  plasma  volume  are  discussed.  The  transport  equation  is 
developed  with  emphasis  on  the  radiative  properties  of  the  plasma.  The 
propagation  of  the  plasma  Into  a  converging  laser  beam  Is  discussed  with 
respect  to  the  Instabilities  associated  with  the  steady  state  plasma 
volume.  Since  the  effort  Is  concerned  with  both  pure  gas  and  gas/metal 
vapor  plasmas,  these  two  cases  are  treated  separately.  The  conduction, 
convection,  and  radiative  transport  to  the  ambient  Is  discussed  next. 
Finally,  the  concept  of  local  thermodynamic  equilibrium  (LTE)  is 
treated.  LTE  Is  the  fundamental  assumption  made  throughout  the  study, 
without  which  plasma  characterization  becomes  unwieldy. 


2.1  Breakdown  Mechanisms  and  Threshold 

Two  mechanisms  contribute  to  the  Initiation  of  a  plasma.  First,  on 
the  basis  of  quantum  mechanics,  there  is  a  finite  probability  of 
multiphoton  absorption  by  a  bound  electron.  The  absorbed  photon  energy 
would  be  sufficient  to  ionize  the  atom,  freeing  the  initial  electron 
"seed".  The  second  mechanism  Is  cascade  or  avalanche  breakdown.  In 
this  case,  it  is  assumed  that  the  initial  free  electron  is  available  for 
photon  absorption  in  the  presence  of  an  atom  or  ion.  The  process  is 
called  Inverse  bremsstrahlung  absorption  and  is  inversely  analogous  to 
"brake  radiation"  emission  of  a  photon  by  the  heavy  particle  In  the 
presence  of  an  electron. 

Multiphoton  absorption  probability  can  be  shown  to  be  low, 
especially  at  10.6  micrometer  radiation.  For  example,  the  ionization 
potential  of  argon  is  15.8eV  and  the  energy  of  a  10.6  micrometer  photon 
is  0.124eV.  Thus,  128  photons  need  to  be  absorbed  simultaneously  at 
this  wavelength,  neglecting  loss  mechanisms,  for  ionization  to  occur  in 
the  cold  gas.  Experiments  have  shown  that  multiphoton  absorption  is 
indeed  negligible,  especially  at  pressures  near  atmospheric. 

In  avalanche  breakdown,  the  first  free  electron,  or,  more  exactly, 
a  large  number  of  free  electrons  are  necessary  to  Initiate  the  plasma 
and  overcome  loss  mechanisms.  In  the  ideal  case,  the  electron 
oscillates  with  the  electromagnetic  field  In  the  laser  focal  volume. 
The  electron  transfers  kinetic  energy  to  the  atoms  leading  to  an 
Increased  amount  of  energy  In  the  random  motion  of  the  heavy 
particles.  The  random  motion  enhances  the  effective  cross-section  of 
the  atom  which  facilitates  Inverse  bremsstrahlung  (IB)  absorption.  The 


net  result  is  that  the  average  electron  energy  remains  the  same. 
However,  if  the  intensity  of  the  laser  is  sufficient  to  overcome  the 
loss  mechanisms,  the  electron  will  eventually  acquire  enough  energy  to 
ionize  an  atom  during  collision.  The  process  results  in  two  slow 
electrons  which  continue  the  process  and  hence  the  avalanche 
phenomena. 

The  loss  mechanisms  consist  of  electron  diffusion,  recombination, 
and  both  elastic  and  Inelastic  collisions.  Electron  diffusion  is  where 
the  electron  moves  out  of  the  focal  volume  prior  to  sufficient  energy 
absorption  for  Ionization.  For  the  relatively  large  focal  volume 
dimensions  with  respect  to  the  electron  mean  free  path  at  atmospheric 
pressures,  macroscopic  diffusion  losses  can  be  shown  to  be  minimal.  At 
the  high  optical  frequencies  in  this  Investigation,  the  electron 
diffusion  Is  further  restricted  by  the  oscillating  field  reversal. 
Discrepancies  in  threshold  measurements  indicate  that  diffusion  may  be  a 
factor.  Zeldovich  and  Raizer  ( 28]  hypothesize  that  local  "hot  spots"  in 
the  focal  volume  due  to  diffraction  and  aberrations  are  the  sources  of 
avalanche.  These  volumes  are  small  compared  to  the  focal  volume  and  may 
depend  heavily  on  electron  diffusion.  However,  the  stochastic  nature  of 
the  isophotes  due  to  laser  mode  changes  make  the  hypothesis  difficult  to 
verify  or  model. 

Recombination  Is  a  highly  probable  process,  especially  in  gases  at 
room  temperature.  The  elastic  collisions  can  be  modeled  as  hard  sphere 
interactions  in  which  the  electron  gives  up  the  mass  ratio,  me/ma, 
energy  per  collision,  where  ma  and  mg  are  the  atom  and  electron  masses, 
respectively.  The  inelastic  collisions  only  become  important  as  the 


electron  energy  approaches  the  ionization  potential.  The  electron 
transfers  energy  to  excite  the  atom.  The  higher  energy  criteria  is  due 
to  the  relatively  large  energies  associated  with  the  first  excited 
states  of  inert  gases.  In  argon,  for  example,  the  first  level  is 
roughly  75  percent  of  the  Ionization  potential. 

The  result  of  the  above  absorption  and  loss  mechanisms  is  that 
there  is  a  threshold  effect.  The  effect  has  been  noted  since  the  early 
studies  of  arc  plasmas.  From  the  net  rate  of  Increase  in  electron 
energy  of  a  cold  gas,  taking  Into  account  elastic  collisions,  the 
maximum  electron  energy  can  be  expressed  classically  as: 

E  (eV)  *  5.8  x  1020  AI(W/cm2)/(u,  +  v2)  (2.1) 

maX  in 


where  Emax  is  the  maximum  energy,  A  is  the  atomic  weight,  I  the  laser 
intensity,  w  the  laser  frequency  (rad/sec),  and  \>m  the  effective  elastic 
collision  frequency.  Therefore,  for  ionization  to  begin,  Emax  must 
exceed  the  ionization  potential.  For  example,  argon  at  atmospheric 


pressure,  C02  laser  radiation,  and  assuming  an  initial  electron  density 
of  lO^cm-3,  Eqn.  (2.1)  shows  that  the  required  laser  intensity  is 
109W/cm2.  In  most  "clean"  gases,  the  lO^cm"3  is  a  number  which  comes 
from  pre-ionization  via  arc  discharge,  otherwise  the  IB  process  does  not 
occur  at  any  laser  Intensity. 

It  Is  seen  in  Eqn.  (2.1)  that  the  threshold  Intensity  (Ith)  is 
proportional  to  the  elastic  collision  frequency  and  hence  pressure.  As 
Is  expected,  Ith  decreases  with  increasing  pressure  (number  density) 
which  Increases  the  collision  frequency.  Other  parameters  play  a  role 
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in  the  threshold  effect.  First,  the  threshold  increases  with  ionization 
potential  as  expected.  An  unexpected  effect  is  the  threshold  peak  in 
the  visible  spectrum  as  one  were  to  change  the  laser  frequency  from  UV 
to  IR.  Again  from  Eqn.  (2.1),  one  expects  the  threshold  lowering  in  the 
infrared,  as  u  is  decreasing.  However,  experiments  have  shown  an 
increase  in  the  threshold  in  the  IR.  The  IR  values  are  too  low  for 
multiphoton  absorption  to  be  a  factor,  but  no  definitive  theory 
currently  exists. 

Another  factor,  as  mentioned  before,  is  the  focal  volume  in  the 
case  where  electron  diffusion  may  be  a  factor.  At  atmospheric 
pressures,  where  this  study's  interests  lie,  the  effect  is  negligible. 
The  effects  of  pulse  duration  are  not  a  major  concern  In  this 
investigation,  since  CW  radiation  is  utilized.  There  are  indications  of 
threshold  lowering  with  increasing  pulse  duration.  The  lowering  is 
marginal  and  was  concluded  by  comparing  the  results  of  several  authors 
under  varying  conditions.  It  should  be  noted  that  other  efforts  have 
reported  no  dependence  on  pulse  duration,  only  on  pulse  power. 


2.2  Bremsstrahlung  Absorption 

Inverse  bremsstrahlung  (IB)  absorption  is  the  basic  mechanism  of 
power  absorption  within  the  plasma  volume.  The  absorption  coefficient 
can  be  derived  via  classical  or  quantum  theory.  The  main  difference  in 
the  two  results  Is  that  the  quantum  derivation  maintains  separate  terms 
accounting  for  IB  absorption  and  losses  due  to  stimulated  emission  at 
the  laser  frequency.  The  classical  derivation  is  a  net  calculation  and 
the  individual  terms  are  lost  in  the  final  form. 
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The  IB  absorption  coefficient  was  first  derived  by  Kramers  [29]  in 
1923.  Kramers  used  classical  theory  in  the  derivation  of  the 
bremsstrahlung  emission  coefficient.  The  author  then  applied  the 
quantum  theory  of  detailed  balancing  to  arrive  at  the  IB  absorption 
coefficient  which  is  of  the  form: 
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where  N^  is  the  Ion  density,  Ng  the  electron  density,  and  Z  the  ion 
charge.  The  derivation  was  hydrogenic  and  was  based  on  calculations 
that  line  emission  energy  offset  continuum  reductions  of  the  actual 
hydrogen  spectrum  compared  to  Plancks  blackbody  radiation  relation. 

Equation  (2.2)  represents  the  electron-ion  coefficient.  The 
electron-atom  coefficient  plays  a  role  at  lower  electron  temperatures 
near  10,000  K.  As  the  ionization  exceeds  a  few  percent,  the  long  range 
Ion  Coulomb  interactions  overwhelm  the  relatively  short  range  atomic 
cross-sections.  Several  authors  have  attempted  to  account  for  atom-ion 
IB,  non- hydrogenic  behavior  of  certain  gases  and  stimulated  emission  by 
developing  effective  coefficients.  Notably,  Wheeler  [30]  made 
corrections  for  argon.  To  account  for  stimulated  emission,  Wheeler 
showed  that  the  net  effect  was  a  reduction  of  a  by  a  factor 
of  [1  -  Exp(-hv/kTg) ] ,  where  v  Is  the  laser  frequency.  To  account  for 
non-hydrogen 1c  behavior,  the  free-free  Gaunt  factor  (gff)  multiplies 
Eqn.  (2.2)  as  well  as  the  effective  charge,  Z.  A  corrected  "effective" 
Gaunt  factor  was  proposed  by  Stallcop  [31].  The  first  correction  of  g^ 
was  to  account  for  electron-atom  IB.  Secondly,  a  term,  hv/kTg,  was 
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added  to  to  account  for  photoionization  effects.  Although 
photoionization  was  shown  to  be  negligible  at  10.6  micrometers  during 
initiation,  the  effects  increase  to  a  non-negligible  value  when  excited 
species  are  present.  The  result  is  an  effective  inverse  bremsstrahlung 
absorption  coefficient: 


a(m_1) 


2ne&Z 

3  4 

/6wkm„  he 
e 


{ l-exp(-h\j/kTe)  ]  (gff+£|-)  (2.3) 


This  absorption  coefficient  has  been  successfully  used  in  a  numerical 
model  [221  for  atmospheric  laser  supported  plasmas  in  pure  flowing 


argon. 

Loss  mechanisms  such  as  electron  and  heavy  particle  diffusion  out 


of  the  laser  geometry  will  also  affect  the  net  absorption.  The 
breakdown  threshold  criteria  is  defined  by  Ne/NQ  *  1CT3  from  Hughes 
[32].  Once  the  plasma  has  been  initiated  (Ne/NQ  approaching  unity)  and 
particularly  at  atmospheric  pressures,  diffusion  losses  should  be 
minimal  due  to  collisional  effects.  Recombination  is  a  high  pressure 


loss  mechanism  which  is  approximately  accounted  for  in 
the  1  -  e~hv^l<T  term  of  Eqn.  (2.3). 

The  question  remains  as  to  whether  the  semi-classical  derivation  is 
valid  for  the  plasmas  In  this  investigation.  Hughes  [32)  notes  that  the 
classical  derivation  Is  valid  for  the  case  when  the  laser  frequency 
exceeds  the  plasma  frequency.  This  Implies  that  binary  collisions 
accurately  account  for  photon  absorption.  This  indirectly  implies  the 
existence  of  local  thermodynamic  equilibrium  (LTE)  as  a  Maxwellian 
distribution  is  prerequisite.  Thus,  if  one  considers  the  oscillating 
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motion  of  an  electron  in  the  E-field  of  the  laser,  the  electron  achieves 
some  maximum  velocity  and  hence  energy.  The  "quivering"  speed,  from 
Hughes  [32],  can  be  shown  to  be: 

Ve  .  |^|  =  25xT^  (2.4) 

1  m„w 1  e 

e 

where  x  is  in  micrometers  and  Te  in  degrees  K.  The  thermal  velocity  of 
the  electron  is  on  the  order  of: 

VT  -  /kym^  (2.5) 

From  this,  the  field  strength  parameter  can  be  defined: 


n « Ve/vr 

If  n  >  1,  the  laser  radiation  field  has  an  effect  on  the  electron  energy 
distribution,  reducing  the  validity  of  LTE.  For  cases  In 
which  n  «  I,  the  laser  radiation  has  no  effect  and  the  classical  IB 
absorption  coefficient,  Eqn  (2.3),  Is  valid.  At  a  Tg  of  20,000  K,  Vy  * 
5.5xl05  and  Ve  *  3.75xl04,  and  n  <  1.  At  a  Te  of  10,000  K,  the  ratio  is 
smaller  still  and  the  classical  approach  appears  valid. 

One  problem  In  the  above  simplified  treatment  is  that  there  Is  a 
finite  probability  of  multiple  photon  absorption  prior  to  stimulated 
emission.  The  result  Is  that  the  criteria  for  the  field  strength 
parameter  Is  relaxed  and  an  average  value  of  Vy  should  be  used.  Due  to 
the  stochastic  nature  of  the  above  process,  the  average  thermal  velocity 
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is  a  nebulous  entity  and  the  rule  of  thumb  is  that  n  is  less  than  or 
equal  to  one. 

An  interesting  artifact  of  the  difference  between  classical  and 
quantum  derivations  Is  that  there  is  an  optimal  intensity  value,  Iopt. 
At  this  value  the  plasma  will  absorb  a  maximum  amount  of  energy  from  the 
Incident  laser  radiation.  The  classical  absorption  coefficient  (a), 
valid  In  the  weak  field  case,  shows  that  power  absorption  increases  with 
intensity.  The  quantum  derivation  accounts  for  non-linear  radiation 
interaction,  the  result  being  that  a  Is  Inversely  proportional  to 
I3/2.  Thus,  the  absorbed  power  decreases  with  increasing  intensity.  In 
the  modified  limit  of  the. field  strength  parameter,  n  3  1,  it  can  be 
shown  that  132): 

1  t  s  3  x  iq1V2  (w/cm2)  {2,7) 

where  Tfi  Is  in  eV.  This  generalized  term  may  be  of  some  value  in  the 
optimization  of  laser  propulsion  systems  thermal  efficiency.  It  will 
also  serve  as  a  value  to  avoid  In  high  energy  laser  material  processing 
applications  as  IQpt  at  20,000  K  and  10.6  micrometers  Is  4.85  x 
1010W/cm2. 

2.3  Transport  Equation 

The  radiative  transport  equation  accounts  for  the  change  in 
Intensity  as  a  ray  passes  through  an  emitting  and  absorbing  medium.  The 
equation  takes  the  form: 
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n2d(%  »  (-1  a  +  j  )dz  (2.8) 

'  2#  v  u  id  u» 
n 

where  n  is  the  refractive  Index,  I  ,  a  ,  and  j  the  frequency  dependent 

ID  ID  ID 

intensity,  absorption  coefficient,  and  emission  coefficient, 
respectively,  and  z  the  direction  of  propagation.  Defining  the  optical 

depth,  D  ,  and  source  function,  S  ,  as: 

or  w 


D  =  f  a  dZ 

“  AZ  “ 

(2.9) 

S  =  j  /a 

ID  ID  ID 

n 

(2.10) 

When  substituted  Into  Eqn.  (2.8),  the  result  Is  the  classical  form  of 
the  equation  of  radiative  transfer: 


I  -I 

q  * _ a  +  s 

dO  1  >  ■ 


id  n 


(2-11) 


If  the  Index  of  refraction  variation  Is  small  enough  over  dD^  (dz),  and 
the  source  function  Is  negligible  at  the  laser  frequency,  w, 
Eqn.  (2.11)  reduces  to: 


^OUT 


‘wIM  exp(-awZ). 


(2.12) 


Equation  (2.12)  Is  the  form  utilized  In  this  investigation  as  detailed 
In  Chapter  3  and  Appendix  B  for  the  calculation  of  absorbed  laser 
energy.  The  complete  radiative  transport  equation  Is  employed  In  the 
system  energy  balance  which  will  be  detailed  In  Chapter  5,  discussing 
the  numerical  modeling  effort. 


32 


The  radiation  terms  that  comprise  the  source  function,  S  ,  are 

01 

bound-bound,  free-bound,  and  free-free  transitions.  A  detailed 

treatment  of  the  above  terms  can  be  found  In  Griem  [33].  The  following 
Is  a  brief  description  of  the  physical  nature  of  the  various 
processes. 

Bound-bound  (b-b)  emission  occurs  when  an  electron  decays  between 
two  well  defined  energy  levels  of  the  system.  The  result  is  a  discrete 
line  spectrum.  A  given  transition  profile  depends  upon  the  upper  level 
spontaneous  lifetime  and  the  velocity  distribution  of  the  atoms  which 
can  be  expressed  via  quantum  physics  as: 

hC  Amn9mN 

h  *  “tf-VV  <2-13> 

Where  I  Is  the  line  Intensity  (W/cm3-sr),  at  wavelength  x,  Amn,  the 
transition  probability,  gm  the  upper  level  degeneracy,  Ze  the  electronic 
partition  function  (-1  for  argon),  Em  the  upper  level  energy,  and  N  the 
total  number  of  emitting  particles.  This  relation  applies  to  radiation 
at  a  pure  wavelength,  x.  In  practice,  the  line  profile  depends  on 

collisions  and  fields  within  the  plasma  volume  as  well  as  self¬ 

absorption.  One  broadening  mechanism  is  due  to  the  thermal  motion  of 
the  heavy  particles  or  Ooppler  broadening.  The  Doppler  effect  results 
In  a  Gaussian  line  profile  centered  at  x  and  the  half-width  is 

proportional  to  the  heavy  particle  temperature,  Tg.  A  second  effect  is 
Stark  or  pressure  broadening.  Stark  broadening  Is  due  to  Coulomb 

Interactions  and  Impresses  a  quadratic  profile  which  Is  proportional  to 
the  electron  number  density,  Ne.  The  two  broadening  mechanisms  are 
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difficult  to  resolve  and  one  treatment  can  be  found  in  Grlem  [331.  A 
general  rule  (321  is  that  Stark  broadening  dominates  the  profile  when 
the  Ionization  percentage  exceeds  0.1  percent.  This  is  due  to  the 
relatively  long-range  effects  of  the  Coulomb  forces.  In  this 
Investigation,  the  plasma  line  radiation  is  mainly  in  the  near  UV  to 
near  IR. 

The  second  radiation  source  Is  the  free-bound  or  recombination 
effect.  This  source  occurs  when  an  electron  is  captured  by  an  ionized 
atom  which  subsequently  decays  with  the  excess  energy  emitted  as  a 
photon.  The  spectrum  Is  a  continuum  with  a  low  frequency  cutoff  called 
the  recombination  limit.  The  limit  corresponds  to  the  minimum  energy 
required  to  Ionize  the  atom  and  was  seen  In  the  calculation  of  the 
bremsstrahlung  absorption  coefficient.  The  process  Is  difficult  to 
model  as  each  bound  state  has  a  unique  continuum  which  overlaps  and  the 
terms  must  be  summed  over  all  bound  states.  Also,  as  the  density 
Increases,  three  body  processes  compete  with  recombination.  Two 
electrons  can  collide  with  an  Ionized  particle  and  the  excess  energy 
carried  off  by  the  non-captured  electron  which  further  complicates 
theory. 

The  third  process  Is  free-free  or  bremsstrahlung  radiation.  The 
form  is  a  continuum  due  to  an  electron-heavy  particle  collision 
resulting  In  a  slow  electron  and  the  emission  of  a  photon.  In  the  hot 
regions  of  the  plasmas  electron- ion  collisions  dominate  and  in  the 
cooler  regions  (<10,000  K)  electron-atom  collisions  dominate  due  to  the 
relative  number  of  particles.  However,  as  shown  by  the  rough 
temperature  boundary  above,  pure  numbers  are  not  the  only  effect. 
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Electron-Ion  bremsstrahlung  Is  aided  by  the  long-range  Coulomb  fields 
which  effectively  Increase  the  collision  cross-section,  similarly  to  the 
case  of  Inverse  bremsstrahlung  absorption.  The  bremsstrahlung  continuum 
dominates  the  total  emission,  particularly  in  the  Infrared 
when  hv  <  0.7  kTfi.  The  relationship  arises  from  the  criteria  that  the 
emitted  photon  energy  must  be  less  than  the  local  electron  energy. 

Summing  the  continuum  terms  over  the  Internal  energy  levels  and 
applying  Kramers-Unsold  theory  1 29 1  yields: 

N  N 

Icx  *  5.44  x  lO"46  eg  Z2ff  -U  ^  (W/cm3-  sr)  (2.14) 

x 

Where  g  Is  the  Gaunt  factor  averaged  over  the  continuum  Interval  ax  and 
Zeff  Is  the  effective  charge.  Equation  (2.14)  was  calculated  using  the 
Planck  radiation  function  based  on  the  assumption  that  line  radiation 
energy  offsets  the  continuum  reduction  between  lines.  A  more  detailed 
calculation  of  the  continuum  emission  was  given  by  Dresvin  (34): 

N  N 

Icx  *  5.44  x  10'46  [exp(hcax/x2kTe)  +  exp(-hc/xkTe) (g  -  l)k(x) 

e  (2.15) 

The  first  exponential  term  accounts  for  recombination  to  widely  spaced 
lower  levels  and  the  second  exponential  describes  the  more  closely 
spaced  upper  level  recombination.  The  factor  c(x)  Is  a  term  derived  by 
Hibernian  (351  which  accounts  for  the  non-hydrogenlc  behavior  of  a  given 
system  and  has  been  calculated  for  argon. 

The  problem  Is  that  both  Eqns.  (2.14  and  2.15)  tend  to  oversmooth 
the  actual  continuum,  however,  the  effect  Is  minimal  In  total  integrated 
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calculations.  The  total  radiation  loss  calculations  are  difficult, 
especially  for  an  optically  thick  plasma  where  radiation  trapping  can 
occur.  In  general,  real  plasmas  are  treated  as  optically  thin  and  the 
optically  thick  transport  is  left  as  a  radiative  conductivity  term,  as 
has  been  done  In  this  effort.  The  spectral  regions  of  interest  are  as 
follows.  In  the  VUV,  the  radiation  is  trapped  as  the  resonant  lines  and 
continuum  are  of  comparable  magnitudes.  The  result  Is  that  this  region 
is  modeled  well  as  a  blackbody  as  long  as  the  pressure  Is  near 
atmospheric  and  the  temperature  distribution  Is  such  that  higher 
Ionization  effects  are  negligible.  The  IR  spectrum  Is  generally  opaque 
as  recombination  of  the  free  electrons  dominates  the  process.  A  general 
rule  In  the  UV  and  visible  spectrums  Is  that  absorption  Is  critical  only 
when  pi  >  10-50  atom-cm  where  p  Is  the  pressure  and  i  Is  a 
characteristic  length  of  the  plasma  volume.  In  the  spectroscopic 
studies  In  this  Investigation,  pi  is  typically  on  the  order  of  1-2  atom- 
cm,  however,  the  effects  of  self  absorption  are  treated  in  more  detail 
In  Chapter  3. 

The  radiative  conductivity  has  been  shown  by  Dresvln  [34]  to  be 
approximately: 

Kr  *  3  x  10'4L  T*  (2.16) 

where  L  Is  the  mean  free  path  of  a  photon  (Rosseland  length).  Equation 
(2.16)  can  be  used  over  small  Intervals  where  the  line  and  continuum 
magnitudes  do  not  differ  by  several  orders.  The  contribution  of  Kr  is 
large  In  the  2,000  to  16,000  K  range  and  agrees  well  with  experimental 
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results.  However,  experimental  results  have  been  used  in  this 
investigation  as  detailed  in  Chapter  5. 

2.4  Plasma  Propagation  and  Equilibrium 

Early  Investigations  reported  the  phenomena  of  plasma  propagation 
into  a  converging  laser  beam.  Propagation  velocities  were  estimated  to 
be  supersonic  on  the  order  of  10^m/sec.  Thus,  the  analogy  to  detonation 
or  blast  wave  phenomena  was  drawn.  Citing  the  analogy  between  laser 
plasma  propagation  and  combustion  wave  theory  is  generally  attributed  to 
Raizer  [36].  The  two  regimes  are  laser  supported  detonation  (LSD)  and 
laser  supported  combustion  (LSC)  waves.  In  LSD  propagation,  the  plasma 
propagates  at  supersonic  velocities  and  the  expansion  Is  radiatively 
driven.  In  LSC  propagation,  the  plasma  expands  at  subsonic 
velocities.  The  LSC  propagation  Is  slow  enough  for  conduction  to 
dominate  the  transport  mechanisms.  As  was  shown  in  the  literature 
survey,  much  of  the  early  work,  including  Raizer' s,  was  interested  in 
defining  the  plasma/laser  parameters  Involved  in  the  transition  regime 
between  LSC  and  LSD  behavior. 

In  LSD  propagation,  a  strong  shock  develops  near  the  plasma  focus 
and  rapid  heating  occurs  within  a  narrow  shock  front.  Nearly  all  the 
laser  energy  is  absorbed  In  the  narrow  expanding  front.  As  the  front 
expands  out  of  the  laser  beam,  the  plasma  abates  due  to  the  lack  of  an 
energy  source.  The  gas  behind  the  shock  also  relaxes  to  ambient 
conditions  for  similar  reasons.  Due  to  the  high  propagation  velocities, 
the  shock  front  Is  radiatively  driven  as  conduction  effects  are  only 
Important  outside  the  laser  beam. 


q 

The  early  experiments  were  performed  with  giant  laser  pulses  l(r  to 
1012  W/cm^.  It  was  soon  realized  that  plasmas  could  be  formed  at  lower 
intensities  via  an  external  source  of  electrons  to  initiate.  In  this, 
the  LSC  regime,  the  shock  expansion  does  not  accurately  explain  the 
propagation  and  plasma  state  behind  the  shock.  The  propagation  velocity 
is  subsonic,  heat  transport  is  primarily  due  to  conduction  and  the 
expansion  occurs  at  close  to  ambient  pressure.  The  subsonic  propagation 
velocity  is  proportional  to  the  local  laser  intensity  and  the  expansion 
front  slows  down  as  it  propagates  Into  the  converging  laser  beam.  The 
front  stops  propagating  at  the  point  where  the  propagation  speed  matches 
the  induced  flow  speed.  It  should  be  noted  that  the  local  mass  velocity 
can  be  co-dlrectlonal  with  the  propagation  velocity  If  the  expansion 
front  Is  strong  enough  to  "snowplow"  the  ambient  gas  to  temperatures 
well  above  ambient.  This  regime  was  thought  to  be  analogous  to  slow 
combustion  wave  theory,  however,  the  measured  expansion  velocities  were 
on  the  order  of  ten  times  larger  than  theory.  This  discrepancy  was 
later  explained  by  noting  that  the  pressure  across  the  expansion  front 
is  nearly  constant.  In  this  case,  the  density  ratio  across  the  front 
is: 


p  /p,  =  T  /T  and  If  Te»  10.000K  and  T  *  300K, 
a  s  s  a  s  a 

the  factor  of  30  accounts  for  the  theory  underestimate,  where  subscripts 
a  and  s  correspond  to  ambient  and  shock  conditions,  respectively. 

In  the  slow  combustion  wave,  upstream  and  downstream  regions  have 
little  effect  on  one  another.  The  downstream  plasma  has  relatively  low 
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thermal  gradients  compared  to  the  upstream  region.  Since  weak 
absorption  occurs  at  any  given  axial  location,  the  transmitted  laser 
energy  Is  sufficient  to  maintain  the  downstream  plasma. 

The  steady-state  energy  equation  is  two-dimensional  and  difficult 
to  solve.  An  analytical  one-dimensional  axial  solution  can  only  be 
determined  with  the  application  of  several  simplifying  assumptions. 
These  assumptions  include  linear  or  constant  thermophysical  properties 
and  a  parallel  (non-focused)  beam.  The  analytical  models  gave  way  to 
numerical  models  incorporating  realistic  beam  geometries  and  non-linear 
properties.  Most  notably  Glumb  and  Krler  [22 ]  developed  a  quasl-two- 
dimenslonal  model  which  accounts  for  real  gas  argon  properties  and  a 
converging  beam.  The  model  accounts  for  two-dimensional  conduction  and 
one-dimensional  (axial)  convection.  More  recently,  Jeng  and  Keefer  [24] 
demonstrated  a  full  two-dimensional  solution  of  the  incompressible 
Navier-Stokes  equation.  The  results  of  the  effort  are  impressive, 
however,  the  need  for  a  full  2-0  solution  In  the  case  of  buoyancy  driven 
plasmas  above  targets  is  questionable. 

The  present  effort  will  follow  the  work  of  Glumb  and  Krier  (22). 
With  the  assumption  of  zero  radial  velocity,  the  continuity  equation 
reduces  to  pu  =  constant.  Assuming  a  constant  pressure  flowfield,  the 
energy  equation  remains  to  be  solved: 

*  i  fc  l(rt£>l  -  .1  Cp  S  ♦  .1  -  PR  -  0  (2.17) 

Where  K  is  the  thermal  conductivity,  p  the  local  density,  u  the  axial 
velocity  (pu  *  constant),  Cp  the  specific  heat,  ol  the  absorbed  power. 
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and  Pr  the  radiative  loss  term.  It  should  be  noted  that  the  assumption 
of  a  constant  pu  value  may  be  unnecessary.  The  recent  effort  of  Jeng 
and  Keefer  (24)  has  utilized  a  full  two-dimensional  solution  of  the 
Navier-Stokes  equations  for  this  system.  It  may  be  possible  to  utilize 
the  results  of  the  radial  velocity  field  solutions  as  a  perturbation  in 
the  quasi-two-dimensional  model.  The  effect  is  treated  in  detail  in 
Chapter  5  as  well  as  the  energy  equation.  It  should  be  noted  that  in 
general  the  various  under/over-estimates  in  properties  mentioned 
previously  tend  to  cancel  one  another.  For  example,  the  underestimate 
in  continuum  loss  is  somewhat  compensated  by  the  VUV  preheating  of  the 
gas  which  is  not  accounted  for  in  Eqn.  (2.17). 

The  equilibrium  solution  of  Eqn.  (2.17)  depends  on  the  thermal 
conductivity  and  specific  heat.  For  an  equilibrium  composition,  three 
relations  for  each  specie  and  each  reaction  are  necessary.  First,  the 
various  reactions  need  to  be  accounted  for.  In  the  case  of  argon  at  Tg 
<  25,000  K,  only  two  reactions  are  applicable: 


ar  :  ar  +  e 


ar  :*r+e 


(2.18) 

(2.19) 


Equation  (2.18)  can  be  expressed  via  Saha  equilibrium  as: 


N  N  2Z,  2wmkT  3/2 

“RT  =  17  (_TTi) 


(E  +  &E) 


TT 


-1 


(2.20) 


Where  N+,  Ne,  and  NQ  are  the  Ion,  electron,  and  neutral  densities, 
respectively,  Z  is  the  electronic  partition  function,  and  E^  is  the 
ionization  potential.  When  combined  with  macroscopic  charge  neutrality 
and  the  law  of  partial  pressures,  the  composition  of  any  plasma  can  be 
determined.  Obviously,  the  number  of  equations  and  solution  becomes 
unwieldy  except  for  simple  monatomic  and  diatomic  systems. 

To  determine  the  specific  heat,  the  enthalpy  must  first  be 
calculated.  The  enthalpy  is  composed  of  several  effects  including 
kinetic,  excitation  and  ionization  for  argon  (monatomic). 


H  *  Hkin  +  "exc  +  H1on 
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Hkin  *  3/2  kT  l  (N^  +  ng)  from  statistical  thermodynamics  where  the 
summation  is  over  all  energy  levels.  Hgxc  *  £  N^,  where  E^  is  the 

excitation  energy  difference  between  states  i  and  j.  H.  „  =  E.  N.  , 
where  E^  is  the  Ionization  potential  summed  over  the  pertinent 
ionization  states.  Finally,  the  specific  heat  is  defined  from  classical 
thermodynamics  as: 


(2.22) 


It  should  be  noted  that  for  systems  Involving  non-monotonic 
species,  there  Is  a  dissociation  component  of  the  enthalpy.  This  factor 
typically  dominates  Cp  in  the  3,000-8,000  K  range.  Thus,  diatomic  gases 
are  sometimes  preferable  in  laser  welding  and  cutting  as  the  effect 
enhances  the  heating  capacity  of  the  laser/plasma  system. 


The  thermal  conductivity  requires  the  knowledge  of  all  pertinent 
elastic  collisional  cross-sections  since  at  low  temperatures  (.5  to  2 
eV)  elastics  dominate.  The  cross-sections  of  interest  in  a  monatomic 
system  are  the  electron-ion  and  electron-atom  terms.  The  e-e,  i-i,  a-a, 
and  i-a  cross-sections  are  typically  orders  of  magnitude  smaller  than 
the  e-i  and  e-a  terms.  From  gas  kinetic  theory: 

K  =  6.3  x  10'4  /TgTH  /<q>  (2.23) 

where  M  is  the  molecular  mass  (gm)  and  <q>  is  the  effective  cross- 
section  for  momentum  transfer  (cm  ).  Ambipolar  and  diffusive  transport 
are  critical  in  the  exact  calculation  of  the  thermal  conductivity: 


and 


K  “  *atom  +  K1on  +  Ke  +  ^reactive 


K  =  K  +  K 

reactive  ionization  excitation 


(2-24) 

(2.25) 


for  the  case  of  a  monatomic  gas.  In  general,  Ke  dominates  once  the 
ionization  is  above  a  few  percent.  An  empirically  based  relation  for  Ke 
was  shown  by  Dresvin  (37)  to  be: 


1.84  x  10'12  T  5/2 

K  =  _ ® _ 

e  1/4  m(A4  +  a^  +  e4) 


(2.26) 


where  Aj  and  a^  are  simple  functions  of  Tfi  and  Ng.  It  should  be  noted 
that  experimental  values  of  thermal  conductivity  are  greater  than  theory 
at  temperatures  above  10,000  K.  This  is  due  to  the  radiation  effects 


(radiative  conductivity)  as  discussed  in  the  previous  section.  The 
exoerimental  measurements  cannot  distinguish  the  effects  of  thermal  and 
radiative  conductivity.  Thus,  numerical  modeling  efforts  utilize  the 
experimental  values  and  could  implement  a  ratio  of  Eqns.  (2.26  and  2.16) 
to  attempt  to  specify  a  relative  magnitude,  if  necessary. 

2.5  Plasmas  in  Pure  Gases 

In  this  investigation,  pure  argon  plasmas  are  initially  investi¬ 
gated.  The  effort  is  in  support  of  a  study  of  the  use  of  laser 
supported  plasmas  for  the  direct  heating  of  a  propellant  for  propul¬ 
sion.  Argon  has  been  chosen  from  the  standpoints  of  well  documented 
thermophysical  properties  and  safety.  Argon  is  also  used  as  an  assist 
gas  in  several  laser  material  processing  applications.  Thus  diagnostics 
developed  in  the  pure  gas  system  can  be  used  in  the  study  of  plasmas 
formed  above  metallic  targets  in  argon  atmospheres.  This  section 
summarizes  the  effects  and  properties  of  the  pure  argon  plasmas  as 
detailed  in  the  previous  sections. 

Figure  1  is  a  schematic  of  the  physics  of  a  flowing  gas  laser 
supported  plasma  (ISP)  system.  The  plasma  is  initiated  by  inserting  a 
zinc  foil  or  tungsten  target  at  the  laser  focus.  The 
vaporization/thermionic  emission  of  the  targets  produces  an  electron 
density  in  the  range  of  10^  to  10^4  cm"3.  These  densities  are 
sufficient  to  lower  the  breakdown  threshold  in  the  argon  but  not  alter 
the  physics  of  the  breakdown  mechanisms.  Upon  initiation,  the  target  is 
removed  and  the  plasma  stabilizes  in  the  converging  beam.  Since  the 
laser  intensity  is  in  the  10®  to  10^  W/cm^  range,  the  breakdown  is 


characterized  as  the  "slow  combustion  wave"  or  laser  supported 
combustion  (LSC)  regime. 

In  the  LSC  regime  the  plasma  expansion  front  propagates  at  subsonic 
velocities  and  the  velocity  Is  Intensity  dependent.  The  plasma  front 
stabilizes  at  a  position  In  the  converging  laser  beam  where  the 
propagation  velocity  matches  the  applied  flowfield  velocity.  The  plasma 
expansion  is  driven  primarily  via  conduction  to  the  gas  upstream  from 
the  plasma  front.  The  front  is  weakly  absorbing  and  sufficient  laser 
energy  is  transmitted  to  the  downstream  regions  to  maintain  the  plasma 
volume. 

The  absorption  mechanism  is  Inverse  bremsstrahlung  In  the  focal 
volume  and  conduction  drives  the  plasma  volume  outside  the  laser  beam 
geometry.  There  are  two  effects  of  radiation.  First,  the  VUV  emitted 
by  the  hot  plasma  core  preheats  the  upstream  gas.  Second,  the  radiative 
conductivity  is  a  significant  percentage  of  the  effective  thermal 
conductivity  which  supports  the  plasma  outside  the  laser  geometry. 

Two  concerns  arise  In  a  laser  propulsion  scheme.  First,  the 
radiation  losses  from  the  hot  plasma  region  to  the  containment  vessel 
walls  may  be  considerable.  Factors  affecting  radiation  losses  are 
pressure,  flow  rate,  and  beam  geometry.  In  general  a  long,  narrow 
plasma  at  high  flow  rates  or  a  larger,  cooler  plasma  is  desirable.  The 
purpose  of  the  LSP  Investigation  is  to  quantify  these  effects  and  verify 
the  trends  predicted  by  a  numerical  model  [22].  The  system  shown  in 
Fig.  1  Is  a  dual  flow  system  In  which  part  of  the  impressed  flowfield  is 
directed  over  the  chamber  walls  for  convective  cooling.  The  main  goal 
of  the  LSP  work  is  to  optimize  thermal  efficiency  defined  as 
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^thermal/Einc1d°nt*  Thus  the  convective  mixing  of  the  cooler  outer  flow 
field  and  the  directly  heated  gas  Is  critical.  Thermal  efficiencies  in 
excess  of  50  percent  are  desired  if  laser  propulsion  is  to  be 
feasible.  An  additional  constraint  Is  that  the  mixing  in  the  downstream 
region  results  in  temperatures  In  excess  of  typical  combustion  product 
temperatures  (>  4,000  K).  The  gas  would  eventually  be  expanded  in  a 
supersonic  nozzle  to  convert  thermal  energy  to  kinetic  energy  or 
thrust.  A  complete  review  of  this  application  can  be  found  in  the  paper 
by  dumb  and  Krier  (20). 

2.6  Plasmas  Formed  Above  Metal  Targets 

To  accurately  model  the  breakdown  of  an  ambient  gas  and  vapor 
system  the  target  physics  of  heating,  melting,  vaporization,  and 
Ionization  must  be  accounted.  Two  extreme  cases  exist.  At  low 
intensities  the  target  surface  temperature  rises  and  conduction 
increases  Internal  temperatures,  with  no  phase  changes.  At  high 
Intensities  multiphoton  Ionization  rapidly  Initiates  the  plasma  above 
the  target.  Between  the  two  limiting  cases,  the  breakdown  physics  must 
Include  the  processes  of  change  of  phase,  vaporization  pressure, 
thermionic  emission,  and  shock  wave  formation.  To  further  burden  the 
effort,  the  temperature  dependent  optical  properties  must  also  be 
accounted. 

In  this  investigation,  the  main  concern  Is  With  the  steady  state 
plasma  volume.  The  initiation  mechanics  can  be  disregarded  as  they  have 
been  shown  to  effectively  be  a  source  of  free  electrons.  The  breakdown 
and  subsequent  plasma  formation  are  independent  of  the  initiation 
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mechanisms.  A  detailed  treatment  of  the  initiation  physics  can  be  found 
in  Ralzer  (381  for  the  CW  case  and  Hughes  [391  for  the  pulsed  case. 

Intensities  are  again  in  the  10^  to  10^  W/cm2  range  and  LSC  theory 

adequately  describes  the  plasma  Initiation  and  propagation.  At  higher 

Intensities  (>I0*°  W/cm2)  the  plasma  formation  is  analogous  to 

detonation  or  blast  wave  theory  and  several  references  can  be  found  In 

the  literature  review.  Figure  2  Is  a  schematic  of  the  physics  of  a 

plasma  formed  In  an  atmosphere  above  a  metallic  target.  When  the  laser 

beam  Is  focused  onto  the  target  rapid  melting  and  vaporization  occur. 

Surface  defect  heating  also  aids  In  the  plasma  Initiation.  The  metallic 

vapor  Is  Initially  contained  In  the  keyhole  region  which  Is  at  high 

pressures  on  the  order  of  a  few  atmospheres. 

♦ 

A  shock  wave  may  Initially  expand,  but  the  shock  is  weak  in  the  LSC 
case  and  does  not  heat  the  ambient  gas  to  temperatures  high  enough  for 
laser  absorption  to  begin.  The  shock  is  followed  by  a  slower  contact 
discontinuity  which  separates  the  initial  nonequilibrium  expansion  of 
the  metallic  vapor  as  It  "snowplows"  through  the  ambient  gas. 
Experimental  results  In  the  LSC  regime  have  Indicated  two  significant 
features.  First,  the  pressure  throughout  the  steady-state  plasma  volume 
including  the  region  at  the  keyhole/ambient  Interface  Is  nearly  constant 
at  the  ambient  value.  Second,  the  Mach  number  just  outside  the  keyhole 
region  has  been  measured  to  be  on  the  order  of  0.05.  Thus  some  effect 
must  restore  equilibrium  to  the  metallic  vapor  across  the  Interface. 
The  currently  accepted  model  was  proposed  by  Knight  (40).  A  dis¬ 
continuity  analogous  to  a  Knudsen  layer  could  be  attached  to  the  target 
surface  to  restore  translational  equilibrium  to  the  metallic  vapor.  The 
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Knudsen  layer  is  a  discontinuity  which  restores  translational  equilib¬ 
rium  to  the  high  pressure  material  vapor  as  It  expands  into  the 
ambient.  The  temperature  decrease  across  the  Knudsen  layer  would  be 
sufficient  to  maintain  the  plasma  by  supporting  a  sufficient  electron 
density  at  the  Interface.  In  this  study.  It  has  been  seen  that  the 
metallic  vapor  line  spectrum  is  evident  in  a  small  region  near  the 
target  surface.  The  metal  vapor  continuum  may  affect  the  continuum 
throughout  the  plasma  and  will  be  analyzed  in  Chapter  4. 

The  remaining  problem  concerns  the  electron  and  heavy  particle 
temperatures  near  the  target  surface.  The  metal  vapor,  as  expected,  is 
a  few  hundred  degrees  above  the  melting  point.  If  total  equilibrium 
were  to  exist,  the  electron  temperatures  would  also  be  at  a  similar 
temperature.  However,  for  the  plasma  to  exist,  the  electron  temper¬ 
atures  in  this  region  must  approach  10,000  K  for  IB  absorption  to 
occur.  This  criteria  also  forms  a  necessary  boundary  condition  for  the 
plasma  modeling.  The  analysis  of  this  region  Is  presented  In  Chapter  5. 

2.7  Energy  Transport  to  the  Ambient 

As  was  seen  in  the  literature  survey  (Chap.  1)  little  or  no  effort 
has  been  made  to  assess  the  role  of  various  heat  transport  mechanisms 
between  the  plasma  and  a  target.  In  a  laser  supported  plasma  the 
primary  mechanisms  are  conduction  and  radiation.  The  magnitudes  of  the 
effects  with  respect  to  the  transmitted  laser  flux  are  of  primary 
concern  if  an  accurate  assessment  of  the  target  flux  Is  to  be  made. 


Most  references  on  plasma  heat  transfer  concern  convective 
transport  to  a  body  In  transverse  flow.  In  this  application,  similarity 


functions  with  appropriate  reference  temperatures  and  temperature 
corrections  have  been  successfully  applied  to  plasma  transport  cal¬ 
culations.  Other  factors  affecting  the  calculations  include  non- 
equilibrium  mechanisms  such  as  dissociation  and  ionization  and  the 
inverse  processes.  Hot  particle  diffusion  to  walls,  plasma  radiation, 
and  generation  of  sheaths  are  other  considerations. 

However,  for  the  plasmas  Investigated  herein,  the  major  concern  is 
of  transport  to  confinement  vessel  walls  and/or  to  a  target.  As 
examples  of  the  magnitudes  of  the  transport  phenomena,  various  figures 
in  the  literature  reported  the  radiative  flux  from  the  plasma  varying 
between  a  few  percent  to  nearly  100  percent  of  the  total  flux  delivered 
to  the  target.  Several  arc  plasma  studies  have  attempted  to  monitor  the 
radiative  and  total  flux  for  the  regions  of  a)  chamber  walls  and 
b)  anode  (target)  surface. 

Lukens  and  Incropera  [41]  showed  that  radiation  accounts  for  as 
much  as  30  percent  of  the  wall  loading.  The  effort  was  an  attempt  to 
model  various  mechanisms  In  a  constricted  argon  arc  jet.  The  30  percent 
figure  corresponds  to  the  fully  developed  flow  region  near  the  anode. 
The  result  of  the  equilibrium  model  was  an  overestimate  of  the  radiative 
flux.  It  appeared  that  a  nonequlllbrlum  model  Incorporating  unique 
electron  and  heavy  particle  temperatures  was  more  realistic.  However, 
this  effect  may  be  an  artifact  explained  by  Ref.  [24]  In  that  the 
diffusion  approximation  overestimates  radiative  conductivity  at 
pressures  below  two  atmospheres.  This  argument  will  be  treated  in 
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Ushlo  and  Hatsuda  (42)  assessed  the  target  transport  in  an  argon 
TIG  welding  assembly.  The  TIG  process  involves  a  forced  convection  flow 
from  the  cathode  to  anode  (target)  and  the  relative  magnitudes  of 
conduction  and  convection  are  not  separable.  However,  the  combined 
effect  of  conduction  and  convection  was  calculated  to  be  8  percent  of 
the  total  arc  power.  Radiative  transport  was  not  mentioned  except  in 
the  problem  formulation  and  one  assumes  this  to  mean  that  radiation 
effects  have  little  significance. 

As  a  first  order  formulation,  from  the  results  presented  in 
Chap.  4,  an  assessment  of  the  radiative  flux  will  be  made.  The  plasma 
is  modeled  as  a  disk  in  the  range  of  2-5mm  radius,  residing  l-3mm  above 
the  target.  A  classical  radiative  view  factor  can  be  calculated  between 
the  disk  and  the  refracted  laser  spot.  Assuming  all  the  absorbed  laser 
energy  Is  lost  as  radiation,  half  of  which  is  emitted  toward  the  plasma, 
the  net  radiative  flux  can  be  shown  to  be  less  than  2  percent  of  the 
transmitted  laser  flux.  In  most  cases,  the  value  Is  less  than  1 
percent.  The  radiative  flux  is  on  the  order  of  10^  W/cm^  in  this 
region,  so  effects  outside  the  laser  geometry  should  also  be 
negligible.  Thus,  It  remains  to  assess  the  conduction  heat  transfer  as 
the  small  radial  velocities  encountered  In  this  Investigation  should 
minimize  the  effect  of  convection.  Conduction  will  be  detailed  in 

4 

Chapter  5. 


2.8  Local  Thermodynamic  Equilibrium  (LTE) 

Several  equilibrium  models  exist  for  high  temperature  ionized  gas 
systems.  All  the  models  are  based  on  the  assumption  that  all  particles 
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have  Maxwellian  velocity  distribution  although  not  necessarily  at  the 
same  characteristic  temperature.  LTE  also  assumes  that  the  system  is  of 
sufficient  density  for  collislonal  processes  to  dominate  the  rate 
equations.  This  Implies  that  all  energy  levels  obey  Boltzmann 
statistics.  Coronal  equilibrium  assumes  that  the  density  is  low  enough 
for  radiative  processes  to  dominate  the  rate  equations.  This  can  be 
interpreted  as  the  gas  temperature  being  low  enough  that  Inverse 
bremsstrahlung  absorption  does  not  occur  In  the  plasma  volume 
(<  10,000  K).  This  condition  is  obviously  Invalid  In  the  case  of  a 
laser  supported  plasma  and  applies  more  to  arc  discharge  plasmas.  The 
gray  area  between  the  two  regimes  Is  covered  by  colllslonal-radlative 
equilibrium  and  various  partial-  or  multi -temperature  local  thermo¬ 
dynamic  equilibrium  models.  These  last  models  are  more  complex  and  make 
It  desirable  to  use  either  of  the  first  two  models.  At  pressures 
approaching  atmospheric,  the  LTE  model  appears  to  be  valid.  In  most 
Investigations,  LTE  has  been  applied  via  statements  such  as  "...has  been 
shown  In  previous  studies..."  or,  "It  is  well  known  that  LTE  prevails  at 
these  conditions.".  The  validity  of  LTE  will  now  be  explored  for  the 
plasma  conditions  of  Tg  greater  than  10,000  K  and  atmospheric  pressures. 

Grlem  [431  makes  several  points  regarding  LTE.  Spatial  variations 
In  the  temperature  field  should  be  small  enough  that  diffusion  of  heavy 
particles  is  negligible  during  equilibrium  relaxation 

i  ? 

times  (-10  sec).  Griem  also  notes  that  If  LTE  is  valid,  the  electron 
and  heavy  particle  temperatures  should  not  differ  significantly.  This 
Is  due  to  the  relatively  efficient  coupling  of  energy  between  electrons 
and  heavy  particles  since  the  heavy  particle  mean  free  paths  are 
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relatively  short,  offsetting  the  effect  of  mass  ratio.  Electron  loss 
processes  such  as  radiation  actually  lower  the  restrictions  on  LTE  since 
the  energy  transfer  requirements  between  particles  have  been 
decreased.  Heavy  particle  conduction  losses  would  restrict  LTE  but  are 
usually  small  compared  to  electron  losses.  Finally,  in  the  outer  radial 
regions  where  temperature  gradients  are  maximum,  diffusion  of  hotter 
heavy  particles  from  the  core  tend  to  offset  the  gradients  and  maintain 
LTE. 

A  schematic  of  the  energy  transport  within  a  monatomic  plasma  is 
shown  in  Fig.  3  [27].  In  general,  the  electron  density  must  be  high 
enough  for  collislonal  de-excitation  to  be  roughly  ten  times  more 
probable  than  radiative  de-excitation.  For  an  optically  thin  plasma, 
considering  the  cross-sections,  the  Saha  equation  can  be  reduced  to 
143]: 


Ne  >  1.6  x  1012  T*  (4E)^ax  cm-3  (2.27) 

Where  aE  is  the  largest  energy  difference  (eV)  between  excited  states. 

For  pure  argon,  the  first  excited  states  are  approximately  70  percent  of 

the  ionization  energy  hence  E^^  is  in  the  range  of  3  to  5  eV.  This 

16  -3 

range  corresponds  to  an  electron  number  density  of  10  cm~J  which 
appears  a  reasonable  value  at  the  prescribed  conditions. 

To  investigate  the  temperature  difference  between  the  electrons  and 
heavy  particles  one  looks  at  the  energy  equation: 

“UCP  37  '  b<V<Te  -V 


(2.28) 
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where  b  is  a  func  jii  of  the  electron  number  density  and  plasma 
frequency,  hence  f(Te).  The  equation  assumes  radial  conduction  is 
negligible  and  has  been  solved  numerically  by  Dresvln  (37). 
At  du  *  .5  g/cm3-sec  and  constant  Tg(=13,000  K),  Tg  *  Tg  within  the 
first  4mm  of  cool  gas  flow  Into  the  plasma  volume.  If  the  energy 
equation  is  solved  assuming  the  axial  (z)  temperature  gradients  are 
zero,  one  finds  that  the  electron  and  heavy  particle  temperatures  are 
within  a  few  percent  over  the  central  2-5cm  radii.  Thus,  at 
temperatures  in  excess  of  10,000  K,  the  equivalence  should  extend  to 
larger  spatial  dimensions. 

Eddy  1 44 1  published  a  classical  review  of  equilibrium  models  and 
proposed  the  concept  of  multlthermal  local  thermodynamic  equilibrium 
(MLTE).  In  the  MLTE  analysis,  the  various  excitation,  electron,  and 
heavy  particle  temperatures  are  accounted  for.  Eddy  also  analyzes  the 
probable  error  of  the  various  spectroscopic  diagnostic  techniques  where 
LTE  Is  assumed.  The  analysis  assumes  a  10  percent  measurement  error. 
The  absolute  line  Intensity  method  appears  to  be  the  minimal  error  at  .3 
to  3  percent  at  Tg  =>  20,000  K.  The  off-axis  peak  method  also  has  a  3 
percent  error.  The  ratio  of  Ionic  to  atomic  argon  lines  resulted  in  a 
5 %  error  in  temperature.  Thus,  the  ion- ion  ratio  should  be  near  this 
magnitude  as  the  energy  differences  are  similar  and  the  ion-ion  ratio 
does  not  depend  on  validity  of  the  Boltzmann  relations  as  much  as  the 
ion-atom  technique.  Line  to  continuum  technique  errors  at  similar 
conditions  appear  to  be  10  percent.  It  will  be  shown  in  this 
Investigations  experimental  results  that  a  ±15  percent  error  in  the 
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calculated  temperature  field  results  in  a  less  than  ±5  percent  deviation 
in  the  global  properties  of  absorption  and  radiation  loss. 

Farmer  and  Haddad  [45]  attempted  to  experimentally  verify  LTE  in  an 
argon  arc  plasma.  The  theory  Is  that  the  measured  emission  coefficient 
(W/cm^-sr)  must  be  identical  for  regions  of  identical  temperatures.  The 
baseline  temperature  was  the  off-axis  peak  or  "norm"  temperature 
location.  The  authors  concluded  that  an  arc  plasma  at  a  pressure  of  1.4 
atm  demonstrated  LTE  behavior  (±5  percent)  over  the  entire  arc  column. 
However,  at  1  atm  pressure,  the  error  is  also  ±5  percent  except  in  the 
region  near  the  cathode  (±20%).  This  is  possibly  due  to  electron  charge 
repulsion  In  the  area  near  the  cathode  and  LTE  again  appears  to  be  valid 
for  an  atmospheric  pressure  laser  supported  plasma. 

Finally,  as  a  first  order  calculation,  theory  Is  combined  with 
experiment  in  order  to  derive  a  function  proportional  to  the  electron 
energy  distribution.  From  gas  kinetic  theory,  the  thermal  conductivity 
can  be  expressed  as  in  the  form  of  Eqn.  (2.23).  Equation  (2.23)  was 
derived  for  <q>  independent  of  temperature.  A  rigorous  derivation  of 
thermal  conductivity  makes  comparison  between  experiment  and  theory 
difficult.  Thus,  as  a  first  order  approximation,  <q>  will  be  allowed  to 
vary  with  temperature  In  Eqn.  (2.23).  It  should  be  noted  that  the  large 
range  of  calculated  values  of  <q>  make  the  variation  with  electron 
temperature  assumption  a  small  error  since  a  discrete  analysis  could  be 
performed  by  adjusting  the  constant  term  of  Eqn.  (2.23). 

If  one  applies  the  conductivity  data  of  Emmons  [46]  for  argon  at 
latm,  4,000-14,000  K,  one  can  calculate  <q(Tg)>  by  inverting 
Eqn.  (2.23).  Using  the  measured  data  of  Golden  and  Bandel  [47]  for  the 
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total  momentum  transfer  cross-section,  q(E),  it  is  possible  to  derive  a 
function  proportional  to  the  "experimental"  distribution  function,  f(E). 

<q(E)>  a  f  q(E)f(E)dE  (2.29a) 

and 

f (E)  a  (d  <q(E)>/dE]/q(E)  (2.29b) 

The  data  and  theoretical  profiles  were  fit  to  fifth  order  polynomials 
via  least  squares.  The  resulting  distribution  function  (not  normalized) 
is  shown  in  Fig.  4. 

One  notes  the  two-temperature  appearance  of  the  distribution  with 
the  break  occurring  at  1.05eV  (~  11,000  K)  This  is  due  to  the  radiative 
conductivity,  as  discussed  previously,  in  the  data  of  Emmons  [47j  at 
temperatures  above  10,000  K.  If  one  uses  theoretical  values  of  thermal 
conductivity  (no  radiation  effects)  the  result  is  the  dashed  line  in 
Fig.  4.  Thus  the  piece-wise  linear  (logarithmic)  behavior  of  the 
distribution  further  supports  the  assumption  of  a  Maxwellian  velocity 
distribution.  Since  the  effort  is  concerned  with  temperatures  above 
10,000  K,  It  appears  that  a  single  characteristic  electron  temperature 
will  describe  the  system  adequately. 


Non-normal i zed  Calculated 
Electron  Energy  Distribution 
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3.  EXPERIMENTAL  APPARATUS 

Two  plasma  chambers  have  been  developed  for  this  research  effort. 
The  first  chamber  is  used  in  the  pure  argon  studies  and  was  designed  to 
Incorporate  multiple  diagnostics.  The  plasma  initiation  and  flow 
chamber  (PIFC)  includes  thermocouple,  calorimeter,  pressure,  and  optical 
diagnostic  ports.  The  PIFC  was  a  collaborative  design  effort  for  the 
previously  mentioned  AFOSR  project.  The  details  of  the  design  can  be 
found  In  the  theses  of  J.  Bender,  M.S.  1985  (481  and  R.  Glumb,  Ph.D. 
1986  [49].  The  second  apparatus,  the  plasma  Initiation  chamber  (PIC), 
Is  used  in  the  study  of  metallic  plasmas  In  controlled  atmospheres.  PIC 
incorporates  a  translation  stage  to  continually  refresh  the  target 
material  as  well  as  pressure  and  optical  diagnostic  ports.  This  chapter 
will  detail  the  two  plasma  chamber  features,  describe  the  diagnostic 
apparatus,  and  pertinent  data  reduction  techniques. 

3.1  Plasma  Initiation  and  Flow  Chamber  (PIFC) 

A  schematic  of  the  PIFC,  support  stand,  and  beam  steering  optics  is 
shown  in  Fig.  5.  The  stand  Is  constructed  of  heavy  steel  uni  strut 
framing  and  Is  mounted  on  four  pneumatic  tires.  The  system  is 
positioned  Into  the  laser  beam  path  and  can  be  located  via  alignment 
pins  Into  holes  In  the  floor. 

The  laser  is  a  10  KW  CW  Avco  Everett  HPL-10  system.  The  laser  is 
an  unstable  oscillator  with  an  annular  beam  similar  to  a  TEMqj*  beam 
profile.  The  beam  Is  focused  to  pass  through  an  aero-window  and 
subsequently  relmaged  to  an  approximate  2.6  in.  diameter  at  the  outer 
1/e  location  (1.98  In.  10).  The  beam  is  delivered  horizontally  across 
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the  laboratory  to  the  flat  copper  mirror  shown  In  the  upper  left-hand 
corner  of  Fig.  5.  The  beam  Is  reflected  90  degrees  downward  onto  a 
second  copper  flat  which  reflects  the  beam  upward  Into  a  6.0  In. 
diameter  convex  copper  mirror  (R  =  12  in.).  The  convex  mirror  expands 
the  bean  Into  a  20  in.  diameter  concave  primary  focusing  mirrors  (R  = 
31  In.). 

The  arrangement  permits  the  beam  to  be  focused  vertically  upward 
into  the  PIFC  which  minimizes  the  effects  of  buoyancy  In  plasma 
stability  and  diagnostics.  The  system  also  permits  a  variable  f/# 
(defined  from  the  primary  mirror  to  focus)  of  f/2  to  f/3.4  which  was  a 
parameter  of  Interest  at  the  onset  of  the  project.  The  f/#  Is 
changeable  by  relative  translation  and  rotation  of  the  lower  flat, 
convex,  and  concave  mirrors  as  shown  In  Fig.  6.  The  low  f/#s  provide 
for  a  large  beam  cross-section,  minimizing  local  intensity,  as  the  beam 
enters  the  PIFC  through  a  NaCl  window,  minimizing  window  damage.  All 
the  mirrors  are  water  cooled  to  prevent  thermal  damage  except  the  20  in. 
mirror  which  has  sufficient  mass  to  conduct  energy  away  from  the  mirror 
surface. 

The  optical  system  Is  off-axis  and  depending  on  the  f/#t  a 
significant  amount  of  aberration  may  be  present  In  the  focal  volume. 
Figure  7  shows  two  orthogonal  Intensity  mappings  of  the  annular  focal 
volume.  The  data  was  acquired  using  a  laser  beam  analyzer.  Model  LBA- 
IA,  from  ALL  GmbH,  West  Germany.  The  Instrument  passes  a  molebdynum 
needle  through  the  laser  beam  and  the  reflected  energy  Is  detected  by 
two  orthogonal  photodiodes.  The  output  Is  stored  on  a  digital  storage 
oscilloscope  where  accurate  time-base  measurements  of  the  Intensity 
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profiles  can  be  made.  The  time-base  is  converted  to  spatial  dimensions 
using  the  calibration  figures  supplied  with  the  LBA. 

The  solid  lines  In  Fig.  7  are  a  third  order  curve  fit  of  the 
Indicated  data  points  which  are  the  1/e  loci  of  the  annular  focal 
volume.  One  sees  the  effect  of  the  off-axis  optical  system  in  Fig.  7. 
The  one  trace  Is  generally  wider  than  the  other  on  the  positive  side  of 
the  focus.  The  trend  Is  reversed  on  the  negative  side  of  the  focus. 
This  implies  that  the  beam  Is  elliptical  and  astigmatism  is  present. 
However,  no  tangential  or  sagital  focus  could  be  found  in  the  IR  or 
visible  (HeNe)  focal  volumes.  Unfortunately  the  detector  saturates  at 
the  Intensities  encountered  In  the  focus  region,  resulting  in  a  narrow 
trace.  The  best  estimate  of  the  circle  of  least  confusion  Is  two 
millimeters  based  on  the  behavior  in  the  outer  thirds  of  the  traces. 
This  figure  agrees  fairly  well  with  the  HeNe  tracer  laser  estimates  but 
further  attempts  to  prevent  LBA  saturation  should  be  made. 

The  PIFC  Is  a  5.0  In.  ID  stainless  steel  cylinder.  Stainless  steel 
was  chosen  for  corrosion  resistance,  relatively  high  melting  point,  and 
strength.  The  large  diameter  was  chosen  to  minimize  radiative  heating 
of  the  walls  and  minimize  wall  effects  in  the  downstream  mixing 
region.  This  type  of  design  is  typically  called  a  dual  flow  chamber. 
Cooler  outer  gas  convectively  cools  the  walls  and  subsequently  mixes 
with  the  directly  heated  gas  from  the  plasma  volume.  Figure  8  Is  a 
schematic  of  the  various  components  of  the  PIFC. 

At  the  bottom  of  the  PIFC  Is  a  flange  which  holds  the  NaCl  inlet 
windows  of  varying  thickness,  depending  on  chamber  pressure.  Inlet  gas 
is  brought  in  radially  via  a  series  of  small  holes  on  the  inside  of  the 
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(A)  Laser  Inlet  Window  Assembly 

(B)  Gas  Inlet  Assembly 

(C)  One-IIalf  Inch  Inlet  Pipes  (2) 

(D)  Main  Chamber  Section 

(E)  Thermocouple  Wire  Retract  Box 
(P)  Viewing  Windows  (2) 

(G)  Target  /  Solenoid  Assembly 

(H)  pressure  Transducer  and  Popoff  Valve 

(I)  One-Half  Inch  Exhaust  Pipes  (4) 

(J)  Thermocouple  Drive  Motor 

(K)  Calorimeter 

(L)  Laser  Inlet  Window 
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Figure  8  PIFC  Components  (49] 
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PIFC.  Some  of  the  gas  is  diverted  onto  the  window  for  cooling,  part  is 
directed  normal  to  the  vertical,  and  the  remainder  directed 
downstream. 

Above  the  inlet  assembly  is  the  plasma  containment  region. 
Diagnostic  windows  provide  direct  optical  viewing  of  the  plasma  event. 
The  thermocouple  box  contains  the  excess  wire  needed  to  permit 
thermocouple  movement  along  the  vertical  via  a  carriage/slider 
assembly.  Pressure  transducers,  pressure  relief  valves,  stationary 
thermocouples,  as  well  as  other  ports,  are  also  located  in  this 
region.  Above  the  main  chamber  the  gas  is  exhausted  through  four  ports, 
recombined  in  a  manifold  and  passed  through  a  tube  in  shell  heat 
exchanger  before  being  exhausted  to  the  ambient.  At  the  top  of  the 
chamber  a  water  cooled  copper  cone  calorimeter  is  mounted.  The 
calorimeter  is  used  to  measure  the  laser  power  transmitted  by  the  plasma 
and  the  global  plasma  absorption  can  be  determined.  When  not  being  used 
as  a  diagnostic  tool  the  calorimeter  serves  as  a  beam  dump.  Particular 
details  of  the  design  and  assembly  can  be  found  in  the  theses  by  Bender 
(48)  and  Glumb  [49]. 

The  experiment  begins  when  the  laser  beam  is  delivered  to  the 
mirror  system  and  focused  Into  the  chamber.  Peak  Intensities  of  the 
laser/optical  system  are  on  the  order  of  10®  W/cm2  which  are 
Insufficient  to  breakdown  the  flowing  argon  without  an  initial  source  of 
free  electrons.  The  electrons  are  created  by  inserting  a  zinc  foil  or 
more  recently  a  reuseable  tungsten  rod,  Into  the  focal  region  via  an  air 
cy Under/solenoid  assembly.  Target  vaporization  Is  sufficient  to  create 
a  high  enough  electron  density  for  cascade  breakdown  in  the  argon  to 
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occur.  The  plasma  then  stabilizes  in  the  converging  beam  as  discussed 
in  Chapter  2  and  the  various  diagnostic  techniques  can  be  applied. 
Figure  9  Is  a  double  exposure  of  the  chamber  and  plasma  during  a  typical 
experiment. 

3.2  Plasma  Initiation  Chamber  (PIC) 

To  facilitate  the  study  of  plasmas  formed  above  metallic  targets,  a 
second  chamber  (PIC)  was  designed.  Figure  10  is  a  drawing  of  the  PIC 
features.  The  system  Is  a  rectangular  aluminum  box  (Inside  dimensions 
of  14  In.  by  5  in.  by  5  in.)  with  0.5  in.  thick  walls.  Aluminum  was 
chosen  for  machlnabillty  as  wall  degradation  was  found  to  be  minimal  in 
the  PIFC  experiments.  One-half  inch  thick  walls  were  chosen  to  add  mass 
to  the  unit  In  consideration  of  heat  transfer  effects.  The  thick  walls 
also  permit  absolute  pressures  up  to  3  atmospheres  to  be  maintained  with 
a  factor  of  safety  of  two,  assuming  proper  window  design  and 
Installation. 

The  beam  is  delivered  to  the  target  via  a  rectangular  or  circular 
NaCl  window  In  the  top  of  the  chamber.  All  windows  are  flange  mounted 
with  Inner  and  outer  0-r1ng$  for  sealing  and  cushioning.  Since  the 
laser  focus  rapidly  begins  drilling  the  target  an  Internal  translation 
stage  was  Incorporated.  The  stage  continually  refreshes  the  target 
material  minimizing  drilling  effects  and  maintaining  the  laser  target 
Interaction  as  a  surface  event.  The  translation  speed  can  be  varied  up 
to  the  point  where  the  plasma  deflects  out  of  the  laser  beam  geometry. 
This  condition  makes  the  diagnostic  techniques  invalid  from  symmetry 
considerations.  Translation  speeds  up  to  60  in. /min  were  found  to  have 
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no  effect  on  plasma  deflection  in  this  study.  The  translation  param¬ 
eters  are  controlled  via  a  Vellmex®  control  unit  configured  as  a  dumb 
terminal  and  connected  to  a  host  terminal.  The  system  uses  a  form  of 
BASIC  language  and  up  to  four  motors  can  be  controlled  independently. 
This  system  will  be  discussed  further  in  Section  3.3. 

The  stage  motor  is  mounted  externally  and  is  sealed  with  an  0-ring 
and  flange  apparatus  to  prevent  Infiltration.  The  stage  Is  12  in.  long 
with  a  net  translation  distance  of  8  in.  end-to-end.  In  general  the 
translation  speed  was  set  to  a  slow  value  to  permit  multiple  data 
acquisition  during  a  single  translation.  Multiple  data  acquisition 
ensures  data  Integrity  and  typically  two  to  four  scans  at  various 
elevations  above  the  target  could  be  performed  during  a  single  trans¬ 
lation.  After  a  given  Internal  translation  of  the  PIC  stage,  the  PIC  Is 
translated  laterally  such  that  subsequent  Internal  translations  are  made 
on  target  material  which  has  not  undergone  laser  damage. 

The  chamber  lid  Is  mounted  via  socket  screws  and  a  teflon  gasket  to 
prevent  leaks  and  Infiltration.  Various  0.25  In.  pipe  thread  ports  are 
available  on  the  top  and  various  positions  about  the  chamber.  The  ports 
typically  have  pressure  transducers  and  inlet  and  exhaust  fittings 
attached.  The  gas  Is  fed  Into  the  chamber  In  a  slow  flow  such  that  the 
plasma  is  not  disturbed  by  any  local  gas  dynamic  effects  as  might  be 
found  In  various  gas  assist  nozzle  systems.  The  pressure  is  monitored 
to  Insure  that  Internal  pressure  Is  slightly  above  atmospheric,  again  to 
minimize  Infiltration.  The  side  walls  have  opposing  3  In.  by  4  In. 
window  ports  to  facilitate  emission,  transmission,  or  interferometric 
diagnostics.  An  end  window  also  permits  simultaneous  application  of 
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multiple  diagnostics.  Figure  11  is  a  double  exposure  photograph  of  the 
PIC  and  plasma  during  an  experiment.  A  major  concern  is  damage  to  the 
NaCl  window  by  the  plasma  and  reflected  laser  beam.  The  5  in.  clearance 
between  the  window  and  target  was  a  compromise  between  the  above  damage 
effect  and  the  desire  to  use  the  PIC  with  a  variety  of  laser  systems. 
The  window  system  has  been  redesigned  with  an  internal  gas  jet  applied 
to  the  window  for  cooling  and  to  prevent  oxide  deposition.  The  system 
has  recently  been  used  to  reduce  window  damage. 

Figure  12  is  the  LBA  intensity  traces  for  the  f/7  Cassegrain 
telescope  focal  volume.  The  beam  is  relatively  symmetrical.  The  spot 
size  is  approximately  1.0  mm  (average  between  the  two  traces)  and  the 
intensity  profile  appeared  to  be  near-Gaussian  at  the  focus.  The 
results  agree  with  the  HeNe  tracer  system  which  also  shows  symmetry. 
The  spectroscopic  data  demonstrated  symmetrical  behavior,  further 
supporting  the  assumption  of  cylindrical  symmetry. 

3.3  Spectroscopy 

The  temperatures  in  excess  of  10,000  K  encountered  in  the  plasma 
core  render  conventional  temperature  diagnostics  Ineffective.  Spectro¬ 
scopy,  as  adapted  from  astrophysics,  is  the  primary  tool  used  in  the 
study  of  laboratory  plasmas.  Spectral  emission  and  absorption  behavior 
can  be  used  to  map  local  properties  such  as  temperature  and  number 
density  (pressure).  The  main  assumption  is  local  thermodynamic 
equilibrium  as  discussed  in  Chapter  2,  without  which  the  relationships 
between  emission  behavior  and  plasma  parameters  are  unwieldy.  Since  the 
spectral  behavior  is  dominated  by  electron  interactions,  the  calculated 
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parameters  correspond  to  the  electrons.  The  heavy  particle  temperatures 
must  be  determined  independently  with  no  reliable  diagnostic  techniques 
currently  available.  Instead,  the  heavy  particle  parameters  are 
typically  calculated  from  the  measured  electron  distributions  as  will  be 
discussed  in  Chapter  5.  Fortunately  the  electron  parameters  dominate 
the  laser  plasma  Interaction  physics  associated  with  this  investigation. 

3.3.1  Review  of  the  Techniques 

Spectroscopy  can  be  divided  into  two  general  categories,  ab¬ 
solute  and  relative.  Each  type  has  been  used  successfully  In  the  diag¬ 
nostics  of  laboratory  plasmas.  Absolute  techniques  typically  Involve 
the  measurement  of  a  narrow  spectral  continuum  band.  The  absolute 
emissive  power  (W/cm3-sr)  Is  calculated  and  compared  to  a  calibrated 
source  such  as  a  tungsten  ribbon  lamp.  By  matching  the  emissive  powers, 
the  local  electron  temperature  can  be  determined. 

Henricksen,  et  al.  [50)  and  Carlhoff,  et  al.  [13)  have  successfully 
demonstrated  absolute  measurements  of  laser  supported  plasmas  in  pure 
gases.  Henricksen  measured  a  10  angstrom  band  of  Ar  continuum  centered 
at  626.5  nanometers.  The  region  Is  Isolated  from  line  spectra  to 
insure  negligible  self-absorption.  Carlhoff  similarly  measured  the 
Intensity  of  a  continuum  band  near  587.6  nanometers  In  helium.  Both 
investigations  reported  peak  temperatures  approaching  15,000  K.  The 
equipment  needed  for  absolute  measurements  is  relatively  modest  but 
inflexible  due  to  the  need  of  optical  filters  centered  about  the  band  of 
Interest.  The  need  to  accurately  account  for  the  filter  transmission 
function,  calibration  source  fluctuation,  and  the  optical  transfer 
function  in  the  above  investigations  can  lead  to  significant  error. 
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The  other  area  of  spectral  diagnostics  of  interest  are  the  relative 
techniques.  Two  spectral  events,  line  and/or  continuum,  as  mapped  from 
the  same  elemental  plasma  volume,  are  used  to  determine  local  electron 
parameters.  Key,  et  al.  [51]  and  Mills  [52]  have  used  the  relative 
intensity  of  two  spectral  lines  In  argon  arc  plasmas.  Both  studies  used 
the  457.9  and  458.0  nanometer  Aril  lines.  The  upper  energy  level 
separation  of  the  lines  allowed  the  determination  of  temperatures 
approaching  15,000  K.  The  general  rule  Is  that  the  upper  level  energies 
of  the  line  events  should  differ  on  the  order  of  the  electron 
temperature  to  be  determined,  l.e.,  in  excess  of  9,000  K.  The  9,000  K 
figure  is  an  arbitrary  value. 

As  will  be  shown  In  Section  3.3.5  of  this  chapter,  the  line  ratios 
are  typically  on  the  order  of  unity.  The  nature  of  the  line  intensity 
equation  Is  that  small  fluctuations  in  the  line  ratio  result  In  small 
fluctuations  about  the  9,000  K  energy  difference.  Thus  the  energy 
difference  should  be  on  the  order  of  the  expected  temperature  field  to 
minimize  error.  The  lines  mentioned  above  meet  this  requirement.  The 
use  of  ionic  argon  lines  In  this  temperature  range  at  atmospheric 
pressures  allow  accurate  measurements  along  the  arc  plasma  centerline 
where  the  current  density  Is  a  maximum.  The  remainder  of  the 
temperature  field  Is  extrapolated  from  the  near-centerline  behavior  of 
the  neighboring  atomic  line  structure.  The  two  studies  above  appear  to 
be  the  best  data  to  date. 


Tsao  and  Pavel  1c  [53]  utilized  the  415.8  and  425.9  nanometer  atomic 
argon  lines  to  determine  electron  temperature  in  an  arc  plasma.  The 
peak  temperatures  below  10,000  K  are  expected  as  the  energy  difference 
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between  the  two  lines  is  relatively  small  ( 2,000  to  3,000  K)  leading  to 
an  underestimation  in  temperature.  The  small  energy  difference  implies 
errors  approaching  50  percent  [441  for  temperatures  which  actually 
exceed  10,000  K.  Kobayashi  and  Suga  [ 54 J  utilized  the  430.01  nanometer 
Arl  and  434.8  nm  Aril  line  pair  to  determine  the  temperature  in  an  arc 
plasma.  Peak  temperatures  of  18,000  K  were  reported  which  appear 
high.  It  was  not  obvious  whether  the  authors  accounted  for  the  more 
explicit  dependence  of  subsequent  ionization  stage  relationships  on 
LTE.  The  lowering  of  the  ionization  potential  also  may  have  had  an 
effect. 

The  ratio  of  a  line  to  continuum  event  or  continuum  to  continuum 
event  can  also  be  used.  Generally,  continuum  to  continuum  errors  are 
large  even  in  the  presence  of  LTE  [33,44].  Line  to  continuum  ratios 
have  been  used  successfully  at  the  University  of  Tennessee  by  Henricksen 
[ 55 ] .  The  technique  was  later  abandoned  as  the  hardware  and  software 
was  developed  for  absolute  measurements.  Predicted  errors  for  line  to 
continuum  diagnostics  are  slightly  higher  than  the  same  as  relative  line 
techniques,  +10  percent  versus  +3  percent,  respectively  (44).  Also 
under  the  pressure  conditions  present  In  this  investigation  LTE  should 
prevail  [551,  minimizing  the  difference  In  error  between  the  various 
techniques. 

Other  spectroscopic  diagnostics  generally  Involve  the  measurement 
of  line  broadening.  The  broadening  mechanisms  compete  and  make 
resolution  of  Individual  contributions  difficult.  Stark  (pressure)  and 
Doppler  broadening  dominate  the  plasma  behavior  at  the  temperatures  and 
pressures  encountered  In  this  Investigation.  Most  authors  agree  that 
Stark  broadening  dominates  once  the  Ionization  exceeds  a  few  percent. 


Stark  broadening  is  due  to  the  forces  exerted  on  the  electrons  by 
the  neighboring  heavy  particles.  The  effect  is  a  quadratic  (sometimes 
linear)  broadening  of  the  emission  line  profile.  The  measured  line 
width  at  half  maximum  is  proportional  to  the  electron  number  density. 
The  technique  has  been  used  successfully  in  systems  where  the  hydrogen 
Balmer  series  Hg  was  resolvable.  Stark  half  widths  could  be  Incorpor¬ 
ated  in  future  work  to  measure  the  electron  number  density  as  a 
continuity  check  on  temperature  measurements.  This  study  was  unable  to 
resolve  the  Balmer  Hg  line  from  the  water  content  in  the  argon.  It 
appears  that  a  small  percentage  of  hydrogen  would  need  to  be  Introduced 
to  the  flow  chambers  in  future  work. 

Doppler  broadening  results  In  a  Gaussian  line  profile  which  is 
difficult  to  resolve  from  the  Stark  effect.  This  Is  unfortunate  as  the 
Doppler  half  width  Is  proportional  to  the  heavy  particle  temperature  and 
is  the  only  diagnostic  available  for  this  purpose.  Griem  [57]  presents 
a  method  of  resolving  Stark  and  Doppler  effects  via  Voigt  profiles. 
However,  the  non-ideal  laboratory  conditions  and  associated  errors  due 
to  Instrumental  broadening  made  measurement  of  line  widths  to  +0.25  A  a 
difficult  task  with  the  current  apparatus  and  laboratory  conditions. 

To  Implement  any  line  broadening  diagnostic  technique  several 
precautions  must  be  taken.  First,  broadening  measurements  are  typically 
made  on  relatively  low  Intensity  glow  discharge  type  events  [43]. 
Extreme  care  should  be  taken  In  reducing  the  laser  plasma  Intensity 
delivered  to  the  detection  equipment,  more  so  than  In  the  spectral 
mappings  of  this  Investigation.  Background  signals  should  be  minimized 
which  was  not  possible  In  all  areas  of  this  investigation.  Secondly,  a 


77 


means  of  accurately  assessing  the  target  responslvlty  and  blooming 
characteristics  using  a  resolution  target  standard  would  have  to  be 
implemented.  Finally  the  optical  transfer  function  of  the  imaging  and 
filtering  system  would  have  to  be  determined  as  well  as  the  spectrometer 
characteristics.  Only  when  each  of  the  above  areas  has  been  quantified 
can  one  begin  to  achieve  the  accuracy  and  resolution  required  for  line 
broadening  measurements. 


3.3.2  Relative  Diagnostics 

Relative  spectroscopic  diagnostics  have  been  utilized  ex¬ 
clusively  in  this  Investigation  for  several  reasons.  First,  the 
necessity  of  a  calibration  source  in  absolute  methods  requires  a  second 
set  of  measurements.  The  calibration  run  should  be  performed  after  each 
data  set  under  identical  optical  conditions.  The  large  number  of 
experimental  conditions  to  be  Investigated  made  the  calibration 
unnecessarily  costly.  The  successful  studies  using  absolute  techniques 
for  two-dimensional  plasma  mappings  [13, SO]  have  invested  a  considerable 
amount  of  time  and  money  into  specialized  hardware  and  software 
systems.  The  systems  far  out  cost  the  system  required  for  relative 
diagnostics.  The  main  advantage  of  absolute  diagnostics  is  that  snap¬ 
shot  Images  of  the  entire  plasma  volume  can  be  made.  This  makes  for 
very  smooth  and  uniform  profiles  In  the  reduced  data.  Relative  diag¬ 
nostics  are  more  flexible  In  that  multiple  spectral  events  can  be  mapped 
in  a  one-dimensional  plasma  cross  section.  The  trade  off  between 
absolute  and  relative  diagnostics  Is  a  second  spatial  dimension  versus  a 
spectral  dimension,  respectively.  It  was  decided  that  the  spectral 
dimension  was  more  Important  In  the  assessment  of  LTE. 
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Figures  13  and  14  show  the  strong  atomic  argon  spectra  resolvable 
in  the  pure  argon  plasmas.  Since  the  upper  energy  level  differences  in 
atomic  argon  are  small,  relative  line  techniques  are  not  accurate.  Thus 
relative  line  to  continuum  techniques  have  been  used  In  this  study.  The 
lines  selected  should  be  optically  thin  to  minimize  the  effects  of  self¬ 
absorption.  As  mentioned  previously,  it  has  been  shown  that  self¬ 
absorption  even  of  the  metastables  is  minimal  in  the  400  to  500  nm 
region  of  argon. 

In  the  metallic  plasma,  the  presence  of  aluminum  vapor  appears  to 
push  the  system  closer  to  LTE  conditions.  This  Is  evidenced  In  Fig.  15 
where  the  ionic  argon  lines  have  begun  to  appear.  The  bottom  (solid) 
trace  Is  the  pure  argon  spectra  of  Fig.  14.  The  dashed  curve  is  the 
aluminum-argon  spectra.  One  notes  the  Aril  lines  at  428,  437, and  437.9 
nm  which  are  not  observable  In  the  pure  argon  results.  However,  the 
upper  energy  level  differences  of  these  lines  are  on  the  order  of 
30,000  K  and  thus  cannot  be  used  to  determine  temperature  in  this 
region. 

The  only  aluminum  structure  observable  was  the  strong  resonant  pair 
at  394.4  and  396.1  nm.  The  aluminum  lines  were  observable  only  within  a 
few  mm  of  the  target  surface  Indicating  that  the  plasma  behavior  is 
dominated  by  the  argon  species.  This  Is  consistent  with  the  results  of 
Ref.  12  wherein  the  plasma  was  found  to  have  little  metal  vapor.  The 
effect  was  attributed  to  the  attenuated  laser  energy  being  insufficient 
to  support  vaporization  although  considerable  melting  was  occurring. 
Also,  In  Ref.  14  It  was  noted  that  the  metallic  number  density  in  the 
vacuum  plasma  was  negligible  after  1  us.  The  authors  concluded  that 
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similar  results  would  occur  In  non- vacuum  atmospheres  and  plasma 
behavior  would  be  dominated  by  the  ambient  species. 

The-  line  radiation  resulting  from  bound-bound  transitions  is 
expressed  via  quantum  physics  as: 

Amn  h  c  9m  N 

h  *  ”.X-Ze"  <3-Z> 

in  W/cm3-sr  where  Amn  is  the  transition  probability  (1/s),  h  Is  Planck's 
constant,  gm  the  upper  level  degeneracy,  Em  is  the  upper  level  energy, 
Ze  the  electronic  portion  function,  and  N  the  atomic  number  density. 
Equation  (3.2)  Is  the  result  of  quantum  theory  for  radiative  emission 
and  absorption  of  a  system  as  an  electron  transfers  between  energy 
levels. 

Several  authors  have  found  the  Kramers-Unsold  theory  adequately 
describes  the  argon  continuum  below  20,000  K  [581 .  The  continuum 
emission  consisting  of  free-free  and  free-bound  transitions  summed  over 
Internal  energy  levels  yields 

I-  *  5.41  x  10"46  c ax  g  N?/(x2  tJ)  (3.1) 

tx  e  e 

In  W/cm^-sr.  C  Is  the  speed  of  light,  g  the  average  Gaunt  factor,  Ne 
the  electron  number  density,  x  the  center  wavelength  of  the 
Interval  ax,  and  T0  the  electron  temperature.  Equation  (3.1)  Is  a  form 
of  the  Kramers-Unsold  continuum  theory  [29,58).  Equation  (3.1)  was 
derived  by  comparing  classical  and  quantum  radiation  theory.  At  long 
wavelengths  In  the  Infrared  quantum  and  classical  theory  agree.  Above 
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some  frequency,  vQ,  the  spectral  nature  of  quantum  theory  differs 
significantly  from  classical  theory.  By  mapping  the  two  theories  in 
this  region,  it  was  shown  that  the  energy  of  the  quantum  line  events  is 
equal  to  the  energy  of  the  surrounding  spectrum  as  predicted  by 
classical  theory.  Thus  classical  theory  was  shown  to  accurately  predict 
continuum  radiation.  Classical  theory  Is  particularly  valid  for 
integrated  measurements  over  a  large  wavelength  interval. 

Using  the  value  of  g  -  2.2  from  Oettlnger  159),  gm  and  Amn  for  the 
415.8  nm  ARI  line  from  the  NBS  Reference  Book  {601  and  Ze  as  derived  by 
Tan  161],  Eqns.  (3.1)  and  (3.2)  are  combined  to  yield 

Ix/Ic  =  1.3357  x  1020(HT^/n|)  ax  (3.3) 


Assuming  LTE  at  one  atmosphere,  N  and  Ne  are  functions  of  Tg  only  and 
Eqn.  (3.3)  can  be  reduced  to  a  function  of  Tg.  Figure  16  is  the 
ratio  Ix/Ic  as  a  function  of  Tg  for  a  continuum  band  (ax)  of  10  A  for 
the  415.86,  420.0,  and  696.5  nm  Arl  lines.  The  ratio  rapidly  decreases 
as  Tg  approaches  17,000  K  with  a  minimum  at  approximately  17,500  K  when 
Aril  contributions  to  the  number  density  become  significant. 

In  the  calculation  of  temperature  from  relative  line  Intensities, 
Eqn.  (3.2)  can  again  be  applied  to  each  line  resulting  in 


A1  gl  x2 
A2  ^2  ^1 


exp  ( 


(3.4) 


which  can  be  rearranged: 


Figure  16  Line  to  Continuum  Ratio 


2/1 


AD-A174  579  THE  ROLE  OF  THE  PLASMA  DURING  LASER-GAS  LASER-METAL  ' 
INTERACTIONS  LU)  ILLINOIS  UNIV  AT  URBANA  DEPT  OF 
MECHANICAL  AND  INDUSTRIAL  ENG  T  J  ROCKSTROH  ET  AL 
UNCLASSIFIED  12  OCT  86  UILU-ENG-86-4011  AFOSR-TR-86-1852  F/C  28/5 


NL 


85 


Te  “  <E2  -  (Ix1A292x1)/(Ix2A1g1^2)J  <3*5> 

As  mentioned  earlier,  the  energy  difference,  E2  -  Ej,  must  be  on  the 
order  of  Te  '  .•  ,  0.9  eV  or  10,000  K)  as  the  denominator  of  Eqn.  (3.5) 
typically  does  not  vary  much  about  unity.  The  values  of  vary  from 
author  to  author,  however  the  variation  from  similar  sources  should  be 
offsetting  In  Eqn.  (3.4).  Figure  17  shows  the  ratio  of  the  457.9  and 
458.9  nm  Aril  lines.  No  Arl  spectra  was  resolvable  near  the  457.9  nm 
and  458.9  nm  pair.  Thus  two  separate  measurements  are  required  if  this 
technique  Is  to  be  used. 

3.3.3  Spectroscopic  Equipment 

The  spectroscopic  equipment  chosen  for  this  Investigation  was 
an  EG&G  PARC  Optical  Multichannel  Analyzer  (0MA)  III.  The  0MA  system 
consists  of  a  console  (CPU),  detector  controller,  two-dimensional 
vidlcon  detector,  monochromator,  optical  system,  and  translation 
stages.  Figure  18  Is  a  schematic  of  the  system  configuration. 

The  console  is  used  to  set  up  the  scanning  sequence  of  the  detector 
and  acquire  data  in  various  modes.  The  console  also  stores  the  data  via 
a  150  K  floppy  disk  and/or  a  15M  Winchester  hard  disk.  The  detector 
controller  Is  a  16-bit  parallel  peripheral  device  which  accepts  the 
scanning  code  from  the  console  and  derives  the  vidlcon  scanning 
sequence. 

The  detector  Is  a  silicon  Intensified  target  (SIT)  detector,  500  by 
512  pixels.  The  pixels  are  25  by  3  micrometers  with  a  total  active  area 
of  12.2  by  12.5  mm.  The  spectral  range  Is  365  to  800  nm  and  the 


Figure  17  Relative  Aril  Line  Ratio 


Figure  18  Optical  Multichannel  Analyzer  (OMA)  Schematic 
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detector  has  been  coated  with  a  UV  scintillator  which  lowers  the  range 
to  200  m.  The  sensitivity  is  uniform  over  the  active  area  to  ±10 
percent.  The  linearity  is  ±2  percent  but  a  function  of  sensitivity. 
The  sensitivity  Is  10  to  20  photons  per  count  and  the  detector  response 
was  found  to  be  linear  as  long  as  the  counts  per  scan  were  below  1000. 
The  detector  can  also  be  gated  down  to  40  ns  if  temporal  resolution  is 
desired.  The  monochromator  Is  utilized  In  a  spectrograph  configur¬ 
ation.  Thus  the  two-dimensions  of  the  detector  correspond  to  the 
spatial  extent  of  the  entrance  slit  length  and  spectral  spread  of  the 
diffraction  grating.  The  monochromator  Is  a  0.32  In.  Czerny-Turner 
configuration  with  kinematically  mounted  components  such  that  the  system 
can  be  mounted  in  any  orientation.  For  the  majority  of  experiments  in 
this  Investigation,  a  1200  line/mm  holographic  grating  was  used.  This 
results  In  a  resolution  of  .5  to  .6  A  at  the  detector  face. 

3.3.4  Data  Acquisition 

As  was  shown  In  Fig.  18,  the  monochromator  Is  aligned  verti¬ 
cally,  that  Is  with  the  entrance  slit  horizontal.  The  plasma  event  in 
all  cases  Is  vertical  to  minimize  the  effects  of  buoyancy.  The  plasma 
Is  Imaged  onto  the  horizontal  entrance  slit  as  shown  In  Fig.  19.  By 
dividing  the  detector  spatial  dimension  Into  N  Individual  tracks,  the 
entrance  slit  and  hence  the  plasma  radial  dimension  Is  effectively 
divided  Into  N  elements.  The  minimum  number  of  pixels  In  a  given  track 
was  five.  This  number  corresponds  to  the  minimum  number  of  pixels  which 
must  be  grouped  together  to  minimize  the  effect  of  thermal  blooming  at 
the  detector  surface.  This  dimension  Is  the  means  of  calculating  the 
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radial  dimension  of  the  plasma  volume.  For  example,  5  elements  out  of 
500  elements  corresponds  to  the  radial  Increment.  Since  the  500 
dimension  Is  12.2  mm  wide,  the  radial  Increment  at  the  detector  is: 
5/500  x  12.2  mm  -  0.122  mm. 

To  calculate  the  magnification  of  the  monochromator,  one  follows 
the  guidelines  of  Lerner  and  Thevenon  (62).  Assuming  the  optical  system 
f/#  was  matched  to  the  monochromator  f/t  at  the  given  wavelength  [62], 
the  magnification  of  the  monochromator  will  be  on  the  order  of  1  to  1.1 
times.  For  example,  at  the  415.8  nm  Arl  line,  the  monochromator  is  an 
f/4.96  system.  Thus  the  radial  Increment  at  the  entrance  slit 

0  122  imt 

Is  - *  .111  mm.  The  actual  radial  Increment  at  the  plasma  Is 

this  value  divided  by  the  magnification  of  the  optical  system.  From 
Fig.  19  the  optical  magnification  Is  -Sg/Sj.  If  this  value  Is  1/3,  then 
the  actual  increment  Is  0.111/(1/3)  *  0.333  mm.  When  designing  an 
experiment,  one  Is  constrained  to  match  the  monochromator  f/l  at  the 
given  wavelength  and  ensure  that  the  optical  magnification  Is  such  that 
the  plasma  Image  underfills  the  entrance  slit  length.  The  axial 
Increment  Is  somewhat  easier  to  calculate  being  the  slit  width  divided 
by  the  optical  magnification.  Thus  for  the  above  system  with  a  0.1  mm 
slit  width,  the  axial  dimension  Is  0.3  mm,  as  shown  In  Fig.  19. 

Once  the  detector  has  been  divided  Into  N  tracks,  each  track  is 
spectrally  spread  across  the  other  detector  dimension.  The  result  is  a 
multitrack  spectral  data  set  as  shown  in  Fig.  20.  The  x-axis 
corresponds  to  wavelength,  y-axls  Is  radial  position,  an  z-axis 
Intensity  (W/cm*  -  sr).  One  sees  that  the  plasma  emission  Is  "hotter" 
near  the  centerline  ,  as  expected,  and  the  continuum  has  also  Increased 
In  this  region. 
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3.3.4. 1  Instrument  Operation 

The  operation  of  the  OMA  will  now  be  discussed.  Data 
can  be  acquired  in  a  number  of  modes.  Initially,  the  plasma  is  scanned 
in  the  liye  mode  where  the  display  Is  updated  continuously.  The  refresh 
rate  is  determined  by  the  total  number  of  elements  scanned  on  the 
detector.  Typically,  a  large  single  area  (track)  Is  scanned  initially 
to  find  the  plasma.  The  horizontal  translation  stage  is  moved  (Fig.  18) 
to  accomplish  this.  Once  the  plasma  has  been  found  the  detector  is 
divided  Into  larger  numbers  of  tracks  with  the  horizontal  stage  being 
moved  to  again  center  the  plasma.  Plasma  centering  Is  not  particularly 
critical  but  the  plasma  should  be  maintained  close  to  the  center  to 
minimize  the  effect  of  detector  response  which  drops  off  rapidly  in  the 
outer  10  percent  regions. 

Once  the  number  of  tracks  has  been  reduced  to  the  desired  size  or 
minimum  (5  elements)  the  Intensity  of  the  plasma  is  monitored.  The 
detector  can  be  damaged  at  Intensities  in  excess  of  14,000  counts  which 
has  been  shown  to  easily  occur  In  laser  supported  plasmas.  The  local 
Intensity  is  reduced  by  reducing  the  channel  dwell  time  (20  ys  minimum), 
the  reduction  of  the  slit  width,  and  the  addition  of  neutral  density 
filters.  Filters  with  optical  densities  ranging  from  0.5  to  5.0  were 
used  In  this  Investigation  depending  on  the  optical  system  and  laser 
power.  The  plasmas  investigated  In  this  report  are  relatively  strong 
emitters  compared  to  the  events  the  OMA  was  designed  to  monitor.  Thus 
another  possibility  for  Intensity  reduction  is  the  addition  of  an 
aperture  at  the  imaging  lens.  The  aperture  will  change  the  effective 
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f /#  of  the  optical  system  but  will  affect  relative  calculations 
minimally  at  the  detector  face. 

In  this  investigation,  the  data  was  accumulated  over  ten  to  twenty 
consecutive  scans  per  axial  location.  Thus,  the  intensity  had  to  be 
reduced  to  less  than  1000  counts  peak  per  scan  to  insure  linearity  in 
the  summation.  The  accumulated  scans  were  selected  to  minimize  local 
fluctuations  in  the  plasmas.  Fluctuations  were  seen  in  the  live  mode  to 
be  minimal  In  most  plasmas  and  were  at  best  a  random  type  feature.  Once 
the  given  radial  scanned  been  acquired,  the  vertical  translation  stage 
(Fig.  18)  was  incremented  to  the  next  highest  location,  az  typically  1 
to  3  mm. 

As  mentioned  previously,  finding  the  plasma  Is  the  most  difficult 
process.  Once  the  focus  position  of  the  laser  is  known,  a  spectral 
lamp  Is  aligned  In  this  region  with  the  laser  axis.  The 
optical /monochromator  system  Is  then  adjusted  via  the  lamp.  The 
vertical  location  of  the  system  is  aligned  to  a  point  where  the  plasma 
Is  known  to  exist.  Thus,  only  horizontal  adjustments  of  the 
optical /monochromator  system  should  be  necessary  to  find  the  plasma. 
Once  the  plasma  was  found,  absolute  vertical  scaling  is  performed.  In 
the  pure  gas,  a  HeNe  laser  Is  aligned  horizontally  through  the  focus  in 
the  PIFC  and  through  the  optical /monochromator  system.  Once  the  plasma 
is  Initiated,  the  optical  system  Is  translated  downward  (upstream)  until 
the  plasma  emission  ceases,  forming  a  datum  relative  to  the  focus.  In 
the  metallic  plasmas,  a  HeNe  laser  beam  Is  split  as  it  Intersects  the 
target  parallel  to  the  surface.  The  transmitted  beam  is  used  to  align 
the  optical  system  and  the  target  surface  forms  the  datum.  A  brief 
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document  on  the  data  acquisition  procedures  can  be  found  in  the  laser 
lab  (ME).  "Guidelines  to  be  Used  for  Implementing  the  OMA  III  System 
into  an  Experiment,"  is  meant  to  be  used  in  conjunction  with  the  various 
OMA  manuals. 

3.3.5  Data  Reduction 

Once  the  multitrack  data  (Fig.  20)  has  been  stored  the  data 
must  be  reduced  via  Eqns.  (3.3)  or  (3.4).  The  data  that  has  been  stored 
includes  whatever  background  radiation  there  is  in  the  laboratory.  This 
effect  can  be  accounted  for  by  acquiring  data  without  the  plasma 
present.  Typically  one  background  file  Is  required  for  each  set  of 
plasma  scans  as  the  variation  in  the  background  signal  Is  small  along 
the  plasma  axis.  The  background  file  will  also  account  for  any  dark 
signal  and/or  "burned"  elements  which  may  be  present  on  the  detector 
face.  The  file  Is  then  subtracted  from  each  of  the  data  files  via  the 
"Curve  Calculations"  mode  In  the  OMA  III  console.  The  resulting  file  is 
now  pure  plasma  Intensity  data  (W/cm^  -  sr). 

Equation  (3.2)  corresponds  to  the  emission  at  a  pure  frequency.  In 
the  case  of  a  laboratory  plasma,  the  energy  has  been  spread  out  into  a 
profile.  There  has  been  much  debate  as  to  whether  the  peak  or 
Integrated  Intensity  values  are  applicable  In  the  analysis  of 
temperature.  The  peak  Intensity  calculations  generally  result  In  lower 
temperatures  since  It  has  to  be  ratloed  to  a  single  element  continuum 
value.  The  procedure  has  been  used  primarily  In  attempts  to  reduce  the 
amount  of  data  processing  required  In  attempts  to  make  real-time 
temperature  measurement  systems.  The  Integrated  Intensity  of  the  line 


profile  Is  the  figure  of  merit  which  corresponds  to  the  theoretical 
value  of  Eqn.  (3.2).  Since  real-time  data  reduction  is  not  possible 
with  the  current  OMA  software,  the  Integrated  profiles  have  been  used 
throughout  this  Investigation. 

To  integrate  the  line  profile  the  "Data  Utilities"  -  "Channel 
Profile"  routine  on  the  OMA  has  been  used.  The  given  file  is  analyzed 
and  a  starting  and  ending  wavelength  Is  selected.  The  utility  then  sums 
the  intensity  values  from  start  to  end  wavelength,  inclusive.  This 
value  Is  now  proportional  to  the  total  line  and  continuum  energy  in  this 
range  and  the  continuum  contribution  must  be  subtracted  for  use  in  Eqns. 
(3.3)  or  (3.4).  To  do  so,  the  continuum  adjacent  to  the  line  profile  is 
Integrated  over  an  equivalent  number  of  channels  (wavelength).  The 
continuum  file  can  then  be  subtracted  from  the  line/continuum  file, 
yielding  the  approximate  line  term.  Approximation  is  used  here  as 
ideally  one  would  attempt  to  subtract  the  actual  continuum,  which  is  not 
possible  In  the  current  OMA  III  software.  In  older  OMA  systems  when 
integrating  line  profiles,  the  outermost  edges  of  the  band  were 
trapezoldally  fit  to  calculate  the  actual  continuum.  This  feature 
should  be  Included  In  future  OMA  software  releases.  The  net  result  is 
typically  a  slight  underestimation  of  the  continuum  in  the  red  wing  and 
overestimation  In  the  blue  wing. 

Now  that  the  actual  line  and  continuum  file  are  ready  It  should  be 
noted  that  the  values  are  llne-of-sight  integrated  Intensity  (W/cm2  -sr) 
values.  To  accurately  use  the  data,  it  must  be  deconvoluted  into  point 
emissive  powers  (W/cm3  -  sr).  This  Is  accomplished  via  the  Abel 
inversion,  as  shown  In  Fig.  21.  A  numerical  inversion  code  was 
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developed  after  the  work  of  Cremers  and  Birkebak  (631.  The  code  listing 
as  well  as  a  brief  review  of  the  extensive  work  in  this  area  is 
presented  In  Appendix  A. 

A  flow  chart  of  the  numerical  Abel  inversion  code  is  shown  in 
Fig.  22.  The  line  of  sight  Intensity  data  Is  read  into  the  code  via  a 
formatted  data  file.  The  radial  positions  of  a  given  scan  (1  axial 
location)  are  normalized  to  unity  for  ease  in  later  curve  fitting  and 
analysis.  The  code  then  curve  fits  the  data  via  cubic  splines,  the 
number  of  knots  being  determined  by  the  number  of  data  points. 
Typically  7  to  20  data  points  on  a  given  radius  were  resolvable  which 
were  subsequently  cubic  spline  fit  with  three,  four,  or  five  knots.  The 
data  was  relatively  smooth  and  cubic  splines  were  found  to  adequately 
follow  the  trends.  The  cubic  splines  are  analyzed  at  the  data  points 
and  an  average  error  file  is  made.  Now  that  the  cubic  coefficients  are 
available,  the  Inversion  Is  begun. 

The  Abel  Integral  Is  of  the  form  (see  APPENDIX  A): 

I(y)  *  2  /  e(r)  r  (r2  -  y2)-1/2dr  (3.6) 

y 

where  I (y)  Is  the  Integrated  intensity  (W/m2  -  sr),  e(r)  Is  the  emissive 
power  (W/m2  -  sr),  and  y  and  r  as  shown  in  Fig.  21.  The  Integral  can  be 
Inverted  yielding: 

e(r)  *  ~  f  (y2  -  rZ)  /  dy  (3.7) 

If  I(y)  has  been  fit  to  the  form: 
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Figure  22  Abel  Inversion  Code:  Flowchart 
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I(y)  *  A  +  8y  +  Cy2+  Oy3  (3.8) 

then  differentiating  I (y)  and  substituting  Into  Eqn.  (3.7)  yields: 

E(r)  =  -  /  8  ±  2Cy.  t  3Dy^  (3.9) 

*  r  2  2 

✓y  -  r 

Where  the  Integration  begins  at  the  radial  point  r  and  proceeds  to  the 
outermost  radius  In  this  case  normalized  to  unity.  The  equation  Is  a 
simple  Integral  easily  evaluated  analytically  on  the  computer.  The 
derivation  assumes  the  outermost  Intensity  to  be  zero  which  Is  typically 
not  the  case  In  experiment.  To  account  for  this  effect  the  finite  end 
value  of  a  given  scan  Is  extrapolated  linearly  to  zero  at  a  fictitious 
radial  point  at  1  +  ar.  The  slope  of  the  resulting  triangle  is  then 
Integrated  yielding  a  DC  offset  which  is  added  to  Eqn.  (3.9)  at  all 
points. 

The  advantage  of  the  above  evaluation,  outside-in  versus  Inside-out 
Is  that  inversion  errors  are  carried  over  to  the  next  calculation  in 
diminishing  margins.  This  makes  the  inversion  less  error-prone  compared 
to  pervious  matrix  and  numerical  techniques.  Since  the  Inversion  is 
analytical,  data  at  fictitious  points  In  between  real  data  can  be 
evaluated.  The  evaluation  does  not  "blow-up"  as  easily  at  the  center- 
line  and  can  Invert  flattened  centerline  profiles,  typical  of  laboratory 
plasmas.  To  account  for  symmetry  at  the  centerline,  the  first  two 
Inverted  data  points  are  extrapolated  from  a  fourth  order  even  poly¬ 
nomial  curve  fit  of  the  third,  fourth  and  fifth  Inverted  points.  This 
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process  should  occur  in  the  Intensity  data,  but  circumvents  the  possi¬ 
bility  of  centerline  "blowup"  when  performed  on  the  emissive  power 
profile. 

A  typical  error  assessment  of  a  radial  profile  Is  shown  In 
Fig.  23.  The  error  was  analyzed  on  the  basis  of  a  ±10  percent  error  in 
the  measured  Intensities.  Since  the  line  and  continuum  are  measured 
adjacent  to  one  another,  one  would  expect  the  ±  local  error  to  be  in  the 
same  direction.  This  type  of  error  resulted  In  a  total  temperature 
error  of  less  than  1  percent  throughout  any  given  radial  profile.  This 
would  be  expected  if  Indeed  the  Abel  inversion  process  was  working 
property. 

The  other  major  source  of  error  would  be  detector  non-linearity. 
The  line  events,  being  relatively  strong  compared  to  the  continuum, 
would  affect  the  detector  response  over  the  accumulation  of  several 
scans.  This  would  bias  the  total  line  Intensity  low  which  would  result 
in  higher  calculated  plasma  temperatures.  Thus  for  the  error  analysis 
shown  in  Fig.  23  the  line  Intensity  was  varied  ±10  percent  with  the 
continuum  held  constant.  The  result  Is  that  a  maximum  error  of  3.3 
percent  is  seen  at  the  centerline  with  a  decreasing  calculated  error  at 
larger  radial  distances.  This  Is  typical  of  errors  seen  In  other 
studies  (511  except  that  the  outermost  temperatures  have  larger  errors 
due  to  the  local  temperature  by  approaching  the  sub  10,000  K  range. 
Thus  a  conservative  error  of  ±4  percent  was  chosen  and  applied  to  all 
the  temperature  fields  that  will  be  presented  In  Chapter  4. 

Nearly  all  Abel  Inversion  techniques  assume  circular  symmetry  of 
the  Intensity  data.  This  was  found  to  be  the  case  In  the  majority  of 


the  plasmas  In  this  Investigation.  Figure  24  shows  the  Integrated  line 
profile  of  the  415.8  mm  ARI  line  as  taken  from  Fig.  20  using  the  "Curve 
Calculations"  and  "Channel  Profile"  utilities  of  the  OMA.  Once  the  line 
or  line  and  continuum  pair  has  been  inverted  the  ratio  of  the  two 
signals  Is  compared  to  the  ratio  as  calculated  In  Eqns.  (3.3)  or 
(3.4).  The  ratios  have  been  curve  fit  and  the  correct  spline  segment 
for  the  given  ratio  is  applied  to  calculate  temperature.  The  temper¬ 
ature  output  file  Is  formatted  for  subsequent  use  in  global  property 
analysis  and/or  plotting. 

The  absorption  code  flow  chart  shown  In  Fig.  25,  calculates  the  net 
power  absorbed  as  the  focussed  laser  beam  propagates  through  the 
measured  plasma  temperature  field.  The  properties  of  beam  refraction 
and  radiated  power  are  also  calculated.  The  code  listing  Is  presented 
in  APPENDIX  B. 

The  absorption  code  reads  the  temperature  file  which  Includes 
spatial  and  datum  dimensions,  laser  power,  and  laser  beam  geometry. 
Once  the  beam  geometry  Is  determined,  currently  f/7  or  f/2.4  TEM01*,  the 
beam  cross-section  Is  broken  Into  twenty  rays.  The  rays  are  propagated 
at  Increasing  angles  into  the  plasma.  Initially  directed  at  some  radial 
position  offset  from  the  laser  axis  at  the  focus.  In  the  case  of  the 
pure  gas  analysis  the  rays  are  mirrored  past  the  focus  to  account  for 
downstream  absorption.  Each  ray  across  the  profile  is  weighted  by  a 
full  Gaussian  with  the  first  and  twentieth  rays  corresponding  to  the 
1/e-polnts.  The  Gaussian  profile  arises  from  the  nature  of  the  beam  and 
the  fact  that  the  absorption  coefficient  operates  on  intensity,  not 
power.  The  Gaussian  weights  are  summed  prior  to  entering  the  plasma 
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Figure  24  Example  of  an  Integrated  Line  Profile 
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Figure  25  Transport  Code:  Flowchart 


volume,  being  proportional  to  the  Incident  power.  As  the  ray  bundle 
Intersects  the  first  axial  position  of  the  plasma  the  radial  temperature 
profile  Is  cubic  spline  curve  fit.  The  radial  location  of  each  ray  is 
calculated  from  geometrical  optics  and  the  local  ray  temperature  is 
evaluated.  The  local  electron  number  density  is  calculated  from  the 
local  temperature  [27],  The  local  absorption  coefficient  is  calculated 
after  the  work  of  Wheeler  ( 301  and  Stallcop  [31]  as  was  discussed  in 
Chapter  2.  The  local  absorption  is  calculated  using  Eqn.  (2.12) 
where  az  corresponds  to  the  optical  path  length  through  the  experimental 
cell  size.  When  the  rays  have  been  propagated  through  the  entire  plasma 
volume,  the  attenuated  ray  weights  are  summed  which  are  proportional  to 
the  transmitted  laser  power.  The  ratio  of  the  difference  between  Input 
and  output  summations  to  the  Input  summation  is  the  percentage  of 
Incident  power  absorbed. 

A  study  was  performed  to  assess  the  effect  of  refraction  and  datum 
errors.  A  datum  shift  of  ±2  mm  typically  changed  the  net  absorption 
results  by  less  than  ±5  percent.  This  type  of  error  magnitude  has  been 
reported  elsewhere  [23].  When  refraction  of  all  rays  was  taken  Into 
account,  the  difference  In  absorption  with  respect  to  the  straight-line 
geometrical  ray  trace  was  ±0.5  percent,  with  a  maximum  error  of  +1.1 
percent.  However,  beam  refraction  will  be  shown  to  be  a  non-negllgible 
effect  In  metallic  plasmas.  The  beam  refraction  results  In  a  larger 
spot  size  which  greatly  affects  target  transport  phenomena.  Thus  to 
Include  this  effect,  the  beam  refraction  of  the  outermost  ray  is 
accounted  for.  Since  the  experimental  axial  data  Increments  are  larger 
and.  In  general,  the  radial  temperature  gradients  are  large  with  respect 
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to  axial  gradients,  the  radial  index  of  refraction  should  have  the 
greatest  effect.  Again,  this  has  been  seen  in  this  effort  and  elsewhere 
(23). 

Figure  26  shows  the  geometry  under  consideration.  The  twentieth 
ray  Is  of  concern  and  hence  a  twenty-first  temperature  point  Is  linearly 
extrapolated  from  the  nineteenth  and  twentieth  points.  The  index  of 
refraction  at  T^g  and  T£0  (last  two  ray  locations)  Is  calculated  after 
Cheng  and  Casperson  [641.  Snell's  law  Is  used  to  calculate  the  new 
propagation  direction  to  the  next  cell.  Since  the  compliment  of  the 
incident  angle  is  of  concern,  Snell's  law  In  the  radial  direction 
reduces  to: 

nJ9  s1n(90-a)  *  n2Q  s1n(90-8') 

or 

s'  *  cos_1,[-^  cos(s)  1  (3.10) 
n20 

The  index  of  refraction  of  argon  as  calculated  is  less  than  unity  for 
temperatures  In  excess  of  6000  K.  This  simply  Implies  that  the  "group" 
velocity  of  the  photons  Is  greater  than  the  speed  of  light.  The  index 
decreases  with  Increasing  temperature  and  thus  the  ray  is  bent  away  from 
the  centerline. 

To  calculate  the  radiative  losses  from  the  plasma,  a  point-by-point 
volume  analysis  Is  made.  The  power  emitted  by  the  continuum  Is 
calculated  in  Eqn.  (2.15).  The  line  emission  losses  are  determined  from 


Figure  26  Geometry  of  Beam  Refraction  In  Plasma  Volume 
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the  calculated  continuum  emission  after  the  work  of  Kozlov  [65).  The 
details  of  the  radiation  losses,  absorption,  and  thermophysical 
properties  are  treated  in  Chapter  5  on  the  numerical  model.  The  result 
Is  that  a  parameter,  thermal  efficiency.  Is  defined  as: 

P  -  P 
_  absorbed  radiated 

^  ^incident 

Thermal  efficiency  is  a  parameter  which  can  be  directly  compared  to 
thermocouple  measurements  in  the  pure  gas  experiments.  The  parameter  is 
a  value  of  the  net  energy  retained  by  the  gas  that  would  be  available 
for  conversion  to  kinetic  energy.  However,  due  to  errors,  the  value  is 
a  qualitative  assessment  of  the  overall  temperature  field.  Thermal 
efficiency  Is  sensitive  to  the  grid  size  and  a  ±1  percent  change  in  a 
grid  dimension  has  resulted  In  excess  of  ±6  percent  variation  in  thermal 
efficiency.  Thus  throughout  this  effort  thermal  efficiency  is  treated 
as  a  figure  of  merit  which  must  be  greater  than  zero  (power  absorbed  > 
power  radiated)  for  the  data  to  be  considered  viable.  The  other 
calculated  property,  global  absorption,  has  been  seen  to  vary  less  than 
tZ  percent  for  similar  variation  in  axial  dimensions  and  less  than  ±1 
percent  for  radial  errors. 

The  absorption  code  also  calculates  the  cell  by  cell  reradiated 
energy  delivered  to  the  target.  During  the  calculation  of  the  total 
radiation  losses  a  view  factor  between  the  lower  face  of  a  given  cell 
and  the  laser  spot  at  the  target  surface.  The  area  of  the  lower  cell 
face  Is  calculated  and  the  radius  an  equivalent  disk  is  calculated.  The 
view  factor  of  a  disk  of  radius  rc,  height  of  above  the  laser  spot  disk 
of  radius  r$  is  (66): 
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F  =  1/2  (Z  -  /Z2  -  4X2Y2) 
c  s 

where  X  =  rc/z,  Y  =  z/r$,  and  Z  =  1  +(1  +  X2)Y2.  The  total  power 
delivered  to  the  spot  is  summed  over  all  cells.  The  results  will  be 
presented  In  Chapter  4  and  will  confirm  the  first-order  calculations 
detailed  In  Chapter  2  that  if  all  of  the  reradiated  energy  was  coupled 
to  the  target,  the  magnitude  is  less  than  2  percent  of  the  transmitted 
laser  power.  Further,  the  total  reradiated  power  is  on  the  order  of  a 
few  watts  and  should  make  little  difference  to  the  target  areas  outside 
the  laser  spot. 


4.  EXPERIMENTAL  RESULTS 


The  previous  chapter  outlined  the  various  data  acquisition  and 
reduction  techniques  utilized  in  this  investigation.  All  techniques  are 
applied  to  each  radial  profile  at  a  given  axial  plasma  location.  All 
axial  locations  are  then  stacked  to  form  a  two-dimensional  temperature 
mapping.  This  chapter  will  present  representative  two-dimensional 
temperature  fields  for  qualitative  comparison.  Quantitative  comparison 
will  be  made  on  the  basis  of  global  absorption  and  radiative  loss  terms. 

The  first  section  will  develop  the  pure  argon  plasma  results.  The 
quantitative  analysis  is  made  as  a  function  of  power  and  flowrate.  The 
results  are  then  compared  with  Independent  measurements  as  a  consistency 
check  of  the  spectroscopic  data.  The  second  section  discusses  the 
applicability  of  relative  line  and  line  to  continuum  diagnostics  in 
metal-gas  plasma  spectra.  Quantitative  comparison  of  plasma  behavior  is 
made  as  a  function  of  power,  translation  speed,  and  surface  treatment. 
The  pure  aluminum  results  have  been  used  in  conjunction  with  a  transport 
model  to  demonstrate  the  effects  of  plasma  formation  on  melt/solid 
interface  growth. 

4.1  Pure  Argon  Plasmas 

The  pure  argon  temperatures  have  been  determined  using  relative 
line  to  continuum  techniques  as  discussed  in  Chapter  3.  The  data  shown 
In  the  following  section  corresponds  to  the  415.8  nm  Arl  line  and 
adjacent  blue  wing  continuum.  Figure  27  is  a  typical  quasi-three- 
dimensional  plot  of  a  5.5  kW,  f/2.4,  argon  plasma  temperature  field  at 
15  cm/sec.  Figure  27a  Is  a  view  from  the  upstream,  27b  Is  a  side-on 
view,  and  27 c  a  downstream  view  of  the  temperature  field. 


One  notes  the  rough  edges  depicted  in  Fig.  27.  The  rippling  at  the 
edges  are  due  to  the  wing  errors  associated  with  the  Abel  inversion. 
The  effect  was  discussed  in  Chapter  3.  The  edges,  however,  are  at  sub- 
10,000  K  temperatures  and  thus  have  little  effect  on  laser  attenu¬ 
ation.  The  edge  effects  are  also  minimized  when  the  geometrical 
propagation  path  of  the  laser  beam  in  Fig.  27  is  taken  into  consider¬ 
ation.  The  main  effect  of  wing  error  is  in  the  radiation  calcu¬ 
lations.  Although  relatively  small  temperatures  are  encountered  in 
these  regions,  the  local  volumes  are  large  compared  to  the  centerline 
regions.  This  makes  the  wing  contributions  to  radiative  losses  non- 
negllgible  and  results  In  the  relatively  large  error  bars  for  thermal 
efficiency  data,  that  will  be  seen  In  subsequent  figures. 

Welle,  Keefer,  and  Peters  1231,  using  absolute  Intensity  calcu¬ 
lations,  showed  off-axis  peaks  in  certain  regions  of  the  pure  argon 
plasmas.  The  off-axis  peaks  are  directly  attributable  to  the  Abel 
Inversion  and  therefore  may  not  be  real  events.  The  Abel  inversion  can 
result  In  a  reduced  emissive  power  near  the  centerline.  This  reduction 
can  be  Interpreted  as  a  lower  or  higher  electron  temperature  depending 
on  what  temperature  the  peak  emission  corresponds  to.  The  off-axis 
peaks  did  not  occur  at  higher  flow  rates  In  Ref.  (23).  Nonetheless,  the 
off-axis  and  centerline  peaks  result  in  negligible  differences  in  the 
calculated  global  properties  when  comparing  the  results  herein  to 
Ref.  (231. 

The  plasma  begins  at  -10  mm  (±0.5  mm)  upstream  of  the  laser 
focus.  The  peak  temperature  at  approximately  -3.5  mm  Is  19,000  K.  The 
core  is  approximately  30  mm  long  and  10  mm  wide  at  the  10,000  K 
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Isotherm,  the  lower  limit  of  the  spectroscopic  analysis.  The  three- 
dimensional  mapping  is  obtained  by  folding  the  radial  data  about  the 
centerline.  In  the  case  depicted  in  Fig.  27,  the  global  absorption  was 
calculated  to  be  82.9  percent  (4560  W).  The  reradiated  power  was  3653  W 
leaving  907  W  retained  by  the  gas  for  a  thermal  efficiency  of  16.5 
percent.  The  values  will  be  shown  to  be  within  experimental  error  if 
values  calculated  from  independent  measurements  and  model  predictions. 
The  results  also  agree  with  the  values  reported  elsewhere  [231 • 

4.1.1  Power  Dependence 

Figure  28  depicts  the  spatial  dependence  of  plasma  extent  as 
a  function  of  power  at  f/2.4  and  15  cm/sec  argon  flow  velocity.  Figure 
28a  Is  for  the  case  of  2.4  kW,  28b  is  4.1  kW,  and  28c  is  5.5  kW  (as  was 
shown  In  Fig.  27).  The  peak  temperature  In  Figs.  28a-c  shows  a  general 
rise  with  power  going  from  17,500  K  to  19,000  K  in  Figs.  28a  and  28c. 
The  peak  temperature  in  Fig.  28b  Is  only  17,000  K  which  is  due  to  the 
error  In  the  Abel  Inversion  near  the  centerline.  The  more  important 
trend  is  that  the  17,000  K  core  region  has  grown  significantly  In  volume 
from  low  power  to  high  power.  The  core  region  and  plasma  front  have 
also  moved  further  into  the  upstream  flow  as  power  Is  increased. 
Finally,  the  plasma  length  has  grown  significantly  between  2.4  kW  and 
marginally  between  4.1  kW  and  5.5  kW.  The  net  result  will  be  seen  to  be 
an  overall  Increase  In  absorption  and  radiative  loss  with  increasing 
laser  power. 

The  calculated  values  of  global  absorption  are  60.0  percent  (1440 
W,  72.9  percent  (2990  W),  and  82.9  percent  (4560  W)  for  2.4,  4.1,  and 
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5.5  kW  Incident  laser  power,  respectively.  Figure  29  compares  the 
results  to  the  values  calculated  by  a  numerical  model  [22]  and  trend 
agreement  within  experimental  error  is  indicated.  The  error  bars 
represent  the  fluctuation  about  the  data  point  given  a  ±4  percent 
temperature  error.  The  temperature  error  was  applied  uniformly  across 
the  temperature  grid  and  thus  represents  a  worst  case  scenario.  The 
highly  non-linear  properties  of  absorption  coefficient  and  radiation 
loss  are  the  main  contributions  for  error  bands  for  absorption  and 
thermal  efficiency,  respectively.  A  more  realistic  approach  to  error 
analysis  would  be  to  segment  the  error  based  on  spatial  location  as  was 
demonstrated  In  Chapter  3.  Finally,  the  data  points  of  Fig.  29 
represent  a  summation  over  400  to  1000  temperature  cells  and  this  type 
of  error  behavior  Is  not  unexpected.  The  Important  point  Is  that  the 
error  bands  Indicate  trends  similar  to  the  data. 


The  thermal  efficiency,  defined  as  the  difference  between  power 
absorbed  and  power  lost  via  radiation  divided  by  Incident  power,  has 
also  been  calculated.  In  Fig.  29,  the  results  again  are  shown  to  agree 
with  the  modeling  effort  within  experimental  effort. 


One  notes  the  relatively  large  error  bars  on  the  efficiency  data. 


Since  the  power  radiated  Is  a  function  of  spatial  extent,  l.e.  radial 


and  axial  increment,  two  more  variables  and  associated  errors  are 


Involved.  In  the  absorption  calculation,  radius  was  shown  to  have 


little  effect  while  axial  dimensions  are  dominated  by  the  absorption 


coefficient  in  a  decaying  exponential  term.  The  two  variables  are 


directly  used  to  calculate  the  unit  volume  In  the  radiation  terms  and 


thus  ±10  percent  of  radial  and  axial  dimensions  significantly  alter  the 
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results.  Thus  It  should  be  mentioned  that  the  thermal  efficiency  and 
radiative  losses  are  used  in  a  qualitative  sense  throughout  this  inves¬ 
tigation.  The  thermal  efficiency  must  be  greater  than  zero  to  avoid 
second  law  violations  that  Imply  radiated  power  to  be  greater  than  ab¬ 
sorbed  power.  The  Important  use  of  the  efficiency  data  Is  that  the 
trends  of  the  numerical  model  have  been  verified.  The  thermal  effi¬ 
ciency  was  calculated  to  be  28.2  percent  (680  W),  28  percent,  (1150  W), 
and  16.5  percent  (910  W)  for  the  cases  of  2.4,  4.1,  and  5.5  kW, 
respectively. 

The  agreement  between  the  spectroscopic  and  numerical  model  are  not 
unexpected.  As  will  be  seen  In  Section  4.1.3,  the  spatial  agreement  of 
plasma  volume  is  reasonable.  Although  minor  differences  exist  between 
the  two  methods  of  calculation,  the  thermophysical  properties  used  were 
identical.  Thus  based  upon  the  spatial  agreement  with  respect  to 
Isothermal  mappings,  global  absorption  and  thermal  efficiency  should  be 
in  close  agreement  as  was  shown  in  Fig.  29.  The  true  consistency  check 
on  the  spectroscopic  data  was  made  by  comparison  with  thermocouple  and 
calorimeter  data. 

The  calorimeter,  as  shown  In  Fig.  5,  Is  mounted  above  the  plasma  in 
the  downstream  region  of  the  PIFC.  The  laser  geometry  Is  carefully 
aligned  to  Insure  that  all  transmitted  laser  energy  Impinges  on  the 
walls  of  the  calorimeter.  When  the  plasma  Is  stabilized,  the  trans¬ 
mitted  laser  energy  is  measured.  A  correction  for  stagnant  gas  heating 
of  the  calorimeter  has  been  made.  Thus,  knowing  the  Incident  and  trans¬ 
mitted  laser  energy,  the  absorbed  energy  can  be  calculated. 
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The  thermal  efficiency  Is  measured  via  locating  the  thermocouple 
carriage  in  the  downstream  region  of  the  plasma.  By  assuming  some 
velocity  profile  In  this  region  whose  average  velocity  equals  the  uni¬ 
form  flow  field  upstream  of  the  plasma,  a  mixing  cup  analysis  can  be 
performed.  The  downstream  enthalpy  flux  Is  calculated  and  compared  with 
the  Inlet  enthalpy  flux  to  determine  the  net  power  retained  by  the 
gas.  Unfortunately,  the  plasma  radiation  effects  on  thermocouple 
measurements  have  yet  to  be  quantified  and  the  thermal  efficiency 
measurements  are  a  qualitative  Indicator  of  trends. 

Figure  30  depicts  the  comparison  between  the  various  experimental 
methods  as  well  as  the  numerical  model.  Reasonable  agreement  between 
calorimetric  and  spectroscopic  data  is  indicated  at  all  powers.  Thermal 
efficiency  agrees  In  trend  and  to  a  lesser  extent  In  magnitude.  The 
result  Is  that  the  spectroscopic  results  have  been  demonstrated  to  be 
reasonable  and  consistent  such  that  the  spectroscopic  data  may  be  used 
independently  of  the  other  diagnostics  as  an  assessment  of  plasma 
behavior.  The  result  was  critical  In  the  evaluation  of  metallic  target 
plasma  as  will  be  detailed  In  Section  4.2. 

4.1.2  Flow  Rate  Dependence 

Figure  31  shows  the  effect  of  Increasing  flow  rate  on  plasma 
spatial  extent  at  a  fixed  laser  power.  The  results  are  for  4.1  kw, 
f/2.4  at  15,  60,  and  120  cm/sec  flowrates  in  31a,  b,  and  c, 
respectively.  It  can  be  seen  that  the  plasma  is  elongated  into  the 
downstream  region  significantly  as  flow  rate  Is  Increased.  Figure  31b 
and  c  have  been  expanded  two  and  three  times,  respectively,  in  the 
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Figure  30  Comparison  of  Spectroscopically  Determined 

Parameters  with  Calorimetry  and  Thermometry  Data 
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radial  dimension  for  clarity.  Increasing  the  flow  rate  from  15  to  60 
cm/sec  nearly  doubles  the  length  of  the  plasma.  The  net  absorption  in 
this  case  Is  reduced  since  the  absorption  length  gains  are  offset  by  an 
overall  reduction  in  temperature.  The  peak  temperature  is  further 
reduced  to  15,400  K  in  the  high  flow  case.  Fig.  31c. 

The  plasma  front  has  been  pushed  closer  to  the  focal  spot  of  least 
confusion  with  Increasing  flow  rate.  The  plasma  is  also  narrower 
approaching  a  cigar  shape  versus  the  tear  drop  shape  at  low  flow 
rates.  One  also  notes  temperature  peak  downstream  of  the  focus.  At 
first  this  was  thought  to  be  an  artifact  of  the  Abel  Inversion.  How¬ 
ever,  nearly  all  of  the  high  flow  cases  (>  60  cm/sec),  at  any  power, 
showed  similar  trends.  A  review  of  the  measured  Intensity  data 
occasionally  showed  local  llne-of-slght  Intensity  Increases  In  this 
region.  This  leads  one  to  speculate  on  the  possibility  of  downstream 
eddy  mixing  regions.  This  Is  a  distinct  possibility  based  upon  the 
recent  results  of  Ref.  [241,  a  full  solution  of  the  Navler  Stokes 
equations.  The  results  of  Ref.  [241  Indicate  that  relatively  little  or 
no  mass  flow  Is  occurring  In  and  about  the  central  regions  of  the  plasma 
near  the  laser  focus.  This  would  support  streamline  formation  with 
components  normal  to  the  plasma  axis  immediately  downstream  of  the 
focus.  The  results  of  Ref.  [241  do  not  indicate  this  flow  feature,  as 
eddy  generation  Is  not  treated  in  the  effort.  The  effect  cannot  be 
verified  until  proposed  flow  measurement  diagnostics  have  been  applied. 

Figure  32  shows  the  flow  velocity  dependence  of  absorption  and 
thermal  efficiency  for  the  spectroscopic  data  at  5.5  kW.  For  comparison 
the  model  [221  results  at  5  kW  are  shown.  One  sees  that  the  total 
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absorption  decreases  with  Increasing  flow  velocity.  The  overall  plasma 
temperature  is  lowered  and  dominates  any  effects  of  increased  plasma 
length.  The  dip  in  the  model  prediction  could  not  be  verified  with  che 
limited  number  of  data.  The  thermal  efficiency  data  indicates  that 
there  is  an  optimal  flow  velocity  for  a  given  power  and  beam  geometry, 
as  calculated  by  the  model.  The  peak  arises  again  from  geometrical 
considerations.  At  low  flow  rates  the  plasma  length  and  temperature 
define  absorption  behavior.  As  the  flow  rate  is  increased  the  lower 
plasma  temperatures  reduce  radiative  losses  and  increasing  thermal 
efficiency  is  the  result.  With  further  increases  in  flow  rate  the 
radiative  losses  are  further  reduced.  The  effect  is  that  absorbed 
energy  roughly  approaches  some  minimum  value,  just  prior  to  blow  out, 
where  radiative  losses  are  negligible.  Unfortunately,  this  regime 
exists  at  very  low  plasma  temperatures  where  the  net  absorption  is 
minimal . 

Figure  33  shows  the  results  for  the  cases  of  2.4  and  4.1  kW  versus 
flow  velocity.  In  general,  the  trends  are  similar  to  Fig.  32. 
Decreasing  absorption  with  flow  velocity  and  indication  of  a  peak  in 
thermal  efficiency.  Another  qualitative  trend  that  the  efficiency  peak 
occurs  at  a  lower  velocity  for  a  lower  power  is  indicated.  This  trend 
was  also  predicted  by  the  modeling  effort. 

4.1.3  Comparison  with  a  Numerical  Model  [22) 

Figure  34  is  a  comparison  of  the  spatial  extent  of  an  argon 
plasma  as  calculated  from  the  spectroscopic  data  (34a)  and  the  model 
(34b).  One  notes  the  spatial  agreement  of  the  isotherms  as  mentioned 


Figure  33  Global  Laser  Power  Absorption  as  a  Function  of 
Flow  Velocity:  Pure  Argon,  f/2.4 
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Figure  34  Comparison  of  Spatial  Extent  of  Spectroscopic 
Data  and  a  Numerical  Model  (22) 
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previously.  Since  the  absorption  and  radiation  calculations  are  treated 
similarly,  the  previously  shown  agreement  between  the  data  and  model  was 
expected.  The  agreement  between  the  spectroscopic  and  calorimeter  and 
thermocouple  data  demonstrated  consistency  in  trends  and  magnitudes. 
This  allows  one  to  apply  the  diagnostics  to  more  complicated  plasmas. 

From  the  data  and  model  predictions,  new  plasma  regimes  of  Interest 
have  been  developed.  First,  the  elongated,  cooler  plasma,  to  a  certain 
extent,  retained  more  useful  energy  in  the  downstream  gas.  This  trend 
was  indicated  by  all  measurements.  Thus,  the  next  step  is  to 
geometrically  elongate  the  plasma  and  operate  at  higher  flowrates  to 
obtain  peak  efficiencies.  This  will  be  accomplished  by  the  use  of 
higher  f/#  optics  to  confine  the  plasma  In  a  longer,  narrower  cone. 
Also  proposed  is  a  means  of  reducing  radiative  losses  via  multiple 
plasma  formation.  The  effect  Is  largely  geometrical  as  the  radiative 
view  factors  to  the  wall  have  been  reduced  by  the  plasma-to-plasma  view 
factors.  Also,  each  of  two  plasmas  operates  at  half  power  resulting  in 
a  lower  percentage  of  radiative  losses  than  in  the  case  of  one  plasma  at 
twice  the  power.  The  other  area  of  interest  is  higher  pressures  where 
radiation  losses  dominate  due  to  Increased  numbers  of  emitters  and 
absorbers.  Once  the  two  proposed  regimes,  f /#  and  multi-plasma,  as  well 
as  Increased  pressure  have  been  quantified  an  accurate  assessment  of 
scaling  laws  from  the  model  can  be  made. 

4.2  Aluminum  Targets  In  Argon  Atmospheres 

Throughout  this  portion  of  the  Investigation,  an  f/7  beam  geometry 
was  used.  This  was  due  to  the  Inability  to  initiate  the  plasma  on 
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aluminum  targets  at  higher  f/#s  with  the  10  kW  AVCO  C02  or  with  the 
Gaussian  beam  of  the  1.5  kW  Photon  Sources  C02  lasers.  The  plasma 
formation  has  been  studied  with  respect  to  surface  treatment  and 
translation  speed. 

4.2.1  Comparison  of  Data  Analysis  Techniques 

Ideally,  one  would  prefer  to  utilize  the  same  data  reduction 
method  for  the  aluminum  data  as  was  used  In  the  pure  argon  studies. 
However,  Is  was  unknown  a  priori,  the  extent  of  aluminum  effects  In  the 
argon  plasma.  Minimally,  the  aluminum  was  expected  to  bias  the 
continuum  emission  rendering  relative  line  to  continuum  techniques 
Invalid.  Thus,  a  pair  of  argon  lines  which  could  be  used  In  a  relative 
technique  were  sought. 

The  spectra,  as  was  seen  In  Fig.  15,  has  resolvable  Aril  lines 
which  meet  the  criteria  for  relative  line  diagnostics.  This  effect  is 
attributed  to  the  aluminum  presence  as  the  Aril  spectra  was  not 
observable  In  the  pure  gas  cases.  Several  authors  had  used  the  457.9 
and  458.9  nm  Aril  lines  successfully.  Two  problems  exist.  First,  the 
Aril  lines  exist  only  In  the  hotter  central  regions  of  the  plasma.  The 
Arl  spectra  must  be  calibrated  to  the  Aril  temperatures  In  order  to 
completely  map  the  plasma.  The  process  Is  relatively  complicated  and 
requires  approximately  twice  the  number  of  data  points  compared  to  the 
relative  line  to  continuum  method.  Secondly,  the  457.9  and  458.9  nm 
Aril  lines  are  separated  from  any  Arl  lines  which  have  a  reasonable 
slgnal-to-nolse  ratio  above  the  continuum.  Thus  two  Independent 
mappings  of  the  plasma  emission  are  required  which  further  Increases 
error  In  the  data. 


128 


The  457.9  and  458.9  nm  Aril  temperature  data  was  compared  to  the 
415.8  nm  Arl  line  to  continuum  data.  If  it  could  be  shown  that  the 
profiles  agreed  within  experimental  error,  or  that  global  properties 
agreed  within  experimental  error,  the  easier  line  to  continuum  method 
could  be  used  for  the  metal  vapor/argon  plasmas. 

The  analysis  was  performed  on  aluminum/argon  plasmas  at  3,  5,  and  7 
kW  Incident  laser  power.  Figure  35  is  a  typical  radial  temperature 
profile.  Figure  35a  is  the  415.8  nm  to  continuum  data,  35b  the 
457.9/458.9  nm  and  35c  the  696.5  nm  Arl  to  continuum  data.  As  can  be 
seen,  fair  agreement  between  all  three  techniques  Is  Indicated  with  most 
of  the  temperature  field  agreeing  within  experimental  error.  Table  1  is 
a  listing  of  the  calculated  values  of  absorption  and  thermal  efficiency 
for  the  case  of  5  kW  plasma.  It  can  be  seen  that  the  global  results  are 
well  within  experimental  error,  further  Indicating  the  accuracy  of  the 
relative  line  to  continuum  technique.  Finally,  the  aluminum  line 
spectra  were  not  resolvable  except  in  the  region  within  2  mm  of  the 
target  surface.  This  further  supports  the  proposal  to  analyze  the 
plasma  as  pure  argon.  The  strongest  aluminum  lines  resolved  were  the 
394.4  and  396.6  nm  resonant  lines  which  extended  a  few  mm  into  the 
plasma,  but  typically  less  than  30  percent  of  the  overall  plasma  length. 


Table  I 


Technique 

pabs<*> 

"th(*> 

Rspot(' 

457.9/458.9  nm  Aril 

16.46 

12.46 

0.506 

696.5  nm  Arl 

14.22 

11.83 

0.504 

415.8  Arl 

16.42 

13.49 

0.508 

Figure  35  Comparison  of  Typical  Radial  Temperature 

Profiles  using  Various  Reduction  Techniques 
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Thus  with  a  reasonable  amount  of  confidence,  the  remainder  of  the 
aluminum/argon  data  will  be  reduced  using  the  relative  line  to  continuum 
technique.  In  light  of  the  large  errors  shown  in  previous  investi¬ 
gations  [ 53,54]  the  above  error  has  at  least  been  demonstrated  to  be 
within  the  experimental  accuracy  of  the  various  techniques.  Also,  when 
the  software  for  the  OMA  Is  fully  operational,  it  may  prove  useful  to 
analyze  the  larger  data  sets  required  for  relative  line  techniques.  The 
Inability  to  transfer  data  from  the  OMA  to  another  computer  system 
except  manually  Is  the  main  reason.  The  Hemenway*  BASIC  language 
supplied  as  an  option  on  the  OMA  does  not  permit  data  reduction  and 
analysis  on  the  OMA.  The  BASIC  system  may  not  permit  the  access  of 
certain  Intrinsic  functions  as  well  as  external  functions  necessary  for 
data  reduction.  The  BASIC  system  is  also  extremely  slow  compared  to  the 
CYBER  or  HP  systems  available. 

In  the  analysis  of  the  aluminum/argon  data  an  additional  parameter 
was  calculated.  The  refraction  of  the  laser  beam  Is  calculated  as  the 
ray  bundle  propagates  through  the  plasma.  In  this  Investigation  and 
elsewhere  [23]  beam  refraction  was  shown  to  have  little  effect  on  global 
absorption  (t0.5  percent).  Thus,  for  simplicity,  the  refractive  index 
of  the  outermost  ray  Is  calculated  as  detailed  In  Chapter  3.  The 
refractive  Index  variation  In  the  plasma  Is  such  that  the  beam  Is 
refracted  away  from  the  laser  axis.  The  effect  will  be  shown  to  be 
considerable  In  Section  4.2.6  when  the  experimental  results  are  utilized 
In  a  numerical  model. 
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4.2.2  Power  Dependence 

Figure  36  is  a  quasi  three-dimensional  mapping  of  the 
temperature  field  of  a  5  kW,  f/7  plasma  above  a  99.9  percent  pure 
aluminum  target  In  an  argon  atmosphere.  The  plasma  extends 
approximately  18  mm  above  the  target  surface  (datum  *  0.00)  which  agrees 
with  visual  observation  (0.5  to  0.75  in.)  the  plasma  Is  nearly 
cylindrical  as  would  be  expected  for  an  f/7  beam  and  an  18  mm  axial 
dimension.  Several  temperature  peaks  along  the  axis  are  noticeable. 
These  are  possibly  due  to  Induced  flow  fields  within  the  plasma  as  the 
argon  atmosphere  was  a  slow  flow  with  Inlet  and  outlet  ports  for  from 
the  plasma  (6  to  8  In.).  Optical  characteristics  of  the  windows  (l.e. 
reflections)  could  not  be  seen  when  using  a  spectral  lamp  source.  The 
likely  possibility  Is  error  In  the  Abel  Inversion  as  all  the  peaks  in 
the  5  to  10  mm  region  are  roughly  within  experimental  error  of  one 
another. 

Figure  37  depicts  the  temperature  fields  for  4,  5,  and  7  kW,  f/7 
aluminum/argon  plasmas.  One  sees  the  general  trend  of  an  Increase  in 
plasma  length  and  volume  as  well  as  overall  temperature.  In  particular 
In  the  5  and  7  kW  cases,  plasma  spreading  well  beyond  the  beam  geometry 
Is  indicated  near  the  target  surface.  This  may  be  indicative  of 
Increased  conduction  effects  in  these  regions  between  the  plasma  and  the 
target. 

Figure  38  shows  the  calculated  values  of  absorption  as  a  function 
of  total  Incident  power.  The  absorbed  power  rises  relatively  rapidly  at 
low  powers  and  again  appears  to  approach  some  asymptotic  value  at  powers 
above  7  kW.  The  thermal  efficiency  Is  a  meaningless  value  In  this 
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portion  of  the  study  but  is  used  as  a  benchmark  to  indicate  the  magni¬ 
tude  of  the  radiative  losses.  The  thermal  efficiency  is  12.5  percent  at 
4  kW,  rises  to  15  percent  at  5  kW  and  flattens  out  above  5  kW.  The 
radiative  transport  was  estimated  via  the  methodology  outlined  In 
Chapter  2  and  has  been  shown  to  account  for  less  than  1.5  percent  of  the 
total  energy  (Including  transmitted  laser  energy)  deliverable  to  the 
target. 

Detailed  calculations  were  made  on  a  point-by-point  basis  in  random 
plasma  fields  and  further  supported  the  above  calculations.  The 
radiative  power  of  each  cell  was  calculated  and  the  solid  angle  view 
factor  between  the  lower  face  of  the  cell  and  the  laser  spot  was 
calculated.  The  net  power  delivered  from  each  cell  to  the  spot  Is 
summed  and  compared  to  the  delivered  laser  power.  The  results  indicated 
that  the  radiation  effects  account  for  less  than  0.5  percent  of  the 
total  delivered  power.  Since  the  view  factors  between  cell  side  faces 
were  not  taken  Into  account,  the  1.5  percent  figure  of  the  first  order 
calculation  can  be  regarded  as  a  conservative  estimate.  At  first  this 
result  was  disconcerting  In  the  light  of  the  almost  universal 
acknowledgment  of  an  enhanced  coupling  regime  when  a  plasma  forms. 

The  effect  of  beam  refraction  was  later  incorporated  into  the 
transport  analysis,  as  detailed  In  Chapter  3.  The  results  are  depicted 
In  Fig.  39.  Figure  39  Indicates  that  spot  size  increases  with 
Increasing  laser  power.  The  effect  is  a  result  of  the  larger  radial 
temperature  gradients  and  Increasing  plasma  length  with  increasing  laser 
power.  The  effect  Is  Intuitively  obvious.  The  index  of  refraction  is 
inversely  proportional  to  temperature.  Higher  temperatures  near  the 
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centerline  thus  tend  to  bend  the  laser  rays  away  from  the  centerline, 
increasing  spot  size.  From  the  temperature  profiles  shown  in  Fig.  37  it 
can  be  seen  that  the  effect  is  continuous  throughout  the  plasma 
volume.  The  error  bars  in  Fig.  39  result  from  the  uncertainty  of  the 
datum  (±0.5  mm).  However,  the  trend  is  indicated  and  will  be  shown  to 
be  significant  in  the  determination  of  the  target  melt-solid  interface 
in  Section  4.2.5. 

4.2.3  Translation  Speed  Dependence 

Figure  40  depicts  the  spatial  dimensions  of  the  plasma  as  a 
function  of  target  translation  speed.  The  effect  of  increasing 
translation  speed  is  a  lowering  of  plasma  temperature  and  volume. 
Figure  41  shows  the  trend  of  the  calculated  power  absorption  as  a 
function  of  translation  speed.  The  60  in/min  speed  was  nearing  the 
limit  where  the  plasma  begins  to  deflect  out  of  the  laser  geometry, 
affecting  the  symmetry  assumption. 

Experimental  observation  indicates  that  the  trend  is  not  continuous 

as  Fig.  41  would  suggest.  There  appeared  to  be  a  threshold  somewhere 

between  the  20  and  60  In/min.  values  where  the  plasma  volume  suddenly 

became  smaller  at  the  given  laser  power.  The  threshold  was  apparent  at 

all  laser  powers  between  4  and  7  kW.  The  nominal  intensity  at  4  kW,  1.0 

fi  7 

mm  beam  diameter,  was  1.6  x  10  W/cm  .  The  corresponding  interaction 
time  was  0.118  and  0.0393  seconds  for  the  20  and  60  in/min.  speeds, 
respectively.  The  data  taken  in  this  region  was  very  noisy,  further 
supporting  a  threshold  condition.  The  noise  did  not  permit  reliable 
assessment  of  plasma  parameters. 
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Two-dimensional  Temperature  Field  of  a  4  kW 
f/7  Aluminum- Argon  Plasma: 

(a)  10,  (b)  20,  and  (c)  60  in/min 


The  absorption  lowering  Is  largely  due  to  the  reduced  temperatures 
between  lower  speeds.  Absorption  lowering  Is  dominated  by  the  smaller 
propagation  path  length  at  the  higher  speeds,  above  the  threshold.  The 
threshold  Is  subjective  however,  since  the  plasma  radiation  Is  Intense 
enough  to  distort  visual  observation. 

4.2.4  Graphite  Coated  Aluminum  Targets 

Graphite  and  other  absorbing  coatings  are  commonly  used  to 
enhance  laser  coupling  to  targets.  A  brief  set  of  experiments  was 
performed  to  determine  whether  a  graphite  coating  affected  plasma 
formation.  Figure  42  Is  a  comparison  of  the  4  kW  plasma  temperature 
fields  for  coated  and  non-coated  aluminum  targets. 

One  sees  that  the  graphite  has  little  or  no  effect  on  the  steady 
state  plasma  volume  and  temperature.  No  carbon  lines  could  be  resolved 
In  any  of  the  spectra  at  any  laser  power.  The  only  effect  that  may 
occur  Is  a  lowering  of  the  plasma  threshold,  but  this  could  not  be 
verified  In  this  set  of  experiments.  The  calculated  laser  power 
absorption  for  the  graphite  coated  targets  was  14.3  percent  which 
compares  well  with  the  pure  aluminum  results  of  Fig.  38. 

The  experiment  was  also  performed  at  5,  6,  and  7  kW,  f/7,  but  the 
Intensity  profiles  were  too  noisy  for  Abel  Inversion  and  were  not 
reduced.  However,  the  visual  observation  of  the  plasma  volume 
correlated  with  pure  aluminum  measurements.  Future  experiments  may 
attempt  to  monitor  the  graphite  aided  plasma  formation  effect.  If  any, 
by  gating  the  OMA  detector  off  the  laser  or  optically  triggering  from 
the  plasma  emission.  At  this  point,  an  Intuitive  guess  would  be  that 


V  V 


142 


graphite,  although  useful  at  sub-plasma  intensities,  serves  little  or  no 
purpose  above  threshold. 

4.2.6  Effect  of  Results  on  Target  Transport  Model 

The  effect  of  the  calculated  beam  attenuation  and  refraction 
were  analyzed  In  a  target  transport  model.  A  three-dimensional 
conduction  model  (67]  was  used  to  evaluate  the  target  surface  temper¬ 
ature  distribution  and  melt-solid  interface.  Nominal  parameters  of  5.0 
kW,  f/7,  1.0  mm  spot  size  and  the  calculated  5  kW  plasma  parameters  of  4 
kW  transmitted,  1.04  mm  spot  size. 

Figure  43  depicts  the  results.  The  higher  peak  temperature  profile 
corresponds  to  the  5  kW  (no  plasma)  case,  which  Is  as  expected. 
However,  both  temperature  profiles  peak  well  above  the  melting  point  of 
aluminum  and  significant  melting  has  occurred  In  the  target.  The 
interesting  result  Is  depicted  In  the  melt/solid  Interface  predic¬ 
tions.  The  plasma  case  actually  results  In  a  larger  melt  volume,  a 
counter-intuitive  result. 

However,  the  larger  melt  volume  Is  Indicative  of  the  enhanced 
coupling  regime,  noted  by  several  references  in  the  literature  survey. 
The  theory  for  the  effect  In  CW  radiation  is  as  follows.  The  plasma  is 
of  the  weakly  absorbing  (LSC)  type.  Therefore  sufficient  power  is 
transmitted  to  the  target  for  vaporization  and  the  plasma  to  be 
supported.  The  main  contribution  of  the  plasma  Is  the  refractive 
bending  of  the  focused  laser  beam  out  of  focus.  Since  significant 
melting  has  occurred  at  the  target  under  plasma  conditions,  the  dominant 
parameter  In  the  process  Is  Interaction  time.  Interaction  time  is 
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defined  as  the  laser  spot  size  divided  by  the  target  translation 
speed.  Thus  a  larger  spot  size  results  in  a  larger  interaction  time  or 
a  longer  beam  duration  on  any  given  point  on  the  target. 

The  effect  should  occur  at  nearly  all  laser  powers  as  the  trend  Is 
for  beam  refraction  to  increase  with  laser  power.  Further  modeling 
efforts  have  been  hindered  to  date  by  the  Inability  of  the  target  model 
to  converge  at  higher  laser  powers.  The  result  Is  that  if  these 
difficulties  can  be  overcome  and  a  plasma  model  can  be  developed, 
complete  a  prior  determination  of  melt  volume  (HAZ)  can  be  achieved. 
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5.  NUMERICAL  MODEL 

Several  efforts  have  been  made  in  modeling  the  steady-state  plasma 
formation  for  CW  laser  radiation  in  pure  gases.  A  complete  review  of 
these  efforts  can  be  found  In  Ref.  20.  Early  models  were  typically  one¬ 
dimensional,  using  cylindrical  (non-focused)  beam  geometries.  Simpli¬ 
fying  assumptions  regarding  the  non-linear  thermophysical  properties  of 
the  high  temperature  gas  were  made  to  form  analytical  solutions.  More 
recently,  converging  beam  geometries  and  real  gas  properties  have  been 
incorporated  In  numerical  solutions.  One  of  the  first  successful  two- 
dimensional,  real  gas,  converging  beam  geometry  models  was  reported  by 
Glumb  and  Krler  (22).  Recently  a  complete  solution  of  the  incom¬ 
pressible  Navler-Stokes  equations  for  a  plasma  system  has  been  reported 
by  Jeng  and  Keefer  [24] .  This  investigation  reports  the  development  of 
a  model  following  the  procedure  of  Glumb  and  Krier  [22].  The  results  of 
Jeng  and  Keefer  [24]  will  be  discussed  a  possible  future  enhancements. 

5.1  Model  Formulation 

The  physics  of  plasma  formation  above  metallic  targets  were 
discussed  in  Chapter  2.  The  main  difference  between  this  type  of  plasma 
and  a  pure  gas  plasma  Is  the  existence  of  an  upstream  boundary  condition 
(target  surface)  which  affects  the  plasma  formation  via  the  vaporization 
characteristics  of  the  target  material.  Figure  44  Is  the  cylindrical 
domain  used  In  the  solution  procedure.  For  the  first  model,  zero  radial 
flow  Is  assumed.  Further,  since  the  axial  Mach  numbers  have  been 
measured  to  be  on  the  order  of  0.05  to  0.1,  axial  convection  has  also 
been  neglected.  The  results  of  Jeng  and  Keefer  [24]  In  flowing  gas 
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Figure  44 


Cylindrical  Domain  of  Metal-Gas  Plasma  Model 
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systems  Indicate  nearly  zero  axial  flow  velocities  In  the  plasma  core, 
supporting  the  zero  axial  flow  assumption. 

The  boundary  conditions  are  as  follows: 

a)  centerline  symmetry  (adiabatic), 

b)  free  boundaries  at  upper-  and  outer-most  cylindrical  grid 
locations,  and 

c)  heavy  particle  temperature  distribution  known  at  the  target 
surface. 

Boundary  condition  a)  Is  self-explanatory.  Boundary  condition  b) 
is  such  that  a  fictitious  final  grid  temperature  Is  extrapolated  from 
the  two  outermost  iterated  grid  temperatures.  The  conduction  losses 
across  the  two  faces  will  then  be  calculated  and  used  in  the  convergence 
test  of  the  energy  equation.  Boundary  condition  c)  results  from  a 
target  transport  model  (671.  The  temperature  distribution  is  close  to 
Gaussian  with  the  peak  temperatures  being  a  few  degrees  above  the 
melting  point  of  the  target  material. 

With  the  above  assumptions  and  boundary  conditions,  the  energy 
equation  takes  the  form: 

l7fKll  +  ?lfClt#>|-‘’“Cpll  +  »I-Pr*0(5a) 

Where  K  Is  the  temperature  dependent  conductivity,  puCp  is  the  convec¬ 
tion  term,  al  is  the  absorbed  laser  power,  and  Pr  Is  the  radiative  loss 


term. 


5.2  Solution  Procedure 

Equation  5.1  can  be  reduced  to  a  difference  form  using  the  grid 
shown  In  Fig.  45. 
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corresponds  to  the  average  conductivity  between  the  cells  In 
question.  Equation  (5.2)  and  the  modified  forms  at  the  various 
boundaries  can  be  solved  for  as  explicit  forms  of  the  surrounding 
cell  temperatures.  The  non-linearity  of  the  gas  properties  prevents 
direct  matrix  reduction  of  the  system  and  a  Gauss-Seldel  Iterative 
procedure  Is  used  as  In  Ref.  22. 

The  laser  beam  Is  focused  Into  the  plasma  grid,  currently  set  at 
f/7.  For  simplicity,  the  beam  Is  modeled  as  a  Gaussian  annular  inten¬ 
sity  profile.  The  Iteration  begins  at  the  upper  central  grid  and 
rastors  radially  and  then  axially,  etc.  The  beam  Is  geometrically 
propagated  Into  the  grid  and  each  cell  along  the  radius  Is  checked  to 
determine  the  percentage  of  the  north  boundary  exposed  to  the  laser 
beam.  The  local  absorption  coefficient  and  absorbed  power  is  then 
calculated  for  use  In  Eqn.  (5.2).  The  average  conductivity  between  N, 
S,  E,  and  W  cells  and  the  center  cell  Is  calculated  as  well  as  the 
radiative  loss  term  at  Tc(0LD).  The  results  are  used  In  Eqn.  (5.2)  to 
determine  th  new  grid  temperature,  T^(HEW).  The  procedure  Is  amenable 
to  the  Incorporation  of  an  over-  or  under-relaxed  Gauss-Seldel  Iteration 
procedure  where  the  new  cell  temperature  Is  determined  by: 
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tnew  s  tnew  +  w<tnew  '  tolo^  ^5*3^ 

W  Is  the  weighting  factor,  0-1  for  the  under-,  and  greater  than  unity 
for  the  over-relaxed  case. 

Once  a  given  Iteration  has  been  performed  the  radiative  losses  and 
laser  absorption  are  summed  over  the  entire  grid  and  an  energy  balance 
of  the  form: 


Pabs  "  Prad  "  PC0N0,  N.S.E  "  PC0NV,N  *  0 
is  checked  for  convergence. 

5.3  Argon  Properties 

Argon  has  been  used  In  this  model  since  it  was  also  used  in  the 
pure  gas  model  and  argon  Is  frequently  used  as  an  assist  gas  In  laser 
material  processing.  The  electron  number  density  was  calculated  from 
LTE  considerations  and  the  equilibrium  after  Dresvln  [27],  The  values 
Include  contributions  from  the  second  Ionization  stages  at  temperatures 
above  16,000  K.  Figure  46  shows  the  variation  of  electron  number 
density  with  electron  temperature.  The  decline  In  population  above 
17,000  K  corresponds  to  the  slight  decrease  in  ion  production  until  the 
second  ionization  stage  contributes  significantly  at  about  20,000  K. 

The  Inverse  bremsstrahlung  absorption  coefficient  is  calculated 
following  the  work  of  Wheeler  [ 30 ]  and  Stallcop  [31)  and  was  shown  as 
Eqn.  (2.3).  The  net  result  Is  that  Wheeler's  corrections  account  for 
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photonization  and  non-hydrogen ic  behavior.  Stal loop's  effective  Gaunt 
factor  accounts  for  electron-neutral  IB  absorption  at  the  lower 
temperature  range.  Figure  47  shows  the  absorption  coefficient  behavior 
as  a  function  of  electron  temperature. 

Thermal  conductivity  as  calculated  using  gas  kinetic  theory  does 
not  agree  with  experimental  measurements.  This  can  be  seen  in 
Fig.  48.  The  discrepancy  arises  from  the  inability  of  experiment  to 
discriminate  between  thermal  conductivity  and  contributions  from 
radiative  absorption.  Thus,  the  experimental  data  can  be  used  in  the 
modeling  effort  and  radiative  transport  within  the  plasma  will  be 
automatically  accounted  for. 

Radiation  losses  are  more  complex  as  both  continuum  and  line 
radiation  terms  are  present.  As  was  shown  in  Chapter  2,  the  continuum 
emission  is  accurately  modeled,  on  an  integrated  basis,  by  the  Kramers- 
Unsold  theory  [58].  However,  line  radiation  must  be  summed  over  all 
energy  levels,  a  difficult  tasks  at  best.  Several  experimental  efforts 
concluded  that  the  line  radiation  was  the  same  magnitude  as  continuum 
losses  and  thus  the  continuum  was  doubled.  The  result  is  valid  up  to 
approximately  13,000  K.  Kozlov  [65]  derived  an  expression  for  the  line 
radiation  as  a  multiplication  factor  of  the  continuum  losses  at  a  given 
electron  temperature.  The  calculation  was  made  on  the  basis  of  a 
smeared  line  spectrum  where  the  line  energy  was  modeled  as  perturbations 
on  the  continuum  centered  at  the  wavelength  in  question.  In  this 
manner,  closely  spaced  discrete  line  spectra  could  be  lumped  into  an 
analytical  expression.  The  total  emission  Is  shown  in  Fig.  49  with  the 
continuum  contribution  (Eqn.  2.15)  shown  for  reference.  More 
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Information  on  the  argon  properties  can  be  found  In  the  thesis  of  Dr. 
Ron  Glumb  [49]. 

5.4  Results:  Comparison  with  Experiment 

Figure  50  depicts  the  flow  diagram  for  the  plasma  model.  The 
spatial  dimensions  of  the  plasma  grid  and  laser  beam  are  read  from  a 
data  file.  The  Initial  temperature  grid  Is  set  to  a  constant  value  of 
10,000  K  or  greater,  otherwise  a  previous  solution  grid  Is  read. 

The  boundary  conditions  are  as  mentioned  before.  The  upper  and 
outer  boundaries  are  currently  set  to  300  K  simply  to  minimize 
convergence  time.  The  current  code  takes  between  1000  and  2000 
Iterations  to  reach  convergence  of  the  energy  equation. 

There  was  a  concern  over  the  target  surface  boundary  condition.  It 
was  desired  to  find  a  relationship  between  the  experimentally  observed 
electron  temperatures  In  excess  of  10,000  K  and  the  predicted  heavy 
particle  temperature  far  below  this  value.  Initially  a  non- 
thermochemlcal  equilibrium  scheme  such  as  proposed  by  Lukens  and 
Incropera  [41]  was  desired.  The  authors  cited  Imbalance  between  ohmic 
heating  and  elastic  electron  energy  transfer  as  a  source  of  thermal  non¬ 
equilibrium  between  electrons  and  heavy  particles.  The  effort  was  for  a 
cylindrical  arc  plasma  In  which  thermal  conductivity  at  the  anode  walls 
could  be  neglected  with  respect  to  ohmic  heating.  Clearly  this  Is  not 
the  case  for  a  laser  supported  plasma  and  the  unwieldy  task  of  perform¬ 
ing  the  thermal  conductivity  balance  was  undertaken.  The  results  of 
predicting  a  temperature  difference  based  on  the  local  E-fleld  magnitude 
were  suspect,  especially  at  E-flelds  above  a  few  hundred  V/cm,  which  Is 
the  case  for  laser  radiation. 
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Figure  50  Metal -Gas  Plasma  Model  Flow  Diagram 
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After  it  was  demonstrated  by  the  modified  target  model  that  the 
heavy  particle  surface  temperature  distribution  varied  significantly 
with  target  material.  It  was  decided  to  use  this  as  the  lower  boundary 
condition.  The  Knudsen  jump  approximation  (40)  is  applied  to  the 
predicted  surface  temperature  distribution,  lowering  the  temperature 
between  5  and  35  percent,  depending  on  keyhole  pressure  conditions. 
This  effect  follows  the  laboratory  observation  that  the  plasma  formed 
above  a  steel  target  Is  larger  than  a  plasma  above  and  aluminum  target, 
and  thus  the  predicted  surface  temperature  distribution  has  been  used  as 
a  first  approximation. 

The  conduction  heat  transfer  to  the  adjacent  nodes,  power 
absorption,  and  radiative  losses  are  calculated  for  each  cell.  The 
algorithms  are  nearly  Identical  to  those  used  in  the  data  reduction 
code.  The  exception  Is  power  absorption.  In  the  data  analysis  code, 
the  laser  beam  was  broken  Into  20  rays,  each  weighted  to  a  Gaussian 
annular  Intensity  profile.  The  Integration  (piece-wise  summation)  of 
the  inlet  and  exiting  ray  bundle  defined  power  absorption.  The 
intensity  Is  necessary  as  the  local  absorption  coefficient  operates  on 
intensity,  not  power,  as  defined  in  Chapter  2.  In  the  model  the  local 
intensity  is  assumed  to  be  an  annular  Gaussian  distribution  and  local 
power  Is  determined  by  Integrating  the  "absorbed"  intensity  over  the 
beam  area. 

The  model  was  initially  run  with  zero  convection.  The  result  is 
that  the  model  overpredicts  absorption  when  compared  to  the  experimental 
results.  Figure  51  Is  a  comparison  of  the  zero  flow  condition  and  the 
experimental  results  that  were  shown  in  Fig.  38.  The  model  predicts 


Plasma  Model  Predictions:  Zero  Convection 
as  a  Function  of  Laser  Power 
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approximately  3  times  the  absorbed  laser  power  with  respect  to  the 
experimental  results.  This  would  be  physically  impossible  as  the  power 
transmitted  to  the  target  would  be  Insufficient  to  support  vaporization 
and  hence  the  plasma.  The  experimental  targets  indicated  that  the  laser 
had  coupled  to  the  target  throughout  the  given  trial,  further  supporting 
the  above  conclusion. 

After  the  mechanics  of  the  model  code  were  checked,  two  possible 
solutions  existed.  First  there  Is  a  large  variation,  in  some  cases  in 
excess  of  ±50  percent.  In  thermophysical  properties  from  author  to 
author  [19,22,241.  One  approach  would  be  to  alter  the  thermophysical 
properties  within  the  model.  The  exercise  was  performed  on  the 
absorption  coefficient  and  the  model  predictions  agreed  with  the  data  to 
within  the  experimental  error.  However,  the  model  has  now  been  made 
empirical  which  Is  not  the  goal.  More  Importantly,  the  full  property 
values  had  been  used  in  the  data  analysis  algorithms  and  thus  a  paradox 
existed. 

The  second  alternative  was  to  Include  convection  In  the  model  and 
was  Incorporated  via  Eqn.  (5.2).  Recently,  the  target  transport  model 
[67]  has  predicted  Mach  numbers  on  the  order  of  0.3  outside  the  Knudsen 
layer.  From  the  predicted  heavy  particle  temperature,  assuming 
atmospheric  pressure  within  the  plasma,  the  gas  density  and  speed  of 
sound  were  calculated.  The  prediction  of  M  *  0.3  appears  high  with 
respect  to  measured  velocities  Implying  M  *  0.01  to  0.005.  Thus,  the 
Mach  number  was  left  as  a  free  parameter  Input  to  the  model. 

The  Mach  number  was  varied  and  by  trial -and-error  a  Mach  number  of 
0.0001  was  found  to  result  In  a  model  solution  which  agreed  with  the 
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experimental  results.  The  manner  in  which  the  Mach  number  is 
incorporated  into  the  plasma,  model,  the  corresponding  gas  velocity  at 
M  =  0.0001  is  approximately  100  mm/sec.  Figure  52  Is  a  comparison  of 
the  M  =  0.0001  predicted  temperature  field  and  the  experimental  results 
as  taken  from  Fig.  37b. 

One  sees  that  fair  spatial  agreement  has  been  obtained.  The  model 
still  overpredicts  absorption  at  the  given  Mach  number.  At  larger  Mach 
numbers  the  model  grid  relaxes  to  300  K  throughout,  which  Is  interpreted 
as  plasma  "blowout".  Figure  53  shows  the  model  predictions  as  a 
function  of  laser  power.  The  experimental  results  are  shown  for 
comparison. 

The  relatively  low  flow  velocity  Is  not  unexpected  In  light  of  the 
experimental  measurements  (67)  and  the  predictions  of  a  full  two- 
dimensional  model  (241.  Both  note  low  axial  flow  velocities  within  the 
plasma  core  at  Mach  numbers  on  the  order  of  0.01  or  less.  The  measured 
flow  velocities  (671  were  between  80  and  120  mm/sec.  Thus  the  model 
solution  at  100  mm/sec  appears  to  be  a  reasonable  first  Iteration.  The 
main  conclusion  is  that  the  plasma  core  flowfleld  must  be  quantified 
prior  to  significant  revisions  of  the  current  plasma  model.  The  surface 
boundary  condition  should  also  be  Investigated  further. 


Plasma  Model  Predictions:  Convection  with 
u  =  100  mm/sec  as  a  Function  of  laser  Power 


6.  CONCLUSIONS 


The  pure  gas  plasma  studies  are  one  of  the  first  treatments  of  CW 
laser-sustained  plasma  at  powers  up  to  10  kW.  At  the  onset  of  the 
project  many  concerns  were  expressed  concerning  plasma  generation. 
First,  Initiation  was  to  be  a  problem  due  to  the  high  laser  intensities 
required  for  breakdown  of  a  room  temperature  gas.  This  was  circumvented 
by  the  use  of  thin  zinc  foils  and  more  recently  a  reuseable  tungsten 
rod.  Initiation  is  now  routine  and  our  techniques  have  been  used  by 
others  [231. 

Another  area  of  concern  was  plasma  maintenance  and  stability,  it 
was  thought  that  the  plasma  may  extinguish  due  to  laser  fluctuations  or 
flowfleld  variations  near  the  focus.  Neither  laser  or  flowfield 
fluctuations  were  seen  to  alter  the  plasma  behavior.  The  effect  may 
become  critical  at  higher  pressures  and  flowrates,  but  the  numerical 
modeling  [22,24)  and  other  experiments  [23]  have  not  Indicated  any 
alteration  in  plasma  stability. 

Radiative  loading  of  the  chamber  was  a  concern.  Experiments  have 
been  limited  to  sub-ten  minute  Intervals  with  cooling  periods  in  between 
to  prevent  chamber  damage.  As  thermal  efficiencies  are  Improved  the 
problem  may  be  circumvented.  Other  engineering  techniques  could  be 
utilized  such  as  reflective  Inner  chamber  walls  and  regenerative  cooling 
of  the  chamber  by  the  propellant  system.  The  argon  systems  have  shown 
promise  for  laser  propulsion.  No  critical  area  has  been  demonstrated 
which  would  be  fatal  to  the  concept  of  laser  propulsion.  Thermal 
efficiencies  on  the  order  of  40  percent  and  absorption  In  excess  of  75 
percent  have  been  demonstrated. 
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The  spectroscopy  work  is  one  of  the  first  two-dimensional  studies 
of  laser  supported  plasmas  over  a  wide  range  of  laser  power  and  flow- 
rate.  The  work  is  the  first  study  in  which  the  spectroscopic  data  was 
checked  for  continuity  against  independent  experimental  diagnostics. 
Agreement  between  the  spectroscopic  data  and  the  experimental  data  [21 
and  the  numerical  model  [22]  have  allowed  the  development  of  new  regimes 
of  plasma  propulsion. 

To  reduce  radiative  losses,  higher  f/#'s  have  been  theoretically 
shown  feasible.  The  concern  for  plasma  stability  at  f/#'s  above  f/10  is 
questionable  in  light  of  the  results  of  the  other  concerns  addressed 
previously.  The  current  system  has  been  limited  to  f/3  or  less  and  the 
current  plans  are  to  Install  f/8  to  approximately  f/16  optics. 
Splitting  of  the  laser  beam  Into  multiple  focli  is  another  proposal  to 
be  tested.  It  is  theorized  that  two  smaller,  more  efficient  plasmas 
will  operate  more  efficiently  as  a  whole  rather  than  a  single  plasma. 
The  radiative  loading  of  the  chamber  walls  would  also  be  reduced  by  the 
energy  absorbed  between  the  plasmas. 

The  laser-supported  plasma  propulsion  concept  appears  feasible  and 
the  types  of  engineering  refinements  previously  mentioned  should  be 
investigated.  In  lieu  of  a  full  scale  laser  system,  a  more  complete 
mapping  of  plasma  characteristics  is  prerequisite  to  defining  scaling 
laws  for  full-scale  systems. 

With  respect  to  the  hardware  and  software  required  to  acquire  and 
reduce  the  spectroscopic  data,  several  Items  are  in  need  of 
improvement.  First,  the  bottleneck  in  the  system  Is  data  transfer 
between  the  acquisition  and  reduction  stages.  As  future  software 
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enhancements  from  EG&G  are  released,  the  in-house  software  should  be 
generated  to  dump  OMA  data  to  the  CYBER  or  HP  computing  systems.  The 
OMA  has  an  RS-232  serial  port  for  this  function,  but  to  date,  the  port 
has  not  been  addressable. 

The  Abel  inversion  used  in  the  data  reduction  is  functional,  not 
sophisticated.  Many  efforts  have  been  devoted  to  the  subject  (see 
Appendix  A)  and  more  sophisticated  algorithms  are  available.  The 
question  remains  as  to  whether  the  university's  computational  software 
will  support  the  required  utilities.  Another  possibility  would  be 
remote  access  of  the  algorithms.  The  hardware  for  the  data  acquisition 
was  designed  to  be  flexible  and  should  support  most  applications  with 
only  minor  adjustments. 

A  first  attempt  at  a  numerical  model  has  also  been  undertaken.  The 
model  currently  predicts  power-dependent  plasma  volumes  that  are 
consistent  with  the  experimental  results.  The  model,  as  run  for 
aluminum,  predicts  no  plasma  formation  for  laser  powers  below  4  kW  at 
f/7,  which  Is  also  consistent  with  experiment.  The  model  currently 
overpredicts  global  absorption  when  compared  with  the  experimental 
results.  The  absorption  trends,  however,  follow  experiment. 

The  two-dimensional  spectroscopic  studies  of  CW  laser-supported 
plasmas  above  metallic  targets  are  the  first  to  be  found  in  the  liter¬ 
ature.  The  conclusions  are  based  on  the  consistency  of  the  spectro¬ 
scopic  data  in  pure  gas  systems  with  respect  to  independent  measurement 
techniques.  The  plasmas  were  found  to  have  little  or  no  metallic 
spectra  in  the  emission  profiles  and  were  treated  as  pure  argon  plasmas, 
which  is  in  agreement  with  other  studies  (10,51,52).  The  plasmas  were 
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experimentally  determined  to  absorb  15  to  35  percent  of  the  incident 
laser  radiation  at  laser  powers  between  4  and  7  kW.  Temporal  measure¬ 
ments  of  the  plasma  event  Indicated  that  the  plasma  is  a  CW  event  and 
does  not  initiate,  propagate,  extinguish,  and  re-initiate  as  proposed  by 
another  Investigation  [171.  OMA  system  was  used  at  50  ys  scanning  times 
with  scan-to-scan  times  of  100  us  for  temporal  resolution.  Future 
experiments  might  try  to  verify  the  results  with  transmission  measure¬ 
ments  using  a  fast  photodetector.  The  plasma  was  also  shown  to  radi- 
atlvely  behave  essentially  as  a  pure  argon  plasma.  This  conclusion  was 
drawn  from  the  inability  to  resolve  aluminum  spectra  except  for  resonant 
structures  within  the  first  few  millimeters  of  the  target  surface.  This 
also  simplified  data  analysis  and  Is  consistent  with  other  studies. 

Future  efforts  might  reconsider  the  near  surface  behavior  of  the  plasma 
to  determine  if  the  metallic  species  significantly  affect  plasma 
absorption.  Intuitively,  the  effect  should  be  at  most  a  few  percent. 

More  Importantly,  the  plasma  has  been  shown,  in  theory,  to 
significantly  refract  the  laser  beam,  altering  the  energy  distribution 
at  the  target.  When  recently  coupled  with  a  three-dimensional  target 
conduction  model  [671,  the  result  was  that  the  HAZ  was  increased  upon 
plasma  formation.  This  supports  the  expected  enhanced  coupling  regime 
noted  by  other  authors  as  reported  In  Chapter  1.  A  metallographlc  study 
of  the  laser/plasma  Irradiated  targets  should  be  performed  in  a  future 
effort.  The  study  would  Indirectly  support  the  refraction  modeling  and 
could  be  used  to  refine  both  the  target  and  plasma  models. 

The  above  result  Is  significant  as  It  alone  appears  to  explain  the  ; 

enhanced  coupling  effect.  Many  theories  have  been  proposed  Including 
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plasma-target  conduction  and  radiative  transport.  In  the  data  analysis 
it  was  seen  that  conduction  was  small  compared  to  the  transmitted  laser 
power.  Radiative  transport,  as  calculated  using  view  factors  on  a 
point-by-point  plasma  volume  analysis  was  also  shown  to  be  negligible 
even  when  assuming  blackbody  behavior  by  the  target. 
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APPENDIX  A:  ABEL  INVERSION 

In  Chapter  3,  the  Abel  inversion  technique  that  is  used  in  this 
Investigation  was  discussed.  This  appendix  presents  a  detailed 
treatment  of  previous  efforts  and  the  details  of  the  Inversion  code 
"NABEL\ 

Early  on  in  spectroscopic  work  it  was  realized  that  line  of  sight 
measured  values  of  an  Inhomogeneous  plasma  volume  did  not  necessarily 
conform  to  local  properties.  Assuming  circular  symmetry,  the  line-of- 
sight  intensity  can  be  written: 

I(y)  -  2  {  ■fclrdr-  (M/cmz  -  sr)  (A.l) 

»  ,r2-y2 

Typically,  the  measured  quantity  is  the  llne-of-slght  Intensity, 
I(W/cm2  -  sr).  Equation  A.l  Is  of  the  form  of  Abel's  Integral  which  can 
be  Inverted  to  yield: 

e(r)  *  I  JR  L4UM41  (W/cm3  _  sr)  (A.2) 

n  r  1  ? 

✓y  -  r 

From  this  basic  set  of  equations  a  large  effort  has  been  made  by  many 
Investigators  or  analytical  and  numerical  solutions. 

One  of  the  first  treatments  was  made  by  Bockasten  [68],  The 
Intensity  distribution  was  substituted  by  a  third  degree  polynomial 
which  Is  constrained  to  zero  slope  at  the  centerline  from  consideration 
of  symmetry.  A  series  of  matrix  coefficients  were  then  calculated, 
centerline  out,  such  that  the  local  emisslvlty 
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£(r)  *  r01  aiK  \  (A*3^ 

where  a1K  are  the  matrix  elements,  NK  the  intensity  data,  and  rQ  the 
outermost  radius  for  normalization.  The  disadvantages  of  the  techniques 
are  that  by  starting  at  the  centerline  the  error  In  a1K  is  a  minimum  at 
the  outer  radii.  This  error  Is  then  propagated  and  summed  inward  when 
the  matrix  reduction  occurs.  However  the  technique  works  reasonably 
well  for  smooth,  noisefree  data,  especially  if  the  radial  profile  is 
near  Gaussian. 

Barr  [69)  derived  a  similar  treatment  using  second  order 
polynomials.  The  main  contribution  was  a  substitution  for  the  dl/dy 
term  In  Eqn.  (A. 2)  such  that  the  Integration  (smoothing)  was  performed 
prior  to  differentiation  (amplification  of  noise).  The  result  was 
reduced  error  compared  to  the  previous  efforts.  However,  both 
techniques  accumulate  errors  from  outside  In.  This  results  in  the 
prediction  of  negative  emissive  powers  near  the  centerline  if  the 
profile  Is  flattened.  This  is  the  case  In  many  laboratory  plasmas  and 
In  a  majority  of  the  plasmas  In  this  Investigation.  Edels,  Hearne,  and 
Young  [70]  performed  an  analysis  similar  to  Barr  [69],  The  authors 
broke  the  profile  Into  equal  sub-ranges  and  the  Integral  at  each  point 
was  replaced  by  the  sum  (Taylor  series)  of  a  series  of  integrals  over 
the  outer  sub-ranges.  The  result  was  that  a  correction  factor  could  be 
assessed,  as  well  as  errors,  such  that  corrections  to  the  Inversion 
could  be  made  numerically.  However,  the  method  still  has  similar 
problems  near  the  centerline. 
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More  recent  efforts  have  included  that  of  Cremers  and  Blrkebak  [63] 
as  used  in  this  investigation,  the  use  of  orthogonal  polynomials  [71], 
and  derivative  free  inversion  processes  [72].  The  last  three  inversion 
techniques  have  ben  demonstrated  to  have  the  ability  to  invert  flattened 
profiles  accurately.  The  Cremers  and  Blrkebak  was  the  simplest  to 
program  and  demonstrated  equivalent  errors  as  the  other  two 
techniques.  The  method  by  Oeutsch  and  Beniaminy  [72]  involves  the 
integration  by  parts  of  Eqn.  (A. 2)  and  an  expansion  of  the  integran 
assuming  a  well  behaved  function.  This  treatment  appears  to  be  the  most 
rigorous  to  date.  The  code  Is  available  on  the  CPC  Library  at  Queen's 
University  of  Bellfast,  N.  Ireland.  The  technique's  results  were  not 
satisfactory  in  this  investigation  and  upon  personal  communication  with 
Dr.  Oeutsch  It  was  concluded  that  the  IMSL  cubic  spline  utilities  were 
not  accurate  enough  for  the  method  to  succeed.  It  may  prove  desirable 
in  the  future  to  investigate  the  routine  further. 

One  final  technique  worth  noting  is  that  of  Yasutomo,  et  al. 
[73).  The  Inversion  process  has  been  extended  in  this  effort  to 
accommodate  asymmetrical  profiles.  The  distribution  is  divided  into 
even  and  odd  (symmetrical  and  asymmetrical)  components  and  a  relative 
weighting  function  defined.  The  weighting  function  is  a  measure  of  the 
asymmetry  of  the  data.  The  even  part  was  Inverted  using  Barr's  [69] 
method.  The  Inverted  signal  was  then  calculated  via  the  product  of  the 
inverted  even  profile  and  the  asymmetry  function.  The  technique  was 
developed  for  the  highly  asymmetrical  profiles  seen  in  torroldal 
plasmatrons.  The  technique  may  prove  useful  If  future  Investigations 
find  plasmas  to  be  asymmetrical,  e.g.,  horizontal  laser  axis  where 
buoyancy  may  affect  the  plasma  volume. 


The  numerical  Able  inversion  code  "NABEL"  Is  currently  supported  on 
the  laser  lab  UNIX®  system.  The  code  Is  accessible  in  the 
laser/todd/abeldata  directory.  Data  file  formatting  Is  listed  in  the 
initial  comments  within  the  code.  The  code  listing  follows  this  text 
and  the  following  Is  a  brief  walk  through  the  process. 

Several  data  statements  at  the  beginning  of  the  code  are  used  In 
the  calculation  of  temperature.  The  data  statements  contain  the 
coefficients  of  the  415.8  nm  llne-to-contlnuum  ratio  as  a  function  of 
electron  temperature.  Once  the  data  file  has  been  read  the  normalized 
position  vector  corresponding  to  the  data  points  Is  calculated.  A  dummy 
vector  Is  set  at  all  data  points  and  one  point  In  between  data  to  be 
used  In  the  Inverted  data  calculations.  The  number  of  data  points  Is 
checked  and  the  number  of  spline  segments  Is  determined  based  on 
previous  error  analysis.  The  radial  profile  (line  and  then 
corresponding  continuum)  Is  the  cubic  spline  fit  using  the  IMSL  utility 
ICSVKU.  Initial  knot  location  estimates  are  made  and  the  utility 
iterates  the  knot  locations  to  minimize  total  curve  error.  Once  the 
coefficients  have  been  calculated,  the  error  Is  calculated  at  each  data 
point  and  the  average  data  error  per  axial  profile  is  printed  In  the 
output  file  "error". 

The  method  of  Cremers  and  Blrkebak  |63]  assumes  that  the  outer  most 
Intensity  values  are  zero.  Since  the  experimental  data  rarely  obeys  the 
rule,  a  correction  Is  made.  The  finite  outermost  intensity  value  Is 
linearly  extrapolated  to  zero  at  the  N  +  1  radial  point.  The  linear 
function  has  been  Integrated  In  Eqn.  (A. 2)  and  the  result  Is  an  offset 
value  added  to  the  Abel  kernel  throughout  the  remaining  evaluation  of 
the  given  profile. 
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The  inversion  begins  at  the  outermost  point  and  segment  of  the 
radial  profile.  Equation  (3.9)  Is  used  to  calculate  the  inverted 
value.  At  the  junction  of  two  spline  segments  the  kernel  is  summed  over 
the  previous  spline  kernel.  The  Integration  continues  between  the 
current  radial  point  and  the  outermost  knot  of  the  given  spline.  The 
process  continues  to  the  centerline.  The  third,  fourth,  and  fifth  data 
points  from  the  centerline  are  fit  by  a  fourth  order  even  polynomial  and 
the  first  and  second  data  points  are  extrapolated.  This  process  insures 
symmetry  at  the  centerline,  one  of  the  Initial  assumptions. 

The  ratio  of  the  Inverted  line  and  continuum  emissive  powers  are 
compared  to  the  Initial  spline  knots  (data  statements).  Once  the  proper 
spline  segment  has  been  chosen  the  electron  temperature  Is  calculated 
from  the  spline  coefficients.  The  ratio  Is  also  checked  to  Insure  it  is 
within  the  temperature  limits.  If  not,  the  data  point  is  flagged  high 
(20,000  K)  or  low  (6,000  K)  which  is  easily  detected  in  the  output 
file.  The  output  file  is  then  formatted  in  a  rectangular  grid.  If  a 
given  radial  profile  has  fewer  points  than  the  largest  radial  profile, 
the  wings  are  strapped  artificially  to  9,000  K.  The  9,000  K  value  was 
chosen  to  provide  a  skirt  for  the  quasi-three-dlmensional  plotting  code 
as  the  spectroscopic  analysis  cannot  be  used  below  10,000  K.  The 
9,000  K  value  Is  not  used  In  the  absorption  code  and  thus  has  no  effect 
on  radiation  calculations. 

The  functions  at  the  end  of  the  code  are  the  analytical  forms  of 
the  first,  second,  and  third  degree  Integrals  of  Eqn  (3.9).  The  outpile 
file  "tempout"  Is  formatted  for  direct  use  in  either  the  plotting  or 
transport  analysis  routines. 
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APPENOIX  B:  TRANSPORT  CODE  FOR  EXPERIMENTAL  DATA 

The  code  "ABSORBS"  was  developed  to  calculate  the  laser  beam 
attenuation,  refraction,  and  radiative  transport  from  a  predefined 
temperature  field  (see  Appendix  A).  The  code  sets  up  an  Incident  laser 
beam  intensity  profile  and  divides  the  beam  Into  20  rays.  Each  ray  is 
propagated  at  a  measured  angle  until  the  first  axial  plasma  temperature 
profile  is  encountered.  The  temperature  profile  is  cubic-spline  fit  (8 
knots)  and  the  temperature  at  each  ray  location  Is  evaluated.  The 
temperature  dependent  absorption  coefficient  Is  evaluated  and  each  ray 
is  attenuated  by 

*0UT  *  *IN  exP(-aA^)  (B.2) 

where  aZ  is  the  optical  path  length  through  the  given  cell. 

The  temperature  dependent  number  density  (27)  and  subsequently  the 
index  of  refraction  are  calculated  for  each  ray  position  where 

n(r)  *  (1  -  up2/w2)  (B.3) 

where  wp  is  the  plasma  (electron)  frequency  and  w  is  the  laser 
frequency.  The  new  propagation  direction  is  calculated  using  Snell's 
Law. 

After  the  plasma  has  been  traversed,  the  ray  bundle  Is  summed  and 
compared  to  the  Incident  ray  sum  to  determine  net  power  transmission, 
power  being  proportional  to  Intensity.  Each  discreet  temperature  cell 
is  used  to  calculate  total  hemispherical  continuum  radiative  power  lost 
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as  derived  In  [591.  The  power  lost  is  modified  by  a  factor  of  2  to 
include  line  radiation  at  temperatures  below  8,000  K.  For  temperature 
above  8,000  K,  the  line  radiation  factor  as  developed  by  Kozlov  [651  had 
been  used.  The  final  output  file  Is  formatted  for  Input  to  one  of  the 
various  plotting  routines. 
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APPENDIX  C:  NUMERICAL  MODEL  CODE  LISTING 
The  following  is  a  listing  of  the  numerical  model  code  for  plasma 
formation  above  a  metallic  target.  The  details  of  the  code  were 
discussed  In  Chap.  5.  The  code  is  operational  on  the  High  Power  Laser 
Lab  HP  9000  computer  system. 
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